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1.0  statement  of  Problem 


1.1  Introduction 


The  ocean  floor  is  one  of  the  last  unexplored  frontiers 
on  earth.  In  recent  years  a  tremendous  interest  has  been  gen¬ 
erated  for  the  discovery  and  exploration  of  this  unknown  por¬ 
tion  of  the  globe  and  its  exploitation,  both  for  military  as 
well  as  for  civilian  purposes. 

Naval  strategists  are  now  openly  considering  the  feasi¬ 
bility  of  deep,  fixed,  underwater  installations  —  manned  and. 
unmanned,  fixed,  underwater  monitoring  and  surveillance  sta¬ 
tions  for  naval  traffic,  marine  benchmarks,  and  navigational 
aids.  In  the  civilian  sector,  the  last  decade  has  seen  an 
almost  unbelievable  expansion  of  ocean  floor  activities.  Off¬ 
shore  drilling,  mostly  for  natural  oil  and  gas,  is  carried  on 
now  in  water  depths  of  up  to  1000  feet.  The  feasibility  of 
extracting  other  raw  materials  from  the  bottom  of  the  sea, 
such  as  sulfur  and  nodules  of  manganese,  is  seriously  being 
studied.  Plans  are  being  made  to  build  manned  stations  off 
the  east  coast  of  Florida  in  water  depths  to  1000  feet  by  mid- 
1972  and  in  6000-foot  water  depths  by  1973. 

One  fact  is  clear:  the  activities  of  man  on  the  ocean 
floor  will  increase  quickly  and  enormously. 

Even  for  present  activities  and  current  problems  our 
knowledge  of  the  ocean  floor  and  its  behavior  and  response  to 
human  activity  are  woefully  inadequate.  We  are  referring  here, 
e.g.,  to  such  problems  as  the  disappearance  of  the  U.S.  sub¬ 
marine  Thresher,  the  extremely  difficult  and  delicate  recovery 
of  a  thermonuclear  bomb  off  Palomares,  Spain,  and  the  failure 
of  several  radar  platforms  off  the  Atlantic  coast. 

It  is  therefore  singularly  appropriate  that  the  soil 
mechanics  aspects  and  behavior  of  ocean  bottom  sediments  be 
explored  and  investigated.  Thus,  a  more  reliable  assessment 
and  prediction  of  the  response  of  ocean  floor  sediments  to 
human  activities  will  be  possible  and  a  better  estimate  of  the 
performance  of  new  installations,  methods,  and  processes  on 
the  ocean  floor  will  result. 


1 . 2  Statement  of  Problem  and  Objective 


Many  problems  of  ocean  floor  technology  either  are 
directly  concerned  with  or  are  reducible  to  the  problem  of 
penetration  of  an  object  through  the  upper  layers  of  the  bot¬ 
tom  sediments.  This,  for  example,  applies  to  the  installation 
of  underwater  foundations,  the  driving  of  piles  for  over-  and 


underwater  platforms  and  for  dolphins,  the  penetration  of 
anchors  for  ships,  buoys,  mines,  and  monitoring  stations,  the 
discovery  and  reclaiming  of  oojects  lost  and  sunk  to  the  ocean 
floor,  the  assessment  of  ocean  sediment  properties  from  the 
dynamic  penetration  record  of  ocean  floor  penetrometers,  and 
even  the  running  aground  and  refloating  of  ships  and  submar¬ 
ines  . 


Thus,  the  topic  of  this  study  is  extremely  broad.  Our 
objective  is  to  perform  a  comprehensive  review  of  the  state  of 
the  art  of  assessing  and  predicting  the  penetration  of  objects 
into  the  ocean  floor,  to  recommend  methods  and  procedures  for 
immediate  use  for  predicting  such  penetration,  and  to  recom¬ 
mend  future  research  and  development  that  would  significantly 
advance  the  state  of  the  art. 

The  rather  broad  subject  of  penetration  of  an  object  into 
the  ocean  bottom  has  been  subdivided  into  several  categories 
according  to  the  phases  or  periods  during  which  driving  forces 
are  acting  on  the  penetrating  missile.  We  are  distinguishing 
between  these  categories  as  follows. 

For  a  free  falling  object,  the  penetration  phenomenon  is 
characterized  by  the  fact  that  some  time  before  as  well  as 
during  the  penetration  process  no  driving  forces  other  than 
gravitational  forces  are  acting  on  the  object.  Therefore,  in 
free  fall  penetration  (ffp)  the  object  follows  a  free  fall 
trajectory  through  the  water  prior  to  entering  the  ocean  bot¬ 
tom  sediments.  Typical  examples  for  ffp  would  be  a  rock,  a 
bomb,  a  spent  torpedo,  or  a  sunken  ship  falling  to  the  ocean 
bottom.  The  physically  significant  feature  here  is  the  fact 
that,  at  some  time,  say  at  tQ,  travel  under  power  has  stopped 
and  in  the  interval  between  tg  and  the  impact  time  t^  the 
missile  moves  along  a  free  fall  trajectory  which  is  determined 
by  the  velocity  Vq  at  time  to  and  the  frictional  resistance  of 
the  medium  through  which  the  missile  travels. 

On  the  other  hand,  a  missile  may  still  be  under  power 
while  making  impact  with  the  ocean  bottom.  This  type  of  pene¬ 
tration  we  have  chosen  to  call  powered  travel  penetration  (ptp) 
and  typical  examples  would  be  a  torpedo  driving  itself  into 
the  ocean  bottom,  a  submarine  or  ship  running  aground,  or  a 
surface-to-surface  or  air-to-surface  missile  impacting  on  the 
ocean  bottom  while  still  in  the  powered  phase  of  its  flight. 

Finally,  there  are  penetration  phenomena  which  are  char¬ 
acterized  by  the  fact  that  the  penetration  of  the  objects  is 
assisted  by  forces  or  power  systems  designed  to  facilitate 
the  piercing  of  the  ocean  bottom  soils.  We  have  chosen  to 
label  this  type  dynamic  penetration  (dp).  Typical  examples 
for  this  type  of  penetration  would  be  the  driving  of  piles, 
the  penetration  of  anchors  and  the  piercing  of  soil  samplers 
by  explosive  charges,  and  the  drilling  of  holes  into  the 
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ocean  bottom.  Thus,  to  summarize,  the  penetration  of  objects 
into  the  ocean  bottom  can  be  subdivided  into  three  categories: 

a)  free  fall  penetration  ffp 

b)  powered  travel  penetration  ptp 

c)  dynamic  penetration  dp 

To  a  large  measure,  the  penetration  of  an  object  into  the 
ocean  bottom  will  depend  upon  the  properties  of  the  ocean  bot¬ 
tom  sediments  in  addition  to  the  mechanical  properties  of 
the  penetrating  object  and  the  dynamical  parameters  of  the 
penetration  process.  It  is  for  this  reason  that  in  Chapter  2 
of  this  report  we  first  review  our  state  of  knowledge  on  the 
nature,  extent  and  distribution  of  ocean  bottom  sediments,  the 
methods  by  which  knowledge  about  these  factors  is  obtained, 
the  various  methods  of  sampling  and  coring  of  ocean  bottom  sed¬ 
iments,  the  methods  used  for  establishing  the  properties  of 
these  sediments  and,  finally,  the  actual  properties  and  the 
range  within  v/hich  they  might  be  expected  to  fall.  It  is  quite 
clear  that  much  less  is  known  about  submarine  soils  than  about 
terrestrial  soils.  Also,  submarine  soils  appear  to  exhibit  ? 
much  variety  and  variability  as  do  terrestrial  soils.  Thus, 
it  appears  that  the  correct  assessment  of  the  properties  of  the 
bottom  sediment  will  be  one  of  the  important  factors  --  if  not 
the  most  critical  one  --  for  the  successful  prediction  of  ocean 
bottom  penetration.  For  these  reasons  we  carefully  reviewed 
the  literature  regarding  ocean  bottom  sediments.  More  than 
250  reports,  papers,  and  articles  were  reviewed  and  of  these 
about  60  were  found  to  merit  closer  study.  The  most  specific 
treatment  in  the  literature  on  ocean  bottom  penetration  is 
that  discussing  the  penetration  of  core  samplers.  However, 
this  is  by  no  means  an  exhaustive  coverage  of  the  topic  and 
the  penetration  problem  has  to  be  studied  on  a  much  broader 
scale,  taking  into  consideration  penetration  of  an  object  into 
any  kind  of  plastic  material  and  applying  the  results  to  the 
problems  at  hand  by  using  the  properties  of  ocean  bottom  sed¬ 
iments  . 


The  most  important  dynamic  feature  in  free  fall  penetra¬ 
tion  is  the  impact  velocity.  In  addition,  the  travel  path  and 
the  orientation  of  the  object  at  impact  are  important.  For 
this  reason,  Chapter  3  presents  the  methods  that  can  be  used 
for  determining  the  impact  velocity  and  trajectory  of  an 
object  falling  through  sea  water. 

In  Chapter  4  we  review  the  mechanics  of  penetration  of 
an  object  striking  a  target  material  with  an  impact  velocity 
Uq  .  The  problem  is  subdivided  into  two  classes.  The  first 
class  is  penetration  with  increasing  penetration  area.  This 
usually  will  apply  most  often  to  large  object  penetrating  only 
partly  into  the  bottom  soil  (shallow  penetration).  The  second 
class  is  penetration  at  constant  area.  While  shallow  penetra¬ 
tion  may  occur  at  constant  area,  deep  penetration  practically 
always  will  be  at  constant  area. 


4 


The  problem  of  dynamic,  driven  penetration  is  discussed 
in  Chapter  5.  For  most  practical  purposes  dynamic  penetration 
is  the  driving  of  piles  or  other  objects  by  a  series  of  imp.-'’.ct 
blows  into  the  bottom  soils  and  the  mechanics  of  this  process 
is  reviewed. 

In  Chapter  6  we  I’eview  the  state  of  our  technology  for 
achieving,  ocean  bottom  penetration  and  report  on  coring  tools, 
anchors  and  pile  driving  technology  presently  existing  and 
■available  for  ocean  bottom  penetration. 

A  rather  comprehensive  and  competent  offshore  construc¬ 
tion  technology  has  developed  during  the  past  two  decades 
mainly  in  connection  with  offshore  exploration  and  production 
of  oil.  All  fixed  platforms  are  anchored  to  the  ocean  bottom 
by  piles  driven  through  template  sleeves.  The  development  of 
this  technology  has  been  traced  through  the  litei^ature  and 
the  present  state  of  the  art  was  assessed. 

The  piles  anchoring  the  platforms  are  usually  driven 
above  water  and,  after  they  have  encountered  refusal  or  the 
required  driving  resistance,  they  are  anchored  to  the  plat¬ 
form  by  welding  and/or  grouting  of  the  space  between  the  pipe 
sleeve  and  the  pile. 

As  of  this  time,  no  record  exists  that  piles  for  an  off¬ 
shore  platform  have  ever  been  driven  under  water  although  at 
least  one  offshore  contractor  has  a  patent  application  pending 
for  an  underwater  pile  hammer. 

It  appears  that  the  most  reliably  successful  method  of 
installing  piles  under  water  and  at  greater  water  depths  (say 
in  excess  of  300-400  feet)  will  be  by  pile  vibrators. 

In  Chapter  7  we  present  recommendations  for  immediate 
use  of  methods  to  assess  and  predict  ocean  bottom  penetration. 
Finally,  we  recommend  future  research  and  development  to 
advance  the  state  of  the  art  of  achieving  ocean  bottom  pene¬ 
tration  as  well  as  predicting  it  (Chapter  8). 
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2 . 0  Ocean  Bottom  Soils  and  Their  Properties 


2 . 1  Ocean  Bottom  Topography 

The  principal  topographic  features  that  characterize  the 
ocean  bottom  are  the  continental  shelf  or  terrace  and  the 
abyssal  plain  of  the  deep  ocean  floor.  The  abyssal  plains 
cover  millions  of  square  miles  and  are  covered  with  deep  sedi¬ 
ment.  The  continental  terrace  is  the  gradual  extension  of  a 
continent  into  the  ocean  down  to  the  continental  slopes  Usu¬ 
ally  the  end  of  the  continental  terrace  and  beginning  of  the 
continental  slope  is  taken  at  a  water  depth  of  600  feet. 

Beyond  that  depth  the  ocean  bottom  usually  drops  off  steeply 
to  the  abyssal  deep.  The  abyssal  deep  is  considered  to  begin 
at  a  depth  from  6,000  to  12,000  feet.  These  main  features  of 
the  ocean  bottom  are  accompanied  by  undersea  mountain  ranges, 
especially  the  midocean  ridges,  isolated,  submerged  mountains 
and  sea  mounts  or  guyots  rising  thousands  of  feet  above  the 
bottom  and  escarpments  a  mile  or  more  in  height  cutting  across 
the  ocean  floor  for  hundreds  of  miles.  In  addition,  the  ocean 
floor  has  been  warped  by  folding  and  faulting,  causing  deep 
submerged  trenches  up  to  20,000  feet  deep  and  resulting  in 
total  water  depths  of  nearly  36,000  feet.  Cutting  into  the 
continental  shelves  are  deep  submarine  canyons  located  usually 
where  surface  rivers  are  discharging  into  the  oceans . 

Thus,  the  ocean  bottom  geology  appears  to  be  almost  as 
varied  as  that  of  the  continents.  The  general  bathymetry  of 
the  ocean  basins  is  provided  by  standard  reference  texts  such 
as  Menard  ( 1964 )  . 


2.2  Ocean  Bottom  Sediments 


The  soils  of  the  ocean  bottom  are  usually,  classified  into 
two  big  groups:  the  terrigenic  soils  and  the  pelagic  soils. 
The  terrigenic  soils  as  their  name  suggests  are  derived  from 
the  continents  and  are  carried  into  the  oceans  by  rivers  and 
currents.  They  are  the  sediments  found  on  the  continental 
shelf  and  are  made  up  of  sands,  silts,  muds  and  clays.  Sub¬ 
marine  sediment  transport  is  controlled  by  the  same  physical 
laws  as  sediment  transport  on  the  continent  and,  hence,  the 
sediments  on  the  continental  shelves  have  similar  particle 
size,  shape  and  characteristics  as  soils  on  land  with  one 
important  exception:  since  they  have  been  deposited  in  deep, 
salt  water,  the  clay  particles  are  flocculated  and,  therefore, 
have  a  very  loose  structure  and,  correspondingly,  very  low 
densities  and  high  porosities  or  void  ratios. 

Pelagic  sediments  blanket  the  abyssal  deep.  These  deep 
sea  sediments  are  clays  and  fine  silts  ground  up  by  the  surf 
and  carried  from  the  shore  line  out  into  the  ocean;  they  also 
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contain  wind-blown  volcanic  and  cosmic  dust  and  cirganic  or 
authigenic  materials  foi’med  in  the  deep  ocean  environment. 
Sediments  containing  more  than  30  per  cent  of  material  of 
organic  origin  in  the  form  of  tests  and  frustules  (or  skele¬ 
tons)  of  microscopic  animals  and  plants  are  called  oozes. 

The  flora  and  fauna  whose  remnants  make  up  the  oozes  may  be 
either  siliceous  or  calcareous.  The  oozes  are  classified  as 
globigerina  ooze,  radiolarian  ooze,  pteropod  ooze,  or  diatom 
ooze,  depending  upon  the  predominant  type  of  animal  or  plant 
from  which  the  organic  portion  of  the  sediment  is  derived. 

Red  clay  is  a  reddish-to-brown  deposit  which  may  be  composed 
partially  of  terrigenous  colloidal  matter,  the  insoluble  por¬ 
tion  of  organism  tests  and  the  products  of  submarine  volcan- 
ism . 

The  floor  of  the  Pacific  Ocean  is  covered  chiefly  by 
these  clastic,  red  clay  deposits.  The  predominant  sediments 
in  the  Atlantic  and  Indian  Oceans  are  calcareous  oozes;  sili¬ 
ceous  oozes  are  less  widely  distributed. 

The  prevalence  of  these  sediments  is  shown  in  Table  2.1. 


Table  2.1  Percentage  of  Sea  Floor  Covered  by  Sediments 


Type  of  Deposit 

Percentage 

of 

Sea  Floor 

Average 

Depth 

(in  meters) 

Shelf  sediments  (terrigenic) 

8.3 

100 

Muds  ( hemipelagic ) 

18 

<2300 

Pelagic  Sediments: 

Globigerina  ooze  (calcareous) 

35.5 

3600 

Pteropod  ooze  (calcareous) 

0.6 

2000 

Diatom  ooze  (siliceous) 

8.5 

3900 

Radiolarian  ooze  (siliceous) 

1.9 

5300 

Red  Clay 

28.3 

5400 

While  the  descriptive  names  suggest 
tion  of  the  various  sediments,  in  reality 
the  ocean  floor  grade  into  one  another  an 
arbitrary  whether  a  specific  sample  is  cl 
an  ooze  and,  if  the  latter,  what  type  of 
the  population  of  the  organic  constituent 
one  although  sometimes  with  one  skeletal 
dominating . 


a  cl 

ear- 

cut  separa- 

the 

sed 

iments  on 

d  it 

is 

somewhat 

assi 

fied 

a  clay  or 

ooze 

sin 

ce ,  again , 

s  is 

usu 

ally  a  mixed 

spec 

ies 

clearly  pre- 

The  type  of  pelagic  sediment  depends  very  much  on  the 
depth,  the  pressure,  and  the  temperature  of  the  water  at  the 
ocean  bottom.  This  is  so  because  the  solubility  of  calcium 
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carbonate  is  a  function  of  pressure,  temperature,  and  carbon 
dioxide  content  of  the  water.  Thus,  calcareous  oozes  are 
deposited  in  more  shallow  basins  where  the  water  temperatures 
are  higher  and  the  pressures  lower.  In  colder  waters,  diatoms 
with  siliceous  skeletons  predominate  over  the  calcareous  fora- 
minifera.  The  red  clays  are  found  in  the  coldest  and  deepest 
parts  of  the  oceans. 

The  rates  of  deposition  of  recent  pelagic  sediments  have 
been  calculated  and  vary  from  0.01  to  6  centimeters  per  1000 
years.  For  red  clays,  measurements  gave  a  deposition  rate  of 
0.07  to  0.2  centimeters  per  1000  years;  for  globigerina  ooze 
in  the  equatorial  Atlantic  Ocean,  about  1.6  centimeters  per 
1000  years  were  found.  These  rates  of  sedimentation  depend 
not  only  on  the  geographic  location  but  also  on  the  time  of 
deposition  because  of  the  variation  of  climate  over  geologic 
time.  According  to  Hamilton  (1959)  the  observed  rates  of  sed¬ 
imentation  as  well  as  geochemical  and  other  calculations  of 
the  total  volume  of  sediment  that  should  be  present  in  deep 
sea  basins,  based  on  materials  released  in  weathering,  indi¬ 
cate  that  the  average  thickness  of  sediment  solids  on  the 
ocean  floor  should  be  1  to  3  kilometers.  Seismic  surveys  in 
all  oceans,  on  the  contrary,  report  a  sediment  thickness  of 
0.1  to  0.5  kilometers,  but  also  reveal  that  in  many  areas  the 
surface  sediments  are  underlain  by  other  layers  whose  elastic 
wave  transmission  velocities  are  intermediate  between  those 
of  sediment  and  the  "basaltic  basement"  rock  of  the  oceanic 
crust.  This  observation  has  led  Hamilton  (1959)  to  the  con¬ 
clusion  that  the  lower  layers  are  composed  of  rocks  formed  by 
the  normal  consolidation  and  lithification  of  the  present 
types  of  deep  sea  sediments. 

Rock  outcrops  occur  over  only  a  few  per  cent  of  the  ocean 
floor.  They  are  found  usually  in  many  isolated  areas,  par¬ 
ticularly  on  the  continental  shelves  and  slopes,  around  islands, 
banks,  seamounts,  submarine  ridges  and  mountain  ranges. 


Properties  of  Ocean  Bottom  Sediments 


In  recent  years  an  increasing  number  of  scientists  have 
explored  the  properties  of  ocean  bottom  sediments.  At  present, 
direct  observation  of  deep  ocean  sediments  is  limited  to  such 
samples  as  may  be  recovered  by  dredging  or  coring  operations. 
These  methods  permit  sampling  of  the  sediment  to  only  shallow 
depths  below  the  sediment  surface.  Hopkins  (1964)  has  prepared 
a  comprehensive  survey  of  marine  bottom  samplers,  and  Richards 
and  Parker  (1967)  discuss  the  problems  and  methods  of  coring. 
Until  recently  drilling  has  been  carried  out  only  in  the  shal¬ 
low  waters  of  the  continental  shelves. 


At  present,  however,  the  National  Science  Foundation  is 
sponsoring  a  program  to  drill  deep  ocean  cores  at  some  thirty 
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locations  in  the  Atlantic  and  thirty-six  locations  in  the 
Pacific.  This  program  promises  to  add  significant  information 
about  the  deeper  layers  of  the  deep  sea  sediments. 

The  topmost  parts  of  the  sedimentary  column  that  could  be 
sampled  earlier  have  been  studied  and  analyzed  in  detail.  Phy 
ical  properties  such  as  water  content,  porosity,  void  ratio, 
density,  particle  size  and  shape,  mineral  composition,  shear 
strength,  consolidation  characteristics,  elastic  wave  veloci¬ 
ties,  attenuation,  etc.  can  and  have  been  determined. 

Acoustic  measurements  play  a  major  role  in  the  study  of 
physical  properties  of  marine  sediments.  This  is  so  because 
acoustic-seismic  data  are  very  easily  obtainable  and  are  plen¬ 
tiful.  Nafe  and  Drake  (1963)  and  Horn  et  al .  (1968)  attempt 
to  correlate  wave  propagation  velocities  with  such  elementary 
properties  as  porosity,  density  and  shear  strength. 

The  results  presented  by  Nafe  and  Drake  show  a  fairly 
good  correlation  between  compression  wave  velocity  and  poros¬ 
ity  with  but  a  moderate  scatter.  However,  their  data,  unfor¬ 
tunately,  show  good  discrimination  only  for  sedimentary  rocks 
at  fairly  low  porosities.  For  unconsolidated  ocean  sediments 
measured  velocities  vary  in  a  very  narrow  range,  largely  from 
1.5  to  1.8  km/sec.  while  the  porosities  range  from  35  to  85 
per  cent.  However,  the  scatter  is  wide  enough  to  make  a  cor¬ 
relation  questionable.  The  correlation  between  wave  velocity 
and  density  was  found  to  be  similar.  Nafe  and  Drake  suggest 
that  much  more  information  may  become  available  in  the  near 
future  from  machine  computations  for  dispersion  analysis  where 
computed  phase  and  group  velocities  for  assumed  sediment 
structures  and  properties  may  be  compared  to  actually  observed 
values . 

Horn  et  al.  attempt  a  correlation  between  acoustical 
properties  and  other  physical  properties  of  sediment  cores. 
They  plot  compressional  wave  velocities  against  porosity, 
moisture  content,  void  ratio,  dry  density  and  wet  density. 

The  maximum  sound  velocities  were  found  in  a  layer  of  silt 
from  the  Norwegian  basin  with  a  porosity  of  35  per  cent  (1732 
m/sec.)  and  in  a  layer  of  fine  grained  sand  from  the  Mediter¬ 
ranean  with  a  porosity  of  34  per  cent  (1832  m/sec.).  The  low¬ 
est  velocities  were  observed  in  clays  with  porosities  varying 
from  70  to  84  per  cent  (1455  to  1480  m/sec.).  However,  again, 
the  scatter  of  the  data  is  wide.  The  conclusions  of  Horn  et 
al.,  perhaps,  deserve  to  be  cited  here:  "The  results  of  this 
study  reveal  that  certain  bulk  and  textural  properties  of 
unconsolidated  marine  sediments  are  indices  of  their  acousti¬ 
cal  characteristics.  Of  the  bulk  properties,  a  decrease  in 
porosity,  moisture  content,  and  void  ratio  is  matched  by  an 
increase  in  sound  velocity.  There  is  no  correlation  between 
sound  velocity  and  either  dry  density  or  carbonate  content. 
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"When  shear  strength  is  plotted  against  sound  velocity, 
the  data  fall  into  distinct  groups  which,  in  general,  reflect* 
mean  grain  size  and  sediment  genesis.  There  is  an  overall 
increase  in  sound  velocity  as  the  shear  strengths  get  larger. 
From  the  limited  amount  of  data  on  hand  concerning  ash  layers, 
these  sediments  have  constant  sound  velocities  regardless  of 
their  shear  strength. 

"The  shear  strength  is  not  considered  to  be  a  reliable 
index  of  the  acoustical  properties  of  unconsolidated  deep  sea 
sediments . 

"A  definite  relationship  exists  between  sound  velocity 
and  mean  grain  size.  The  coarser  a  sediment  is,  the  higher 
is  its  sound  velocity...." 

Faas  (1968)  recently  reported  a  relationship  between 
the  reflection  coefficients  and  the  porosity  of  ocean  bottom 
sediments  and  claims  that  it  can  be  used  for  predictive  pur¬ 
poses.  In  addition,  he  observes  that  areas  of  sand  could  be* 
classified  as  acoustically  "hard,"  whereas  areas  of  clay,, 
from  which  reflections  were  greatly  subdued,  were  c.l.3ss-i#ied 
as  acoustically  "soft."  At  any  rate,  it  appears  that  acoustic 
measurements  are  capable  of  providing  much  more  information 
than  currently  is  obtained  from  their  records. 

Extensive  investigations  on  mass  physical  properties- of 
ocean  bottom  sediments  were  carried  out  by  Richards  (1961, 
1962),  Moore  (1961,  1962,  1964),  Harrison  et  al.  (1964)  and 
Hamilton  (1956,  1964).  Tabid  2.2  shows  the  data  collected  by 
Moore  on  more  than  80  core  samples  from  the  North  Pacific. 

A  comprehensive  summary  of  most  available  data  was  pre-  * 
sented  by  Keller  (1967)  who  gave  the  results  of  tests  on  some 
500  sediment  cores  (Atlantic,  300;  Pacific,  200).  He  pre¬ 
sented  in  particular  the  ranges  of  sediment  types,  shear 
strengths,  water  content,  and  dry  and  wet  densities  that  can 
be  expected  in  both  ocean  basins.  For  the  problem  of  penetra¬ 
tion  the  shear  strength  is  the  most  critical  property.  There¬ 
fore,  Figure  2.1  and  Figure  2.2  reproduce  the  data  presented 
by  Keller  on  shear  strength  in  the  Atlantic  and  Pacific 
respectively.  The  shear  strengths  reported  were  measured  on 
marine  sediment  cores  by  either  a  laboratory  vane  shear  test 
or  an  unconfined  compression  test.  The  cores  used  varied  in 
length  from  12  inches  to  20  feet,  with  an  average  of  about 
seven  feet. 

Average  shear  strengths  found  in  these  cores  range  from 
less  than  0.5  psi  to  2.5  psi  for  the  upper  few  feat  of  sea 
floor  sediments.  These  extremely  low  values  are  caused  by 
the  extremely  low  values  of  dry  density  and  the  corresponding 
high  values  for  porosity  and  water  content  of  the  sediments. 

An  inspection  of  Figures  2.1  and  2.2  shows  that  sediments 
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from  the  North  Pacific  (after  Moore). 
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with  a  shear  strength  from  0.5  to  1.0  psi  predominate  in  the 
North  Atlantic.  The  highest  strengths  found  in  the  basin 
range  from  1.0  to  1.5  psi  and  are  associated  with  calcareous 
deposits.  Values  of  less  than  0.5  psi  are  often  found  in 
coastal  areas  where  local  conditions  of  drainage  or  current 
strongly  influence  the  depositional  environment.  Another 
area  of  low  shear  strength  is  the  red  clay  -^deposit.  The  red 
clay  deposits  predominate  and  cover  a  lerge  portion  of  the 
area  in  the  North  Pacific  basin.  It  is  for  this  reason  that, 
in  general,  the  sediments  in  that  basin  are  weaker  than  those 
in  the  North  Atlantic.  As  far  as  general  deviations  from 
these  patterns  go,  Keller  writes: 

"In  contrast  to  the  North  Atlantic,  large  portions  of 
the  North  Pacific  sea  floor  are  covered  with  sediments  whose 
average  shear  strength  is  less  than  0.5  psi.  These  areas 
coincide  closely  with  the  distribution  of  *red  clay'  and  com¬ 
prise  a  major  portion  of  the  sea  floor.  Localities  of  higher 
shear  strength  occur  within  the  area  of  low  shear  strength  as 
a  result  of  changes  in  bottom  topography  influencing  the 
depositional  environment.  Local  currents  in  and  around  topo¬ 
graphic  features  can  account  for  changes  in  the  distribution 
of  certain  sediment  properties.  It  has  been  found  that  on 
topographic  'highs'  shear  strengths  are  slightly  higher  than 
in  the  surrounding  areas.  This  may  be  attributed  to  the 
winnowing  effect  of  currents  which  tend  to  keep  the  'highs' 
free  of  softer  sediments. 

"Shear  strength  values  are  generally  higher  along  the 
margins  of  the  basin  and  in  the  lower  latitudes.  Coarser 
sediments  and  shallower  water  depths,  normally  found  closer 
to  land,  account  for  the  general  pattern  shown  in  Figure  2.2. 
In  areas  of  increased  calcium  carbonate,  such  as  the  low  lat¬ 
itudes,  shear  strength  is  also  found  to  increase.  The  high¬ 
est  range  of  values  2.0  to  2.5  psi  observed  thus  far  occur  in 
the  calcareous  oozes  of  the  Pacific  basin  (Figure  2.2). 

"A  comparison  of  the  overall  data  presented  in  Figures 
2.1  and  2.2  indicates  that  North  Atlantic  sediments  possess 
relatively  higher  shear  strengths  than  do  those  of  the  North 
Pacific.  It  is  also  evident  from  these  data,  that  the  North 
Pacific  basin  can  be  divided  into  two  sedimentary  provinces, 
each  distinctly  different  from  the  other.  The  northern  por¬ 
tion  consists  of  sediment  possessing  strengths  ranging  from 
0.25  to  0.5  psi;  in  the  lower  latitudes,  values  of  1.0  to  1.5 
predominate . . . . " 

In  addition  to  the  shear  strength,  the  porosity  and  the 
water  content  of  the  ocean  bottom  sediments  also  appear  to 
have  some  albeit  reduced  and  indirect  influence  on  the  pene¬ 
tration  problem.  The  former  because  of  its  measure  of  solid 
content  and  hence,  its  point-to-point  contacts  and  bonds  in 
the  sediment,  the  latter  because  of  the  pore  water  pressures 
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that  result  in  any  straining  or  deformation  process  in  a  fine 
grained  sediment.  In  general,  the  water  content  of  ocean  bot¬ 
tom  samples  ranges  from  40  per  cent  to  200  per  cent  with  the 
majority  of  samples  in  the  range  of  80  to  150  per  cent;  the 
porosity  ranges  from  40  to  85  per  cent  with  the  majority  of 
the  samples  in  the  range  of  60  per  cent  to  80  per  cent. 


2 . 4  Variation  of  Sediment  Strength  with  Depth 

Most  seismic  and  acoustic  data  on  wave  propagation  indi¬ 
cates  that  wave  velocities  increase  with  depth  in  the  sediment. 
Also,  since  the  overburden  pressure  in  the  sediment  increases 
linearly  with  depth,  one  would  assume  that  gradual  consolida¬ 
tion  takes  place  in  the  sediment  leading  to  a  gradual  decrease 
in  porosity,  water  content  and  void  ratio  and  a  corresponding 
increase  in  shear  strength  with  depth.  Actual  measurements  of 
the  variation  of  these  properties  with  depth  show  rather  inter¬ 
esting  and  unexpected  results  (Hamilton,  1964;  Moore,  1964). 

Contrary  to  what  one  would  expect,  the  porosities  and 
void  ratios  do  not  consistently  decrease  with  depth.  Except 
for  some  decrease  close  to  the  sediment  surface,  the  porosity 
was  found  to  remain  constant  to  depths  of  up  to  140  meters 
(460  feet).  Since  the  overburden  pressure  clearly  increases 
with  depth,  this  constancy  of  the  porosity  and,  hence,  of  the 
soil  skeleton  is  only  explicable  by  the  fact  that,  simultan¬ 
eously  with  the  pressure,  the  strength  also  increases.  This 
indeed  was  found  to  be  true.  Various  attempts  to  correlate 
sediment  strength  with  other  mass  properties  were  rather  neg¬ 
ative  except  for  an  apparent  correlation  of  increased  shear 
strength  with  relative  percentage  of  montmorillonite  and 
rather  good  agreement  between  carbonate  content  and  strength. 
Moore  suggests  the  good  agreement  between  carbonate  content 
and  strength  to  be  due  to  incipient  cementation  by  crystali- 
zation  of  carbonate  solutions.  Many  investigators  have 
observed  that  slow  deposition  and  great  age  can  contribute  to 
sediment  strength.  This  strength  may  be  attributable  both  to 
chemical  cementation  or  to  rigid  hydrogen  bonds  of  adsorbed 
water  molecules.  Both  have  been  suggested  to  result  from  age 
and  extended  exposure  to  chemical  solutions . 

Figure  2.3  shows  the  variation  of  sediment  strength  with 
depth  as  a  function  of  the  rate  of  deposition.  It  is  quite 
evident  that  the  rate  of  deposition  and,  hence,  the  age  of 
the  sediment  has  a  most  significant  influence  on  the  variation 
of  strength  with  depth.  The  data  show  a  linear  increase  of 
strength  with  depth.  Similarly,  Hironaka  (1966)  reported 
strength  variations  close  to  the  surface  by  performing  vane 
shear  tests  at  5-inch  intervals  in  2  to  3-foot  long  cores 
taken  from  the  sediment  surface  near  San  Miguel  Island.  By  a 
stepwise,  linear  regression  analysis,  he  derived  the  following 
equation  to  fit  his  data: 
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VSS  =  k^(D)  +  k2(LL)  +  kgCMD)  +  k4  (2.1) 

where 

VSS  =  vane  shear  strength,  psi 
D  =  depth  in  sediment,  inch 
MD  =  mean  grain  diameter,  mm 
LL  =  liquid  limit,  per  cent 
kjj  =  constants  (n  =  1,  2,  3,  4) 

The  first  term  of  equation  (2.1)  again  suggests  a- linear 
increase  of  strength  with  depth.  Although  Hironaka's  analyses 
were  based  on  laboratory  measurements,  he  believes  that' these 
relationships  are  also  valid  for  in  situ  strength  values. 

Data  reported  by  Bayles  (1965)  from  cores  taken  at  the  site  of 
the  Thresher  disaster  also  reveal  a  linear  increase  of  the 
shear  strength  with  depth  from  a  value  of  0.3  psi  at  the  sed¬ 
iment  surface  to  1.7  psi  at  a  depth  of  200  feet. 


2 . 5  Environmental  Effects  on  Properties  of  Ocean  Sediments 

Most  of  the  data  available  on  ocean  bottom  sediment 
strength  are  from  laboratory  tests  by  vane-shear,  direct-shear 
box  or  triaxial  compression  tests. 

This  means  that  the  test  specimens  were  removed  from- 
their  in  situ  environment  and  were  tested  for  strength,  and'  • 
other  properties  in  the  laboratory.  It  is  evident  that  .samp¬ 
ling  and  removing  from  the  ocean  bottom  environment  exposes 
the  test  samples  to  an  entirely  different  regime  than  that 
experienced  by  the  ocean  bottom  soil  in  situ. 

The  question  is  if  and  how  much  such  a  change  in  the 
environment  will  affect  the  properties,  particularly  the  shear 
strength.  Crisp  (1968)  presents  an  analysis  according  to 
which  the  strength  should  decrease  by  10  to  15  per  cent  due 
to  the  change  in  particle  spacing  from  the  expansion  of  the 
pore  water  when  a  sample  is  brought  from  a  water  depth  of 
12,000  feet  to  the  surface.  However,  in  his  calculations  he 
does  not  consider  the  effect  of  the  shear  distortion  that  is 
caused  in  a  sample.  Since  the  radial  expansion  in  the  core 
sample  is  effectively  prevented  by  the  sampler,  all  expansion 
necessarily  must  take  place  in  the  axial  direction. 

Crisp  also  analyzes  the  effect  of  dissolved  gases  (oxygen 
and  carbondioxide )  and  concludes  that  it  is  negligible.  Vey 
and  Nelson  (1967)  performed  direct  shear  tests  and  vane  shear 
tests  in  a  seawater  environment  with  a  hydrostatic  environ¬ 
mental  pressure  chamber  capable  of  providing  up  to  10,000  psi 
pressure  in  a  temperature  environment  of  1  to  3°C.  The 


results  are  somewhat  inconclusive  due  to  a  wide  scatter  of  the 
data,  but  their  tests  indicate  that  the  direct  shear  tests 
showed  a  decrease  in  shear  strength  with  increasing  environ¬ 
mental  pressures,  while  the  vane  shear  test  showed  an  increase 
in  shear  strength  with  environmental  pressures  (as  is  suggested 
by  Crisp’s  analysis)  for  the  more  plastic  soil  but  a  decrease 
for  the  less  plastic  soil  tested.  Whatever  the  effect,  in  this 
writer's  opinion^ the  data  reported  by  Vey  and  Nelson  have  such 
a  wide  scatter  that  it  is  difficult  to  make  any  conclusions 
based  on  these  experiments.  Certainly,  variations  in  strength 
of  10  to  15  per  cent  would  not  be  revealed  by  the  amount  and 
range  of  data  reported  by  them.  This  should  not  be  considered 
as  a  criticism  of  their  work  .  The  problem  of  making  shear 
strength  measurements  at  the.  high  void  ratios  and  the  low 
strength  properties  of  ocean  bottom  sediments  is  an  enormously 
difficult  one. 

Another  factor  that  enters  the  discussion  here  is  the 
problem  of  sample  disturbance .  In  addition  to  the  stress 
release,  the  sampling  procedure  may  also  cause  considerable 
disturbance  or  remolding.  The  problems  and  requirements  for 
proper  coring  and  minimum  disturbance  were  discussed  by 
Richards  and  Parker  (1967). 

One  might  conclude  that  the  extremely  low,  measured  values 
of  the  shear  strength  of  ocean  bottom  sediments  must  be  caused 
by  a  combination  of  stress  release,  pore  water  expansion,  and 
sample  disturbance,  and  that  in  situ  strength  is  considerably 
higher.  However,  the  limited  data  available  on  in  situ  obser¬ 
vations  (Hironaka  and  Smith,  1967)  as  well  as  the  fairly  high 
porosities  of  the  sediments  which,  really,  are  not  in  doubt 
preclude  a  gross  divergence  from  the  measurements  made  on  cores. 
Undoubtedly,  all  the  factors  cited  above  introduce  errors  into 
the  measurement  of  shear  strength  and  in  situ  strength  should 
be  higher  than  the  one  obtained  from  cores.  However,  a  drastic 
increase  such  as  one  order  of  magnitude  cannot  be  expected,  but 
rather  an  increase  of  perhaps  30  to  50  or,  at  most,  100  per 
cent . 

The  best  solution  for  overcoming  this  difficulty  would  be 
a  relatively  simple  device  capable  of  measuring  the  in  situ 
strength  of  a  sediment  quickly,  reliably,  and  cheaply.  The 
most  promising  approach  for  developing  such  an  instrument  is 
believed  to  be  a  penetrometer-type  instrument  mounted  with  an 
accelerometer  as  developed  by  Scott  (1967)  or  Schmid  (1966). 

Data  from  such  instruments  could  be  collected  easily  and 
quickly.  It  would  have  the  additional  advantage  that  the 
results  are  immediately  available  and  the  strength  test  is 
performed  in  the  natural  environment  of  the  ocean  bottom  soil. 
Impact  deceleration  data  could  be  transformed  into  shear 
strength  data  by  relations  such  as  given  by  equations  (4.21), 
(4.22),  (4.23),  (4.28)  or  (4.36)  of  Chapter  4. 


2 . 6  Summary 


While  we  do  know,  in  a  general  way,  the  nature  and  char¬ 
acteristics  of  ocean  bottom  sediments,  our  knowledge  in  detail 
and  at  specific  locations  is  extremely  scant  or  nil.  For 
example,  the  information  presented  by  Keller  (1967)  is  based 
on  the  results  of  some  300  sediment  cores  in  the  North  Atlanti 
This  amounts  to  an  average  coverage  of  one  sample  per  30,000 
square  miles.  In  the  Pacific  the  situation  is  even  worse.  It 
is  quite  evident,  therefore,  that  the  extrapolation  of  these 
data  has  been  carried  extremely  far.  While  the  variation  of 
soil  properties  of  the  ocean  bottom  in  a  horizontal  direction 
is  not  likely  to  be  as  great  as  it  is  on  land,  nevertheless 
such  variations  do  exist  and  our  present  knowledge  of  ocean 
bottom  sediments,  their  properties,  and  their  topography  will 
present  a  challenge  to  the  best  minds  and  skills  for  some  time 
to  come. 
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3 .0  Free  Fall  Velocity  and  Trajectory 

As  stated  before,  the  most  important  parameter  determin¬ 
ing  the  penetration  characteristics  is  the  velocity  av  impact. 
In  addition,  the  direction  of  the  impact  as  well  as  the  orien¬ 
tation  of  the  object  at  first  contact  are  controlling  factors. 
Thus,  it  becomes  necessary  to  evaluate  in  detail  the  velocity, 
attitude,  and  trajectory  of  an  object  falling  to  the  ocean 
bottom. 


3 . 1  Hydrodynamic  Drag 

When  an  object  falls  through  a  fluid  at  rest  the  gravita¬ 
tional  pull  tends  to  accelerate  the  object  in  the  direction  of 
the  gravitational  field.  At  the  same  time,  the  fluid  exerts 
a  force  on  the  object  in  a  direction  opposite  to  that  of  the 
motion.  This  force  is  called  the  profile  drag  Fp,  or  "resist¬ 
ance  . " 


Fundamentally,  the  total  force  on  each  element  of  the 
object's  surface  can  be  resolved  into  a  formal  and  a  tangential 
component  as  shown  in  Figure  3.1. 


The  normal  components  are  pressure  forces  and  their  re¬ 
sultant  against  the  direction  of  motion  is  the  pressure  drag 

This  pressure  drag  is  given  by: 


(Fp) 


P 


sin  0  dA 


(3.1) 


The  tangential  components  are  viscous  or  frictional 
resistances  and  their  resultant  against  the  direction  of  the 
motion  is  called  the  friction  drag  (Fj))f. 


The  friction  drag  is  given  by: 


(FD)f 


cos  9  dA 


(3.2) 


The  total  drag,  then, 
Fp  =  (Fj))p  +  (Fd)^  = 


is  the  sum  of  both  components 


sin  6  dA  + 


cos  6  dA 


(3.3) 


The  relative  magnitude  of  the  two  drag  components  depends 
very  much  on  the  shape  and  orientation  of  the  falling  object. 
For  a  streamlined  body  with  its  axis  parallel  to  the  motion, 
the  stream  lines  conform  very  well  to  the  surface  configuration 
on  the  body.  Here,  the  pressure  drag  is  very  small  and  is  of 
the  same  order  of  magnitude  as  the  friction  drag. 
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For  a  plate  placed  perpendicular  to  the  motion  (see  Figure 
3.2),  the  total  drag  is  almost  entirely  due  to  pressure  differ¬ 
ences  between  the  upstream  and  downstream,  sides,  while  the 
friction  drag  is  practically  negligible.  On  the  other  hand,  if 
the  plate  is  placed  parallel  to  the  motion,  the  pressure  drag 
is  practically  zero  and  the  total  drag  is,  therefore,  essen¬ 
tially  equal  to  the  friction  drag. 

In  general,  the  total  resistance  or  drag  is  measured  exper¬ 
imentally  and  is  then  given  in  terms  of  the  density  of  the  fluid, 
the  velocity  of  motion,  the  projected  area  of  the  object  in  the 
direction  of  motion,  and  a  drag  coefficient  Cj^: 


F 


D 


2 


(3.4) 


where  Fq 
Pf 

Yf 

u 

A 


drag  force  (pounds)^  /  o 

.  Yf  /^Ib  sec 

density  of  fluid  =  —  | - JT — 

g  \  ft"* 

unit  weight  of  the  fluid  (Ib/ft^) 

velocity  of  object  (ft/sec) 

projected  area  normal  to  the  direction  of  flow 
or  motion  (ft^) 

drag  coefficient 


3 . 2  Drag  Coefficients 

The  drag  force  against  an  object  depends  on  its  shape,  its 
orientation,  the  Reynolds  number,  and  also,  on  the  nature  of 
the  incident  flow.  Because  of  the  infinite  variety  of  shapes 
it  is  possible  to  consider  only  a  few  typical  forms. 

A)  Smooth  Flat  Plate  of  Length  1  at  Zero  Angle  of  Attack 
in  Zero  Pressure  Gradient.  The  drag  is  due  to  the  shear 
stresses  on  the  plate  and  is  said  to  be  frictional.  If  the 
boundary  layer  is  laminar  from  the  leading  edge  to  a  distance 
X  along  the  plate  and  the  total  drag  is  given  by: 


F 


D 


Pf  u' 


( 3.4a) 


the  drag  coefficient  is  given  by  Blasius  as: 


1.328 


(3.4b) 


where  =  ux/v  is  the  Reynolds  number  and  v  is  the  kinematic 

viscosity  of  the  fluid  (ft^/sec). 
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B)  Streamline  Forms.  Streamline  forms  are  defined  as 
those  along  which  there  is  no  separation  of  flow.  Since  a  form 
on  which  separation  occurs  at  low  Reynolds  number^  or  in  a  low- 
turbulence  stream^ may  become  separation-free  at  higher  Reynolds 
numbers,  or  in  a  stream  of  higher  turbulence,  whether  or  not  a 
form  is  streamlined  depends  upon  these  conditions . 

The  drag  of  a  streamline  form  is  principally  frictional, 
i.e.,  due  to  its  shear  stresses.  The  boundary  layer  is  laminar 
up  to  some  transition  point  x<p  (which  fluctuates  even  under 
steady  conditions  and  can  be  defined  only  in  a  statistical 
sense),  beyond  which  it  becomes  turbulent.  At  Reynolds  numbers 
R  =  uc/v  <  10^,  in  which  c  is  the  chord  of  the  profile,  laminar 
separation  occurs  unless  the  thickness-chord  ratio  d/c  is  less 
than  0.1  (see  Figure  3.4).  The  experimental  results  of  tests 
of  various  streamlined  shapes,  including  elliptical  sections, 
can  be  summarized  in  the  empirical  formula: 


C 


D 


2Cf 


R  <  10^, 


0  .5 


(3.5) 


in  which  is  the  drag  coefficient  for  a  flat  plate  with  a 
laminar  boundary  layer  (equation  3.4b),  and  Cp  is  the  drag 
coefficient  for  the  profile 


1  2 
-  PfU  c 


At  higher  Reynolds  numbers  turbulent  separation  occurs  except 
on  certain  forms  having  a  thickness  ratio  d/c  <0.4.  An 
approximate  expression  for  the  drag  coefficient  (except  for 
"laminar-flow"  forms)  is  then 

R  >  10^,  -  <  0.4  (3.6) 

c 

in  which  the  value  of  k  depends  primarily  upon  the  position  of 
the  point  of  maximum  thickness  xj^ ,  varying  approximately  from 
k  =  2.0  when  x^/c  =  0.3  to  about  k  =  1  .2  when  Xj^/c  =  0.5,  and 
Cj=  is  now  the  drag  coefficient  for  a  flat  plate  with  a  turbu¬ 
lent  boundary  layer. 

In  the  intermediate  range  of  Reynolds  numbers,  10  <  R  ' 

10^,  the  drag  coefficient  is  sensitive  to  the  pressure  distri¬ 
bution  along  a  particular  profile  and  to  disturbances  in  the 
boundary  layer.  The  turbulent  shear  stresses  downstream  from 
x™  tend  to  increase  the  orag  coefficient;  on  the  other  hand, 
if  laminar  separation  had  been  occurring  at  lower  Reynolds 
numbers,  it  would  now  be  replaced  by  turbulent  separation  or 
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no  separation,  which  would  tend  to  decrease  the  drag  coefficient. 
For  most  forms  of  thickness  ratio  d/c  >  0.10,  the  drag  coeffi¬ 
cient  in  this  range  is  less  than  that  given  by  equation  (3.5). 

A  typical  curve  showing  a  composite  of  the  measured  drag  coeffi¬ 
cients  vs .  Reynolds  number  of  various  forms  of  thickness  ratio 
d/c  =  0.25  is  shown  in  Figure  3.3. 

C)  Laminar  Flow  Profiles.  A  class  of  forms,  the  laminar 
profiles,  designed  to  maintain  a  laminar  boundary  layer  over  a 
large  part  of  the  surface  up  to  high  Reynolds  numbers  (R  =  10^), 
has  drag  coefficients  differing  from  equations  (3.5)  and  (3.6). 
These  forms  have  their  maximum  thickness  occurring  between  40 
and  65  per  cent  of  the  chord  from  the  leading  edge,  a  small  nose 
radius  of  curvature,  and  negative  curvature  near  the  sharp  trail¬ 
ing  edge  (Figure  3.4).  Such  a  shape  results  in  a  negative  pres¬ 
sure  gradient  up  to  about  the  point  of  maximum  thickness,  which 
consequently  stabilizes  the  laminar  boundary  layer  to  much 
higher  Reynolds  numbers  than  for  more  normal  shapes.  At  Rey¬ 
nolds  numbers  R  <  10^,  the  drag  coefficient  would  be  somewhat 
larger  than  that  given  by  equation  (3.5)  because  of  more  vio¬ 
lent  laminar  separation,  but  at  higher  Reynolds  numbers,  exceed¬ 
ing  R  =  10',  turbulent  separation  would  occur,  but  with  the 
boundary  layer  still  laminar  ahead  of  the  point  of  maximum 
thickness.  In  their  range  of  effectiveness  these  forms  have 
about  half  of  the  drag  of  conventional  profiles  of  the  same 
thickness  ratio. 

D)  Blunt  Forms .  The  drag  of  a  blunt  form  is  principally 
due  to  its  pressure  distribution,  which  differs  radically  from 
that  for  inviscid  flow  because  much  of  the  after  end  of  such  a 
body  is  immersed  in  a  zone  of  separation.  Separation  on  a  body 
occurs  either  at  slope  discontinuities  in  the  body  profile 
(sharp  edges)  or  at  positions  in  the  boundary  layer  where  a 
positive  pressure  gradient  is  beginning  to  cause  reverse  flow 
of  the  slowly  moving  fluid  near  the  wall.  In  the  latter  case 
separation  occurs  farther  downstream  in  a  turbulent  boundary 
layer  than  in  a  laminar  one  because  fresh  momentum  is  trans¬ 
ported  to  the  wall  by  the  process  of  the  random  turbulent 
motions . 

A  phenomenon  observed^in  the  flow  about  blunt,  two- 
dimensional  forms  is  the  periodic,  alternate  shedding  of  vor¬ 
tices  from  the  opposite  sides  at  the  rear  of  the  body,  forming 
in  the  wake  two  staggered  rows  of  uniformly  spaced  vortices, 
called  a  vortex  street.  It  was  shown  theoretically  by  von 
KdrmAn ,  by  analyzing  the  stability  of  such  an  array,  that  the 
ratio  of  distance  between  the  rows  to  the  distance  between  suc¬ 
cessive  vortices  in  a  row,  h/b,  satisfies  the  equation 


(3.7) 
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which  gives  tlje  approximate  solution  h/b  =  0,2806,  a  value  in 
very  good  agreement  with  experiment.  The  theory  also  yields 
an  expression  for  the  resistance  coefficient  of  the  body. 


in  which  d  is  the  maximum  width  of  the  profile,  and  u^  is  the 
velocity  of  the  vortices,  a  result  of  the  correct  order  of 
magnitude.  It  cannot  be  used  to  predict  drag,  however,  since 
the  values  of  b/d  and  Uy/u  cannot  yet  be  obtained  bv  theory. 

The  frequency  n  with  which  the  vortices  are  shed  is  of  great 
practical  interest  since  the  body  experiences  an  oscillating 
lift  as  the  circulation  about  it  varies  from  positive  to  nega¬ 
tive.  Both  theory  and  experiment  show  a  relationship  between 
the  drag  coefficient  of  a  body  and  its  dimensionless  frequency, 
S  =  nd/u,  called  a  Strouhal  number.  A  composi te,  empirical  fit 
to  the  data  from  a  wide  range  of  blunt  bodies  is  given  by  the 
formula 


S 


0  .21 

r  3/4 


(3.9) 
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In  the  range  of  Reynolds  numbers  10  <  R  <  10  ,  blunt 

bodies  with  sharp  edges  show  little  variation  in  drag  coeffi¬ 
cient.  The  coefficients  for  some  typical  cases  are  given  in 
Figure  3.5. 

Blunt  but  rounded  forms  show  a  remarkable  reduction  in 
drag  coefficient  over  a  small  range  of  Reynolds  numbers  in 
which  transition  from  laminar  to  turbulent  separation  occurs. 
The  designations  "s ubcritical"  or  "supercritical"  are  used  to 
indicate  the  Reynolds-number  ranges  in  which  laminar  or  turbu¬ 
lent  separation  occurs.  Curves  showing  the  variation  of  C’p 
with  R  for  a  circle  and  a  2:1  ellipse  are  shown  in  Figure  3.6. 
The  variation  of  the  drag  coefficient  with  thickness  ratio  of 
ellipses  is  given  for  both  the  subcritical  and  supercritical 
ranges  of  Reynolds  numbers  in  Table  3.1. 


Table  3.1  Drag  Coefficients  for  Elliptical  Sections 


Subcritical _ Supercritical 


1 

1.17 

0 .40 

2 

0.88 

0.20 

3 

0.42 

0.12 

4 

0.37 

0  . 0 

6 

0 . 30 

0  .  ;i  0 

8 

0.28 

0,12 
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If  the  vortices  shed  by  a  blunt  cylinder  of  width  d  cause 
it  to  oscillate,  the  width  of  the  resulting  vortex  street  is 
increased  approximately  by  the  double  amplitude  of  the  oscilla¬ 
tion,  2a,  and  the  resulting  drag  is  also  increased  approximate¬ 
ly  by  the  factor  1  +  2a/d.  Such  vibrations  commonly  occur  in 
cables  situated  at  an  angle  to  a  current. 

The  basic  principles  recited  in  this  Section  3.2  follow 
closely  the  presentation  by  L.  Landweber  in  the  Handbook  of 
Fluid  Dynamics,  Section  13  (Streeter,  1961). 


3 . 3  Terminal  Velocity 

If  an  object  begins  its  fall  through  sea  water  with  zero 
velocity  it  starts  accelerating  and  develops  speed  in  the 
direction  of  the  gravitational  force.  Simultaneously  with  the 
build-up  of  speed,  it  experiences  the  drag  force  which,  as 
equation  (3.4)  shows,  grows  with  the  square  of  the  velocity. 
Thus,  very  soon  a  point  will  be  reached  where  the  drag  force 
has  grown  as  large  as  the  gravitational  pull  and,  therefore, 
no  further  increase  in  the  velocity  is  possible.  The  velocity 
at  which  this  occurs  is  called  the  terminal  velocity  because 
this  is  the  highest  velocity  the  object  will  attain  and  main¬ 
tain  throughout  its  fall.  Stokes  calculated  a  terminal  vel¬ 
ocity  for  a  spherical  object  falling  in  a  viscous  fluid  when 
the  viscous  forces  predominate  over  the  inertia  forces .  The 
result  is  given  by  the  well-known  Stokes'  law  which  equates 
the  drag  force  Fq  encountered  by  a  sphere  cf  radius  r  with  the 
net  downward  force  exerted  by  the  gravitational  field: 
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where 


2 

y  is  the  dynamic  viscosity  of  the  fluid  (lb. sec/ft  ) 

Pf  is  the  density  of  the  fluid  ( lb . sec^ /ft^ ) 

Pg  is  the  density  of  the  sphere  (Ib.sec  /ft**) 
u  is  the  terminal  velocity  (ft/sec) 


Rearranging  terms,  we  find  the  terminal  velocity  u  given  by: 


u 


(3.10a) 


Equation  (3.10)  was  obtained  by  neglecting  the  inertia 
forces  in  the  Navier-Stokes  equation  and,  therefore,  is  only 
valid  for  very  small  Reynolds  numbers  (R  <  0.5).  Since  R  is 
given  by  ur/v,  where  v  is  again  the  kinematic  viscosity  (ft^/ 
sec),  Stokes'  law  is  applicable  only  for  very  small  velocities 
or  verj’  small  particles  falling  in  relatively  viscous  fluids. 
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In  general,  the  Reynolds  numbers  for  the  phenomena  we  are 
investigating  here  will  be  much  higher  (our  range  of  interest 
will  extend  from  R  =  1000  to  R  =  500,000).  Thus,  Stokes’  law 
is  not  applicable  to  the  problem  at  hand. 

However,  the  basic  approach  of  Stokes  for  calculating  the 
terminal  velocity  is  still  valid,  provided  we  introduce  the 
proper  driving  and  resisting  forces.  The  drag  or  resisting 
force  for  any  object  is  given  by  equation  (3.4).  The  driving 
force  for  any  free-falling  object  through  a  fluid  is  given  by 
its  weight  W  minus  its  buoyant  force  B.  Thus,  the  general 
equation  of  equilibrium  at  the  terminal  velocity  u  is: 

-  Cd  A  Pf  u^  =  ¥  -  B  =  V(Yo  -  Yf)  (3.11) 


where  V  is  the  volume  of  the  object,  Yq  its  unit  weight,  and 
all  other  terms  have  been  defined  before. 

Since  the  drag  coefficient  depends  on  the  shape  of  the 
falling  body,  it  will  be  convenient,  for  the  sake  of  illustra¬ 
tion,  to  select  a  shape  for  the  purpose  of  discussion.  Also, 
it  is  necessary  to  select  the  attitude  or  orientation  of  the 
object  during  fall,  and  the  Reynolds  number. 

The  procedure,  then,  is  one  of  successive  approximations: 
we  first  guess  the  terminal  velocity  u  and  calculate  the  re¬ 
sulting  Reynolds  number:  R  =  uc/v.  For  this  Reynolds  number 
we  can  find  the  drag  coefficient  from  a  diagram  such  as  Figure 
3.6,  provided  we  know  the  shape,  attitude,  and  trajectory  of 
the  falling  object.  Knowing  also  the  unit  weight  and  volume 
of  the  object,  we  have  all  the  terms  necessary  to  calculate 
the  velocity  u  from  equation  (3.11).  If  this  calculated  vel¬ 
ocity  is  different  from  the  initially  assumed  one,  successive 
approximations  are  made  and  the  cycle  is  repeated  until  assumed 
and  calculated  velocities  are  essentially  in  agreement. 

The  shape  of  the  object  of  most  interest  for  our  consider¬ 
ations  here  may  be  a  cylindrical  body  with  hemispherical  ends 
as  shown  in  Figure  3.7.  If  the  length  of  the  cylinder  is  given 
by  L  and  its  diameter  by  D,  the  overall  length  of  the  body  will 
be  L  +  D  and  the  slenderness  ratio  L/D  will  describe  the  geo¬ 
metric  configuration  of  the  object. 

Values  of  the  drag  coefficients  Cq  for  several  slenderness 
ratios  L/D  and  through  a  range  of  Reynolds  numbers  from  R  = 
1,000  to  R  =  500,000  now  have  to  be  calculated. 
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3 . 4  Cylinder  Falling  with  Axis  Horizontal 

If  a  cylinder  such  as  snown  in  Figure  3.7,  but  without 
tail  fins,  and  properly  balanced,  falls  through  a  fluid,  it 
will  maintain  its  axis  in  a  horizontal  position.  The  drag 
coefficients  for  a  cylinder  falling  vertically,  while  its  axis 
remains  horizontal,  are  taken  from  data  given  by  Prandtl  and 
Tietjens  (1957)  for  an  infinitely  long  cylinder.  For  a  finite 
length  cylinder,  the  drag  coefficient  will  be  smaller  because 
of  side  leakage.  This  is  so  because  flow  around  the  ends  into 
the  reduced  pressure  zone  above  the  cylinder  occurs .  This 
causes  a  different  pressure  distribution  which  results  in  a 
reduced  pressure  drag.  The  smaller  the  L/D  ratio,  the  smaller 
will  be  the  pressure  drag  and,  hence,  the  drag  coefficient. 

For  L/D  =  0  the  pressure  coefficient  must  be  that  of  a  sphere . 

Prandtl  and  Tietjens  also  evaluated  the  reduction  of  the 
drag  coefficient  due  to  the  finite  length  of  a  cylinder  for 
various  L/D  ratios  at  a  Reynolds  number  of  100,000.  Assuming 
the  reduction  of  the  drag  coefficient  at  other  Reynolds  num¬ 
bers  to  be  the  same,  we  calculated  values  of  the  drag  coeffi¬ 
cients  Cj)  for  various  slenderness  ratios  and  Reynolds  numbers. 
This  assumption  is  substantiated  by  other  data  given  by 
Schlichting  (1968).  The  values  thus  obtained  ar.e  given  in 
Table  3.2. 

Experimental  verification  of  the  values  in  Table  3.2  has 
been  carried  out  by  several  investigators  and  they  have  been 
generally  confirmed.  Table  3.2  summarizes  the  most  important 
results.  Thus,  using  equation  (3.4)  and  the  coefficients  of 


Table  3.2 

Drag  Coefficients  Cq  for  Cylinders  Falling  with  Axis  Horizontal 
Reynolds  Slenderness  Ratio  L/D 


Number 

L/D=«o 

L/D=50 

L/D=20 

L/D=10 

L/D  =  5 

L/D  =  3 

L/D=2 

L/D  =  l  .0 

L/1);=0 

1 

10.0 

8.15 

7.60 

6.7 

6.0 

5.9 

5  .6 

5 

.2 

10 

2.7 

2.20 

2.05 

1.81 

1.62 

1.59 

1.51 

1 

.40 

100 

1  .45 

1.18 

1.10 

.97 

.87 

.85 

.81 

.76 

1.1 

1,000 

1.00 

.82 

.76 

.67 

.60 

.59 

.56 

.52 

.48 

10,000 

1.10 

.90 

.84 

.74 

.66 

.65 

.62 

.57 

.40 

. 50,000 

1.20 

.98 

.91 

.80 

.72 

.71 

,67 

.62 

.48 

100,000 

1.20 

.98 

.91 

.80 

.72 

.71 

.67 

.62 

.44 

200,000 

1.15 

.94 

<.87 

.77 

.69 

.68 

.65 

.60 

.41 

300,000 

0.80 

.65 

.61 

.54 

.48 

.47 

.45 

.42 

.095 

500,000 

0.30 

.24 

.23 

.20 

.18 

.18 

.17 

.16 

.090 

1,000,000 

0.35 

.29 

.27 

.24 

.21 

.21 

.20 

.18 

.15 

'  i 


Table  3.2,  the  total  drag  force  on  a  cylinder  as  depicted  in 
Figure  3.7,  but  without  fins,  sinking  horizontally  is  given  by: 


-  Cq  D(L  +  -  itD)  pf  u' 


(3.12) 


3 . 5  Cylinder  Falling  with  Axis  Vertical 

If  tail  fins  are  added  to  the  cylinder  as  shown  in  Figure 
3.7,  or,  if  the  center  of  gravity  and  the  center  of  buoyancy  in 
the  cylinder  do  not  coincide,  the  cylinder  will  eventually  fall 
with  its  cylindrical  axis  vertical. 

We  assume  the  front  and  back  ends  of  the  cylinder  to  be 
hemispherical.  The  total  drag  for  the  cylinder  then  will  be 
equal  to  the  drag  on  a  sphere  of  diameter  D  plus  the  skin  fric¬ 
tion  drag  on  the  cylindrical  surface  and  on  the  tail  fins.  This 
skin  friction  drag  may  be  taken  to  be  equal  to  that  for  a  flat 
plate  having  an  area  and  length  equal  to  that  of  the  cylinder 
bodyr  It  is  taken  from  the  classical  Blasius  solution  for  vis¬ 
cous  drag  on  a  flat  plate  at  zero  incidence.  The  total  drag, 
hence,  of  the  cylinder  with  its  axis  parallel  to  the  velocity 
vector  is  given  by: 


^total  "  ^sphere  ^cylinder  ^fins 
The  drag  components  for  the  sphere  are  given  by: 


(3.13) 


i.  2  2 

^sphere  ~  8  ^sphere  ^  ®  Pf  u 


(3.14) 


Again,  the  dimensionless  drag  coefficient  Cgp^gpg  is  a  function 
of  the  Reynolds  number  R.  Some  values  for  Csphere  for  a  range 
of  Reynolds  numbers  1  <  R  <  10®  are  given  in  Table  3.3. 


Table  3.3 

Drag  Coefficients  Cp  for  Spheres  as  a 
Function  of  Reynolds  Numbers 


1 

10 

100 

1,000 

10,000 

50,000 

100,000 

200,000 

300,000 

500,000 

1,000,000 


'sphere 

28 

4.2 

1.1 

.48 
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Experimental  results  given  by  Schlichting  (1968)  are  shown 
in  Figure  3.8.  The  sudden  drop  in  the  drag  coefficient  from 
R  =  2  X  10^  to  R  =  5  X  10^  is  caused  by  the  change  in  the 
regime  of  flow  in  the  boundar*y  layer  going  from  laminar  to  tur¬ 
bulent  . 

The  viscous  skin  friction  drag  on  the  cylinder  surface  is 
given  by: 

^cylinder  "  ^cyl  'ff  D  L  Pf  u^  (3.15) 


The  drag  coefficient  is  the  same  as  that  for  a  flat  plate 

at  zero  incidence  and  is  given  by: 


Cgyi  =  1.328  Rjj 


(3.16) 


It  should  be  noted  that  the  Reynolds  number  Rl  for  this  case  is 
given  by: 


where  v  is  again  the  kinematic  viscosity. 

Since  the  drag  contribution  of  the  fins  also  may  be  con¬ 
sidered  to  come  from  the  skin  friction  on  a  flat  plate,  the 
effect  of  the  friction  drag  on  both  the  cylinder  surface  and 
the  fins  is  given  by: 


2 

Fcyi  =  .664  Rl  V  d  L  Pf  u  (3.17) 

and 

—V  9 

Ffins  =  .664  Rp  ^  Ap  pf  u  (3.18) 

It  should  be  noted  that  Rp  =  uZ/v  where  Z  is  the  length  of 
the  fins  and  Ap  is  the  total,  wetted  area  of  all  fins. 

* 

Thus,  the  total  drag  for  the  cylindrical  object  falling 
with  its  axis  vertical  is: 

^total  "  [TTD('g‘  Csphere^" .  664L(  Rl  )  "^)+.664(Rp)  ^Ap]p^u^  (3.19) 

Another  approximate  solution  for  calculating  the  drag  of  a 
cylindrical  body  in  three-dimensional  flow  parallel  to  the 
cylinder  axis  has  been  given  by  Glauert  and  Lighthill  (3955). 
They  present  their  results  in  terms  of  a  table  of  logarithmic 
values  of  Reynolds  number  vs.  drag  coefficient.  These  are 
given  in  Table  3.4. 
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Table  3.4 


Viscous  Drag  Coefficient  vs.  Reynolds  Number  for 
Finite  Cylinder  (after  Glauert  and  Lighthill) 


logic 

loglO^^d/Pf’^^ 

7ra^ ) 

-3.0 

2.64 

-2.5 

2.89 

-2.0 

1.16 

-1.5 

1.44 

-1.0 

1.73 

-  .5 

0.04 

0 

0.37 

0.5 

0.73 

1  .0 

1.11 

1.5 

1.50 

2.0 

1.91 

where 

2.5 

2.33 

3.0 

2.76 

a  =  D/2 

3 . 6  Protuberances  and  Rough  Surfaces 


If  the  falling  object  has  a  surface  that  is  hydraulically 
"rough"  or  has  protuberances,  the  drag  coefficient  may  or  may 
not  be  affected.  In  general,  bodies  with  sharp  edges  such  as, 
e.g.,a  flat  plate  at  ninety-degree  incidence,  are  quite  insen¬ 
sitive  to  surface  roughness  or  protuberances.  Bluff  bodies  on 
the  other  hand,  such  as  circular  cylinders  or  a  flat  plate  at 
zero  incidence,  are  very  sensitive  to  surface  roughness.  Page 
and  Warsap  (1930)  reported  the  effect  of  varying  surface  rough¬ 
ness  k/d  for  infinite  cylinders  with  their  axis  normal  to  the 
direction  of  flow.  Their  results  are  reproduced  in  Figure  3.9, 
giving  the  drag  coefficient  Cp  as  a  function  of  the  Reynolds 
number  and  the  roughness  k/d  where  k  is  th.e  height  of  the  pro¬ 
tuberances  and  d  is  the  diameter  of  the  cylinder. 


3 . 7  Correction  of  Drag  Coefficient  for  Yaw 

If,  for  example,  a  finite  cylinder  should  fall  vertically 
while  its  axis  has  a  constant  inclination  to  the  horizontal 
(yaw  angle)  of  magnitude  (Ji,  the  total  drag  may  be  roughly 
approximated  in  the  same  manner  as  is  done  for  swept-back  wings 
in  aircraft.  If  u  is  the  fall  velocity  and  c  the  chord  length 
measured  in  the  direction  of  the  velocity,  the  drag  coefficient 
Cp  corrected  for  yaw  is  given  by; 
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Cjj(R,(}i)  =  Cq(R  cos^4i,0)cos  (p 
where  R  =  uc/v. 


(3.20) 


3 . 8  Other  Shapes 


If  the  falling  object  has  a  shape  different  from  that  of 
a  sphere,  a  cylinder,  or  a  flat  plate,  the  corresponding  drag 
will  have  to  be  established.  For  complex  shapes  the  drag 
coefficients  will  have  to  be  found  either  by  a  series  of  model 
experiments  running  through  the  entire  range  of  Reynolds  num¬ 
bers  of  interest,  or  it  will  have  to  be  approximated  by  break¬ 
ing  up  the  complex  shape  into  several  simple  shapes  for  which 
drag  coefficients  are  known  and  given  in  Section  3.2  or  are 
computed  analogous  to  the  method  demonstrated  for  the  cylinder 
in  Section  3.4.  A  few  more  drag  coefficients  for  simple  shapes 
are  given  below. 

Detailed  procedures  for  determining  the  profile  drag  are 
given  by  Schlichting  in  Chapter  25  of  Boundary  Layer  Theory . 

The  drag  coefficients  for  flat  plates  at  90°  incidence  are 
constant  through  a  wide  range  of  Reynolds  numbers. 

For  a  circular  disk,  according  to  Wieselsberger : 

Cd  =  1.10,  for  10^  <  R  <  10® 


For  an  infinite  plate  strip: 

Cd  =  1.85,  for  10^  <  R  <  10® 

For  a  Joukowsky  aerofoil,  with  a  thickness  to  chord  ratio 
=  .05: 

Cd  =  3.6  R"^ 

For  an  airship  hull  with  flow  parallel  to  its  axis,  the 
drag  coefficient  varies  very  little  with  the  Reynolds  number, 
as  is  shown  in  Table  3.5.  This  Table  also  shows  the  drag 
coefficients  for  a  streamlined  strut  with  flow  normal  to  its 
axis  . 


Table  3.5 


Drag  Coefficient  for  an  Air  Ship  Hull 
and  a  Streamlined  Strut 


R 

Cd 

(airship  hull) 

Cd 

( strut ) 

10,000 

.12 

CD 

40,000 

.09 

.08 

100,000 

.07 

.06 

1,000,000 

.06 

.05 
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3 . 9  Terminal  Velocity  and  Fall  Trajectory 


Assuming  the  drag  coefficient 
the  terminal  velocity  can  be  found 
(3.11): 


for  the  object  to  be  known, 
by  rearranging  equation 


u 


o 


Vg(p  -  pf) 

1.41 - 

Cq  a  pf 


(3.21) 


where  p^  is  the  average  density  of  the  object;  p^  ■=  Yo/g. 

Since  in  most  cases  Cp  is  a  function  of  the  Reynolds  num¬ 
ber,  that  is  to  say,  a  function  of  u,  a  first  estimate  of  u  has 
to  be  made  for  picking  the  appropriate  drag  coefficient. 

By  successive  approximations,  the  correct  drag  coefficient 
will  eventually  be  found  corresponding  to  the  Reynolds  number 
determined  by  the  terminal  velocity.  Fortunately,  Cp  is 
almost  constant  for  a  wide  range  of  R. 

If  the  initial  velocity  is  zero  or  has  a  vertical  com¬ 
ponent  only,  and  if  ocean  currents  are  absent,  the  trajectory 
is  a  straight,  vertical  line  to  the  ocean  bottom.  If,  on  the 
other  hand,  the  initial  velocity  has  a  horizontal  component, 
or  if  horizontal  ocean  currents  are  present,  the  trajectory 
will  drift  horizontally.  In  general,  for  free  fall,  the  effect 
of  an  initial  velocity  will  be  relatively  small  because  the 
drag  forces  are  quite  high  and  they  rapidly  reduce  any  hori¬ 
zontal  velocity  component  to  a  negligible  value .  Ocean  cur¬ 
rents,  however,  may  cause  a  considerable  drift  of  any  long, 
free  fall  trajectory.  In  both  cases  the  amount  of  lateral 
drift  can  be  calculated. 


3.9.1  Drift  Due  to  Horizontal  Velocity  Component 

ime  t  =  0  at  the 
the  (horizontal) 
in  the  y-direction, 
law,  the  equation 


312  =  «  =  -  7=“y  .  3t 

where  Ay  is  the  projected  area  of  the  object  onto  the  plane 
normal  to  the  y-axis .  Assuming,  as  a  first  approximation,  the 
drag  coefficient  Cpy  to  be  independent  of  the  Reynolds  number, 
that  is  to  say,  independent  of  dy/dt,  equation  (3.22)  becomes 
a  second  order  differential  equation; 


Assuming  an  object  of  mass  m  is  at  t 
point  y  =  0  and  has  a  velocity  component  in 
y-direction  of  Vq.  The  profile  drag,  also 
is  assumed  to  be  Fpy.  Then,  using  Newton's 
of  motion  of  the  object  is  given  by: 

^  A„  2 
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For  a  typical  missile,  toroedo  or  instrument  package,  K  might 
vary  from  .005  to  .02  (ft“^). 

Having  the  solution  (3.25),  it  is  interesting  to  find  the 
distance  and  the  time  after  which  the  velocity  decreases  to  a 
negligible  value. 

Assuming  K  =  .01  (ft  )  the  excursion  values  y  and  the 
times  (in  terms  of  Vq)  for  a  decrease  of  the  velocity  by 
increasing  orders  of  magnitude  have  been  calculated  and  are 
given  in  Table  3.6.  We  note  that  the  drag  coefficient  was 
assumed  to  be  constant  which  is  the  case,  roughly,  for  Reynolds 
numbers  R  =  1000  to  R  =  200,000.  But  at  the  lower  velocities, 
the  drag  coefficient  drastically  increases.  Hence,  the  excur¬ 
sions  calculated  by  equation  (3.25)  are  higher  than  those  that 
actually  v;ould  occur.  The  values  of  Table  3.6  may,  therefore, 
be  considered  to  be  an  upper  bound.  Note  that  the  second 
column  gives  the  distance  (in  feet)  that  the  object  would  have 
traveled  if  the  drag  resistance  would  have  been  zero.  It  is 
quite  evident  that  horizontal  excursions  due  to  initial,  lateral 
velocity  components  will  be  quite  small. 


t 
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Table  h  6 


Typical 

Excursion  y 

and  Residual  Velocl 

t  y 

of  Ob)e 

ct  Falling  t 

hrough  Sea  Water  wi 

th  an 

Initial  Horizontal 

Velocity  Component 

Vo 

Residual 

Velocity  Y 

Vot 

y{  f t ) 

Vo 

1 

O 

900 

230 

Vo 

■ 

9 , 900 

400 

Vo 

*  10*^ 

99 , 900 

690 

Vo 

■  10”^* 

999 , 900 

920 

v„ 

O 

t 

9,999,900 

1  ,1S0 

3.9.2  Drift  Due  to  Ocean  Currents 


If  the  object  falls  through  water  having  a  horizontal 
current  component  V,  it  will  drift  off  a  Etraijtbt  trajectory 
due  to  the  profile  drag  produced  by  the  current. 

Assuming  again  an  object  of  mass  m  failing  through  water 
having  a  horizontal  current  component  i  7,  »he  equation  of 
motion  is  given  by: 


y  +  K{y)^  -  2KVy  =  -  KV^ 


If  the  current  velocity  is  constant  V  -  V^,  equation  f?.28) 
can  be  integrated. 

From  the  nature  of  the  problem,  it  is  evident  that  after 
some  time  the  falling  object  assumes  the  constant  horizontal 
velocity  component  V^.  Since  the  vertical  velocity  is  also 
constant  and  equal  to  the  terminal  velocity  Uq  ,  the  trajectory 
will  be  a  straight,  inclined  line.  Typical  trajectories  for  an 
object  falling  through  water  with  a  constant  horizontal  current 
are  shown  in  Figure  3.10, 

From  Figure  3.10  it  is  evident  that  the  continuously  act¬ 
ing  effect  of  a  horizontal  drag  component  can  cause  large  excur¬ 
sions  y  for  long  trajectories,  i.e.  for  the  object  falling 
through  deep  water. 
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3.10  Calculation  of  Terminal  Velocities  for  a  Typical  Object 

Assuming  a  typical  slenderness  ratio  for  the  cylinder  of 
Figure  3.7  of  L/D  =  10,  the  drag  coefficients  for  fall  with 
axis  horizontal  and  vertical  are  given  in  Table  3.7  and  shown 
in  Figure  3.11,  Curve  B.  Curve  A  shows  the  drag  coefficients 
for  an  infinitely  long  cylinder.  The  drag  coefficients  of 
Table  3.7  were  obtained  by  evaluating  equations  (3.12)  of 
Section  3.4  and  (3.19)  of  Section  3.5  respectively. 


Table  3.7 


Drag  Coefficients  for  Cylinder  with  L/D  =  10 


Reynolds  No. 

Cn  (axis  horizontal) 

CD(axis  vert.) 

Drag  Force 

Fd  =  I"  CdD(L  + 

n  .  2 

ij-  D)pfU 

1  ttD^  2 

2  Cd  —  Pf'^ 

10 

1.81 

20.8 

100 

.97 

6  .4 

1,000 

.67 

2.1 

10,000 

.74 

.93 

100,000 

o 

cx> 

.57 

200,000 

.77 

.52 

300,000 

.54 

.22 

500,000 

.20 

.17 

1  ,000,000 

.24 

.20 

We  can  now  calculate  the  terminal  velocities  of  a  cylinder 
as  shown  in  Figure  3.7  having  a  slenderness  ratio  L/D  =  10  both 
for  fall  with  axis  horizontal  and  vertical. 

The  terminal  velocity  according  to  equation  (3.21)  is: 


Hence,  the  terminal  fall  velocity  is  also  a  function  of  the 
volume  to  area  ratio  V/A,  and  the  ratio  of  the  relative  object 
density: 

PQ  ~  Pf 
Pf 

The  terminal  velocities  for  a  cylinder  with  a  slenderness 
ratio  of  L/D  =  10  have  been  calculated  and  are  shown  in  Figures 
3.12  and  3.13  for  the  cylinder  falling  with  its  longitudinal 
axis  horizontal  or  vertical  respectively. 
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4 . 0  Free  Fall  and  Powered  Trav  Penetration 

If  the  velocity  of  the  arrival  at  the  ocean  bottom  is 
known  or  has  been  calculated  by  the  method  discussed  in  Chapter 
3,  the  penetration  into  the  ocean  bottom  can  be  determined. 

The  problem  of  an  object  penetrating  into  a  target  mater¬ 
ial  is  the  classical  problem  of  terminal  ballistics.  An  im¬ 
mense  amount  of  literature  both  theoretical  as  well  as  experi¬ 
mental  exists  on  this  very  subject.  However,  most  of  these 
investigations  are  concerned  with  projectiles  striking  a  rela¬ 
tively  hard  target  at  considerably  higher  ballistic  velocities. 

The  earliest  known  penetration  equation  is  due  to  Robins 
and  Euler  who  assumed  that  the  resistance  of  a  medium  to  pene¬ 
tration  is  a  constant.  This  results  in  a  constant  deceleration 
of  the  projectile  and  an  increase  of  the  depth  of  penetration 
with  the  square  of  the  impact  velocity.  In  the  most  general 
formulation  one  may  assume  that  the  resistance  of  a  target 
material  to  penetration  is  the  combined  result  of  three  com¬ 
ponents,  fi(v°),  f2(v^),  f3(v2),  which  are  functions  of  the 
instantaneous  velocity  v.  Thus,  according  to  Newton's  second 
law: 

-  m  ^  =  f^(v°)  +  f2(v)  +  (4.1) 

Equation  (4.1)  essentially  postulates  that  during  the 
first  phase  of  penetration,  at  high  velocity,  the  motion  is 
governed  by  the  resistive  forces  proportional  to  the  square  of 
the  velocity.  These  forces  would  be  analogous  to  the  drag 
resistance  in  fluid  flow.  Subsequently,  at  "moderate"  veloci¬ 
ties,  the  resistance  will  be  proportional  to  the  first  power 
of  the  velocity  and,  hence,  during  this  phase  the  resistance 
is  analogous  to  the  viscous  resistance  in  fluid  flow.  In  the 
last  phase  resistance  is  primarily  caused  by  static  forces. 

If  the  functions  in  equation  (4.1)  are  assumed  to  be  con¬ 
stant  coefficients,  the  statement  is  simplified  to: 

dv  9 

-  m  =  a  +  3v  +  yv  (4.2) 

Equation  (4.2)  is  a  general  formulation  for  various  well- 
known  penetration  formulae.  If,  for  example,  one  takes  g  = 

Y  =  0,  one  obtains  the  Robins-Euler  formula.  For  the  case  of 
6  =  0 ,  we  obtain  Poncelet's  equation: 


m 


dv 

dt 


(4.3) 
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For  the  case  of  o  =  0,  we  obtain  the  equation  credited  to  Resal: 

dv  o  ,  . 

”  Tv  “  +  yv  (4.4) 

The  problem  with  all  these  equations  is  that  the  material 
constants  a,  3  and  y  have  to  be  established  experimentally.  If 
equation  (4.2)  or  any  of  its  special  cases  is  integrated  twice 
with  the  initial  conditions,  t  =  Ot  x  =  0,  v  =  v^,  the  final, 
maximum  penetration  can  be  calculated  by  setting  the  velocity 
V  =  0. 


The  resulting  ultimate  penetrations  for  the  Robins-Euler , 
Poncelet  and  Resal  equations  respectively  are: 


^max 

^max 

^max 


m  V 


2a 


m 


Yv, 


^Jind  + 


m 


Yv 


0.2 


***  *  i*  ^  \ 


2y 


e 


Robins-Euler 

Poncelet 

Resal 


(4.5) 

(4.6) 

(4.7) 


The  solution  for  the  general  equation  (4.2)  is  somewhat 
more  involved  and  depends  on  the  relative  magnitude  of  the 
coefficients  a,  6  and  Y* 

2 

For  the  case  where  6  <  4aY ,  integration  of  (4.2)  leads  to: 


m 

2Y 


£n 


YVo"+3Vo+a 


m3 


“  T”  A  (4aY  -  3  ) 
Yv^+3v+a  Y 


2  s-h 


t  =  2  A  m  (4aY  -  3^)"^ 


(4.8a) 

(4.8b) 


where  x  is  the  instantaneous  penetration  at  time  t  and  velocity 
V  and 

A  =  tan“^(2YvQ+3)(4aY-3^)”’^-  tan“^  (  2yv+3  )  ( 4aY-3^  )  ‘ '^ 

The  maximum  penetration  can  be  obtained  by  setting  v  =  0 . 


If  3'^  >  4aY,  integration  yields: 


m 


YVq  +3Vo+a 


X  -  2y  Yv2+8v+a 
t  =  2  B  m  (3^  -  4aY)"^ 


m3 


B(3^  -  4aY)“^ 


(4.8c) 

(4.8d) 


where 


B  =  x  X'n 


t2Yv+3+( 62-4aY)^] [ 2 YVq+B - ( 3^ - 4aY )^] 
2  [2yv+3-( 3^-4aY)^] [2YVo+6+( 6^-4oY)^] 
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For  the  special  case  where  =  4aY,  the  integration  yields: 


-  JL  0  6  . 

^  ■  2Y  g+2Yv  "  Y  ^ 


(4.8e) 


t  =  2  C  m 


(4.8f) 


where 


C  =  (2Yv+6)“^  -  (2YvQ+g)“^- 


For  this  same  case  g  can  also  be  expressed  as  g  =  2  Y^, 

and  equation  (4.2)  can  be  written  as: 


-  m  i:  (a^  +  Y^v)^ 
df^ 


(4.9) 


Then,  the  integration  with  the  given  initial  condition  leads  to: 
^  a^+Y^v_  .  . 


_  m  _ o  V  _ o _ 

^  Y  "  ot^+Y^v  “  ^  To^+Y^oTo^+Y^^ 


(4.10) 


and  the  maximum  penetration  is  given  by: 


m  Y 

=  r  And  +  “^  v^)  - 


(a'^  +  Y'^v^) 


(4.11) 


Also: 

t 


m 


V  -  V 

o 


(4.12) 


As  stated  before,  the  material  constants  a,  g,  y  have  to 
be  determined  for  each  material  and  impacting  object.  While 
these  formulae  have  been  used  in  expressing  a  wide  variety  of 
test  results ,  the->  T>  validity  has  been  tested  only  within 
limited  ranges  of  impact  velocities  and  often,  by  the  very 
nature  of  the  experiments,  the  data  sho;?  a  considerable  scatter. 


Also,  the  coefficients  a,  g,  Y  do  not  appear  to  have  any  clear 
relationships  to  the  basic,  mechanical  properties  of  the  target 
material.  At  any  rate,  such  a  relationship  has  not  yet  been 
sufficiently  explored. 


For  this  reason  we  consider  it  risky  to  apply  these  formu¬ 
lae  directly  to  the  ocean  bottom  penetration  problem.  Normally, 
it  is  highly  improbable  that  in  situ  tests  for  a,  g,  Y  can  be 
performed  at  reasonable  cost.  Before  this  is  attempted,  an 
analysis  should  be  made  how  the  coefficients  o,  g,  Y  are 
related  to  the  basic  material  properties  of  the  ocean  bottom 
soils . 
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Apart  from  terminal  ballistics  investigations,  a  group  of 
impact  penetration  studies  in  the  literature  are  concerned  with 
meteorite  impacts.  These  are  definitely  beyond  the  practical 
range  of  ocean  bottom  penetrations.  Meteorite  impact  usually 
occurs  at  hypervelocities  resulting  in  liquefaction  and  vapori¬ 
zation  of  the  target  material,  as  well  as  scattering  and  splat¬ 
tering  of  the  impacting  object.  These  phenomena  are  beyond  the 
scope  of  this  study. 

Other  investigations  reported  in  the  literature  were  con¬ 
cerned  with  the  soft-landing  problem  of  a  space  vehicle  on  an 
unexplored  soil  surface  . 

While  the  theoretical  aspects  of  an  object  striking  a  tar¬ 
get  material  do  not  appear  to  be  drastically  different,  cer¬ 
tainly  the  experimental  investigations  are  greatly  dependent 
upon  projectile  and  target  material  and  also,  on  the  impact 
velocity  and  energy. 

All  impact  phenomena  are  controlled  by  the  ratio  of  spe¬ 
cific  impact  energy  of  the  projectile  to  the  potential  strain 
energy  of  the  target  material.  Therefore,  even  purely  analyti¬ 
cal  studies  have  to  concentrate  on  those  phenomena  that  predom¬ 
inate  the  particular  problem  at  that  ratio.  Below  we  are  pre¬ 
senting  those  methods  that  appear  most  promising  and  applicable 
within  the  present  state  of  the  art  for  analyzing  ocean  bottom 
penetration  within  the  range  of  probable  parameters.  First, 
however,  we  must  point  out  two  fundamentally  different  phases 
of  the  penetration  process.  The  first  one  is  the  penetration 
phase  with  increasing  penetration  area.  (The  penetration  area 
is  defined  as  the  orthogonal  projection  of  the  contact  surface 
on  a  plane  normal  to  the  penetration  trajectory.)  This  phase 
occurs,  for  example,  during  the  initial  penetration  when  a 
sphere,  a  conical-  or  an  ogive-nosed  projectile  penetrates  into 
a  target  material.  During  this  initial  penetration,  the  pro¬ 
jected,  cross-sectional  contact  area  continuously  increases. 

At  a  penetration  of  p  =  D/2  for  a  sphere,  or  equal  to  the  length 
of  the  nose  cone  for  a  projectile,  the  maximum  projected  contact 
area  has  been  reached  and  any  penetration  thereafter  occurs  at 
a  constant  area.  Also,  a  flat,  cylindrical  punch  (such  as  a 
flat-ended  pile)  will  always  have  a  constant  penetration  area. 
Since  the  mechanics  of  the  two  problems  is  somewhat  different, 
we  have  to  make  a  distinction  between  increasing  area-  and  con¬ 
stant  area-penetration. 


4 . 1  Penetration  with  Increasing  Area 


Increasing  area  penetration  usually  occurs  as  the  initial 
phase  of  deep  penetration  or  for  shallow  penetrations  with  the 
penetrating  object  having  a  geometry  such  that  any  increase  in 
penetration  Ax  causes  an  increase  in  the  penetration  area  AA . 
This  condition  can  be  expressed  by: 
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Ax 

Studies  investigating  increasing  area  penetration  were 
carried  out  by  the  Aeronutronic  Division  of  Philco  for  develop¬ 
ing  a  lunar  penetrometer  (Anon.  1966)  and  by  Schmid  (1966)  and 
Tsai  and  Schmid  (1969)  for  determining  soil  properties  in  situ 
by  an  impact  penetrometer. 

These  investigations  analyze  the  penetration  of  a  spheri¬ 
cal  f  inetrometer  into  a  soil  material.  The  Aeronutronic  study 
integrated  the  resulting  equation  of  motion,  a  nonlinear  dif¬ 
ferential  equation  utilizing  a  Taylor  series  expansion.  The 
solution  was  given  in  the  following  form: 


where  v  is  the  instantaneous  velocity,  p  and  o  are  the  density 
and  compressive  strength  of  the  target  material  respectively. 
The  term  v^  represents  the  initial  or  impact  velocity  of  the 
penetrometer,  while  D  is  the  penetrometer  diameter  and  m  is  its 
mass  . 


Equation  (4.13)  is  valid  for  the  entrance  phase  of  the 
penetration  process .  That  is  for  x  <  D/2 . 

A  considerably  simpler  solution  was  obtained  for  motion 
within  the  second  phase,  i.e.,  after  the  nose  of  the  penetrome¬ 
ter  is  completely  embedded  and  the  surface  area  of  contact  is 
no  longer  varying  with  depth  of  penetration.  This  solution  is: 


(4.14) 


where  Vj^  is  the  velocity  at  the  end  of  the  entrance  phase 
obtained  from  equation  (4.13)  for  x  =  D/2.  Thus,  equations 
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(4.13)  and  (4.14)  define  the  instantaneous  velocity  of  the  pene¬ 
trometer  at  any  depth  of  penetration  x.  This  theoretical  pre¬ 
diction  of  penetration  behavior  has  been  verified  experimentally 
for  some  range  of  projectile  geometries,  impact  velocities,  and 
types  of  target  materials . 

Very  good  agreement  is  claimed  for  the  velocity  range  from 
500  to  6,000  ft/sec  and  for  target  materials  ranging  from  soft 
polyethylene  to  glass  fiber  reinforced  epoxy  resins.  However, 
for  the  impact  velocity  range  of  interest  in  ocean  bottom  pene¬ 
tration,  no  experimental  confirmation  of  the  validity  of  equa¬ 
tions  (4.13)  and  (4.14)  exists  and,  therefore,  they  must  be  used 
with  caution. 

Even  though  equation  (4.13)  is  rather  lengthy,  its  evalua¬ 
tion  does  not  present  any  problem  with  the  help  of  a  computer. 
Previous  work  by  Schmid  (1966)  also  dealt  with  the  penetration 
of  a  sphere  into  soils  . 

For  the  penetration  of  a  rigid  sphere  into  a  soft,  plastic 
target  material,  one  can  make  the  assumption  that  the  resistance 
to  penetration  offered  by  the  target  material  is  equal  to  a  con¬ 
stant  yield  stress  or  flow  pressure  p©  times  the  projected  area 
of  contact.  Hence,  using  the  terms  defined  in  Figure  4.1'; 

P  =  irp^a^  =  7rpQ(2Rx  -  x^)  (4.15) 

Thus,  if  m  is  the  mass  of  the  penetrating  object,  applying 
Newton's  law  furnishes: 

X  +  -jjj—  (2Rx  -  x^)  =  0  (4.16) 


If  the  penetration  is  small,  say  x  <  R/4,  the  second  term  in 
the  parentheses  becomes  relatively  small  and  can  be  neglected 
compared  to  the  first  term.  Then  the  solution  of  equation  (4.16) 
is  given  by: 


X 


Vq3  sin 


t 

3 


(4.17) 


where 


The  maximum  depth  of  penetration  occurs  at  the  end  of  the 
impact  period,  when 


.  t 

X  =  Vo  cos  j  =  0 


which  is  satisfied  if 
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£ 

2 


Thus,  the  impact  period  T  is: 

T  =  1.57  8 

and  the  maximum  penetration 

«max  =  ^o^  ^ 

2 

If  the  penetration  is  larger  than  R/4,  then  the  x  term  in 
equation  (4.16)  cannot  be  neglected  and  the  solution  of  the 
equation 

X  +  Ax  ~  Bx^  =  0  (4.18) 


leads  to  an  integral  that  can  be  evaluated  only  numerically: 

dx 


(Vq^  -  Ax^  +  j  Bx^)^ 


+  C 


where 


and 


A  =  2  RB 
ifPo 


B  = 


m 


If,  on  the  other  hand,  the  penetration  is  larger  than  R, 
we  have  the  case  of  constant  area  penetration  and  equations 
(4.15)  through  (4.17)  respectively  are  modified  to: 


11 

ttPqR^ 

(4.15a) 

o 

o 

(4.16a) 

X  + 

R"  =  0 
m 

It 

i  «2  2 

-  2  ’^Po^  +  '^o't 

(4.17a) 

and  penetration  period  T  and  maximum  penetration  respectively 
are  given  by; 


T  = 


m 


2 


and 


1  m 


^max  "  2  TrpQR2  ^o 


39 


We  note  that  these  results  now  correspond  to  the  solution  of 
Robins-Euler  (equation  4.5)  with  a  =  hPqR^  . 

Rather  than  assuming  a  constant  flow  pressure  type  of 
resistance  (i.e.,  an  ideally  plastic  target  material),  one  can 
also  assume  the  target  to  have  the  rheological  properties  of  a 
four-parameter  model  as  shown  in  Figure  4.2. 

This  approach  also  was  explored  by  Schmid  (1966)  who  cal¬ 
culated  a  penetration  quantity  b(t)  by  using  the  principle  of 
elastic-viscoelastic  analogy  using  Laplace  transforms .  The 
results  are  given  in  terms  of  a  deformation  parameter  b(t). 

This  parameter  had  to  be  chosen  because  the  elastic-viscoelastic 
analogy  requires  a  linear  relationship  between  the  load  and  the 
deformation  parameter.  Since  for  the  problem  of  penetration  of 
a  sphere  into  a  target  neither  the  penetration  *'x"  nor  the  con¬ 
tact  radius  "a"  (see  Figure  4.1)  is  directly  proportional  to 
the  load,  the  quantity  b  =  x*a  was  chosen  which  is  a  linear 
function  of  the  load  as  long  as  x  <  R,  i.e.  as  long  as  the 
penetration  is  an  "increasing  area"  penetration. 


The  resulting  differential  equation  for  the  four-parameter 
model  is: 


0  ®1  ®1 
ni  42 


+  — )Dt  + 


^1^2 

4^42 


T  =  G^D^Y  +  6;^  ^  DY  (4.19) 


Here  the  moduli  Gi  ,-  G2 ,  '41  and  42  are  the  spring  and  dashpot 
constants  respectively  of  the  four-parameter  model  target  mater¬ 
ial  (see  Figure  4.2),  t  and  y  represent  the  stress  and  strain 
tensors  respectively,  and  the  operator  D  represents  differenti¬ 
ation  with  respect  to  time:  D  =  d/dt. 

Equation  (4.19)  was  solved  for  an  impact  pulse  having  a 
sinusoidal  variation  with  time: 


P  =  Pj^  sin  wt 


(4 .20) 


Experimental  accelerometer  measurements  performed  by  a 
penetrometer  showed  that  for  most  soils  (e.g.  for  clays,  silts 
and  even  loose  sands)  the  impact  force  pulse  could  indeed  be 
approximated  by  a  sinusoidal  curve.  Figure  4.3  shows  such  a 
typical  accelerometer  curve. 

The  solution,  obtained  by  a  Laplace  transform,  is  given 
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where 


P  - 

^  6K  ’  ^  ^2 


OJh  =  «  Wo 

^  ni  ^ 


42 


•and 


1  \  1  /if  1  2 

(1)3  =  J  <  a)2+Y:{:y(  w^+ka)2  )C  ±  ^  <  tw2+j::jrp'(a)j^+ka)2  )]  ■ 

K  is  the  bulk  modulus  of  the  target  material. 
Also : 
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The  resulting  equation  (4.2,1)  is  rather  complicated  and  cumber¬ 
some  to  evaluate  . 

Later  Tsai  and  Schmid  (1969)  considerably  simplified  this 
solution  by  separating  volumetric  and  deviatoric  components  of 
the  deformation. 

Under  high  energy  impacL,  soils  behave  essentially  like  a 
liquid  material  and  consolidation  is  relatively  small.  It 
appears  reasonable  then  to  assume  that  the  deviatoric  deforma¬ 
tions  predominate.  This  means  that  the  volumetric  change  is 
negligible,  or  K  =  »,  during  impact.  This  assumption  yields 
r  =  0  also,  and  the  solution  of  equation  (4.21)  is  then  simpli¬ 
fied  to: 

,  s  ,  TTt  ’  ,  ^  nt  -^2^ , 

b(t)  =  — ^  [C^sin  —  -  (C2+C3)cos  —  t  C2  +  C3  e  n2J 

(4.22) 

where 


n 


Go 

+  (-^) 
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To  evaluate  soil  parameters  from  an  impact  penetrometer 
test,  the  theoretical  equation  (4.22)  and  the  experimental 
curves  were  matched.  By  using  the  principle  of  the  least 
square  method,  soil  parameters  were  obtained  for  the  theoret¬ 
ical  curve  giving  the  best  fit.  Since  these  soil  parameters 
can  have  physical  significance  only  for  non-negative  values, 
the  occurrence  of  negative  values  were  excluded  in  the  process 
of  curve  fitting. 

The  results  thus  obtained  showed  that  the  four-parameter 
model  (Figure  4.2)  could  be  simplified  to  a  Maxwell  model, 
since  G2  was  nearly  equal  to  zero  and  112  could  be  combined 
with  to  represent  a  single  dashpot. 

For  the  simplified  soil  model,  equation  (4.22)  simplifies 
further  to: 

^  r  i  .Tit  T  iTt.^  ,, 

b(t)  =  — ^  Sin  -f  +  7^  (i  -  cos  "Y")}  (4.23) 


where  T  is  the  penetration  period  (sec.).  This  simplified 
theoretical  curve  (equation  4.23)  was  then  matched  with  the 
experimental  curve  by  using  the  same  approach  that  was  dis¬ 
cussed  previously.  The  calculated  soil  parameters  are  shown 
in  Table  4.1. 

The  problem  when  using  any  one  of  these  rheological  solu¬ 
tions  invariably  is  the  determination  of  the  model  constants 
G  and  n*  The  rheological  parameters  G  and  n ,  in  essence,  may 
be  regarded  here  as  material  properties  corresponding,  in  con¬ 
ventional  soil  mechanics  terminology,  to  a  modulus  of  elasticity 
and  a  shear  strength  at  a  constant  strain  rate. 


4 . 2  Constant  Area  Penetration 

If  an  object  has  penetrated  into  a  target  material  up  to 
the  maximum  cross-sectional  area,  any  further  penetration  there¬ 
after  will  take  place  at  a  constant  area.  Formulating  the 
resistance  to  penetration  then,  according  to  equation  (4.2), 
and  assuming  the  coefficients  a,  B,  andyto  be  constants,  we 
see  that  the  total  resistance  is  made  up  of  a  static  component 
a,  a  viscous  component  Bv  and  a  drag  component  yv^ .  Such  a 
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Table  4.1  Calculated  Soil  Parameters  for  Soil  as  a  Maxwell  Model 

(after  Tsai  and  Schmid,  1969) 


11 

Unconfined 

Target  Material 

Remarks 

compressive 

(psi) 

( psi-sec ) 

strength  (psi) 

Sandy  Princeton 

Natural  soil 

262.5 

0.5885 

clay 

after  two  days 

264.3 

0.5721 

soaked  in 

251  .4 

0.5357 

heavy  rain, 

W  =  27.42% 

266.2 

0 .5436 

Princeton  red 

Dry,  relatively 

156.6 

0 .4402 

18.4 

clay 

loose  material. 

201.7 

0.4987 

W  =  2.65% 

182.8 

0.4587 

Silty  sand  at 

Natural  soil 

665.3 

0.8938 

31  .2 

soil  surface 

882.0 

1.0219 

Silty  sand. 

Natural  soil 

669.3 

0.9133 

42.5 

20"  below 

626  .9 

0.8810 

surface 

610  .0 

0.8287 

Sandy  silt 

Natural  soil. 

1178 

1.222 

62  .2 

and  clay 

dry 

1192 

1.234 

1331 

1.271 

Sandy  clay 

Field  recently 

150.1 

0.4093 

39  .7 

plowed  and 

225  .9 

0.5021 

disked 

201.5 

0.4726 

Crushed  rock 
h-1" 

Compacted 

1709 

1.5010 

CBR^  =  18.5 

Fine  Sand 

Natural  soil 

373.6 

0  .6466 

CBR^  =  5.8 

373.6 

0.6466 

Wax:  Cambar 

Material  used 

155.6 

0 .4249 

M-348  at  3lOC 

for  comparison 

151.1 

0.4111 

purposes 

164.1 

0.4288 

161.7 

0.4227 

1 


CBR 


California  Bearing  Ratio 


formulation  means  implicitly  that,  during  the  first  phase  of 
penetration  (the  object  having  high  velocity),  the  motion  is 
governed  mainly  by  resistive  forces  proportional  to  the  square 
of  the  velocity.  Subsequently,  at  ''moderate"  velocities,  the 
resistive  forces  will  be  proportional  mainly  to  the  first  power 
of  the  velocity.  Finally,  towards  the  end  of  the  impact  period, 
at  very  low  velocities,  the  resistance  is  caused  mainly  by 
static  forces. 

As  evident  from  Chapter  3,  the  free  fall  terminal  velocities 
are  so  low  that  the  drag  component  of  the  total  resistance  to 
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penetration  will  be  extremely  small  and  can  be  neglected.  This 
component  is  of  importance  for  the  penetration  of  targets  at 
very  high,  ballistic  velocities  which  occur  rarely,  if  ever,  in 
ocean  bottom  penetration.  Thus,  we  will  confine  ourselves  here 
to  the  evaluation  of  the  resistance  equation. 

P  =  a  +  3v  (4.24) 


4.2.1  Constant  Area  Penetration  with  Static  Penetration 
Resistance 


Assuming  the  target  material  to  be  made  up  of  an  ideally 
plastic  material,  the  resistance  to  penetration  is  given  by: 

P  =  A  •  Pq  (4.25) 

where  A  is  the  constant  area  of  the  penetrating  object  and  p^ 
is  the  bearing  capacity  of  the  target  material.  If  the  impact 
is  ideally  plastic  (no  rebound)  the  maximum  penetration  Xj^gj^ 
can  be  calculated  from  the  energy  relationship: 


T,  ^  2 

^  ^max  “  2 


(4.26) 


where  m  is  the  mass  of  the  impacting  object.  Thus, 


The 


The 


m  2 

^max  ■  2APq  '^o 

impact  period  can  be  calculated  from  the  equation 

=  APq 

solution  of  this  simple  differential  equation  is 


(4.27) 

of  motion: 

(4.28) 
given  by: 


A  2 

^  ^  Po  t  +  '^ot 


(4.29) 


which 


for  X 


0,  yields: 


and 


T 


Vom 

"  RPo 


’^max 


(4.30a) 


(4.30b) 
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The  results  (4.26)  through  (4.29)  are  exactly  as  those 
given,  by  (4.15a)  through  (4.17a)  except  that  here  the  penetra¬ 
tion  area  does  not  have  to  be  circular.  Thus,  the  maximum 
penetration  for  constant  area  penetration  into  perfectly  plas¬ 
tic  soil  is  directly  proportional  to  the  mass,  the  square  of 
the  impact  velocity,  and  is  inversely  proportional  to  the  pene¬ 
tration  area  and  the  bearing  capacity  of  the  ocean  bottom. 

The  determination  of  the  ultimate  bearing  capacity  requires 
information  about  the  shear  strength  characteristics  of  the 
ocean  bottom  (see  Chapter  2). 

According  to  Terzaghi  (1943),  the  bearing  capacity  of 
shallow  foundations  is  given  by: 

Po  =  1.3  cNjj  +  0.4  YBN^  +  YDfNq  (4.31a) 

for  a  square  footing,  and 

Po  =  1.3  cNc  +0.6  YRNy  +  YDfNq  (4.31b) 

for  a  circular  footing,  where  c  is  the  cohesion  and  Y  the  unit 
weight  of  the  soil.  The  values  Ng,  Ny  and  N_  are  influence 
values  depending  on  the  angle  of  friction.  B,  R  and  Df  are  the 
dimension  of  the  footing  or  the  depth  of  foundation  respective¬ 
ly. 

In  most  instances  the  friction  angle  of  ocean  bottom  soils 
will  be  zero.  Hence  the  bearing  capacity  can  be  expressed  as 
given  Ly  Tschebotariof f  (1951): 


=  5.52  c(l  +  0.38 


or  by  Skempton; 


=  5.0  c(l  +  0 


+  0.44  2-) 
B  L 


.2  pa  + 


0  .2 


(4.32a) 


(4.32b) 


Moore  (1962)  has  calculated  bearing  capacity  values  from 
shearing  strength  data  of  ocean  bottom  cores .  These  are  given 
in  Table  2.2. 


4.2.2  Constant  Area  Penetration  with  Viscous  Penetration 
Resistance 

Assuming  the  target  material  to  be  an  ideally  viscous 
substance,  the  resistance  to  penetration  for  a,  say,  circular 
penetration  area  is  given  by: 
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P  =  2  Tiayv  (4.33) 

where  a  is  the  radius  of  the  penetrating  area,  y  (poise)  is  the 
viscosity  of  the  target  material,  and  v  is  the  penetration  vel¬ 
ocity.  The  equation  of  motion  then  furnishes: 

-  mx  =  2  Trayx  (4.34) 

or  X  +  ax  =  0 

1 

where  a  =  2  Ttay  (4.35) 

The  solution  of  (4.35)  is  given  by: 

^o  -at 

x=~(l~e)  (4.36) 


By  inspection  of  equation  (4.36)  one  recognizes  that  the 
maximum  penetration  is  given  by: 


^max 


2'rray 


(4.37) 


which  occurs  for  a  value  of  t  =  “.  Similar  to  the  case  of  the 
foui -parameter  rheological  model,  the  problem  will  be  how  to 
assess  or  determine  the  viscosity  y  (g  sec/cm^)  of  the  ocean 
bottom  soil.  We  note  that  in  this  formulation  of  the  problem 
the  penetration  period  T  at  which  maximum  penetration  is 
achieved  is  infinite  which  means  that  the  deformation  really 
never  stops,  but  asymtotically  creeps  up  to  the  ultimate  value 
Vo/a. 

This  comes  from  the  very  formulation  of  the  problem  as  a 
purely  viscous  phenomenon.  Often  for  such  problems  the  retar¬ 
dation  time  T  is  a  more  concrete  measure  of  the  time  scale  of 
a  process.  Assume,  for  example,  we  were  interested  in  the  time 
at  which  t  =  t  =  1/a.  At  that  instant 

X  =  —  (1  -  e“^)  =  0.63  — 
a  a 


and,  from  (4.35) 


m 

^  “  2nay 

The  retardation  time  is,  hence,  that  time  after  which  the  major 
portion  (namely  63  per  cent)  of  the  deformation  in  a  purely 
viscous  process  has  taken  place. 
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•  2  •  3  Empirical  Penetration  Formulae 

Quite  a  large  number  of  empirical  penetration  formulae 
can  be  found  in  the  literature.  These  are,  however,  mostly 
concerned  with  high  impact  velocities.  For  example.  Goldsmith 
(1960)  reports,  in  addition  to  equation  (4.27),  also  the  fol¬ 
lowing: 


Im^.l  2 
^max  ~  2  +  1) 


(4.38) 


V  V 

^max  =  DB(^)''  for  0.1  <  ^ 
o  o 


<  1.0 


(4.39) 


where  D  is  the  diameter  of  a  projectile  and  c^  is  the  wave  vel¬ 
ocity  in  the  target  material.  The  values  ,  B^,  and  B  are 
values  determined  from  test  data.  The  value  B  xs  given  as 
approximately  2.5  and  n  =  1.4. 

Kornhauser  (1964)  reports  several  penetration  formulae 
found  in  the  literature.  For  penetration  into  rock,  Maurer 
and  Rinehart  (1960)  find: 


’^max  ”  ^1  n2 


(Vq  -  K2) 


(4.40) 


where  is  the  total  penetration  in  inches,  W  is  the  weight 

of  the  impacting  object  in  pounds,  D  is  the  diameter  of  the 
object  in  inches,  and  v^  is  the  impact  velocity  in  ft/sec.  The 
coefficients  Kj^  and  K2  are  material  parameters. 

Tolch  and  Bushkovitch  (1947)  recommend  for  small  pro¬ 
jectiles  striking  soft  rock: 


x„^^  =  4.6 


,1.83 


0.001  V, 


(4  .41) 


and  for  large  projectiles: 


=  ol  .53  (0.001  v^)' 


(4.42) 


Lang  (1956)  on  the  other  hand  reports: 


xmax  =  ^2  ^o- 


1.96/s 


0.065 


(4.43) 


47 


where  s  is  the  unconfined  compressive  strength  of  the  target 
rock  or  soil. 

Rinehart  (1960)  and  Palmore  (1961)  suggest  for  softer 
materials  such  as  soil,  the  following  type  of  equation: 

K ,  ^ 

^max  “  jj2  ^  ^2^0  ^  (4.44) 

which  is  very  similar  to  equation  (4.38)  given  by  Goldsmith. 

Kornhauser  reduced  the  data  reported  by  several  investiga¬ 
tors  to  the  common  basis  of  impact  by  steel  spheres  .  His  com¬ 
parison  is  shown  in  Figure  4.4. 

As  can  be  seen,  these  data  are  all  for  impact  velocities 
of  more  than  100  feet  per  second.  For  this  reason,  they  appear 
to  have  limited  applicability  to  ocean  bottom  penetration. 

Comparing  equations  (4.38)  and  (4.44),  we  note  them  to  be 
identical,  except  for  the  use  of  a  different  nomenclature. 

Also,  recalling  the  classical  equations  of  ballistics  of  Sec¬ 
tion  4.1,  we  find  that  these  two  equations  are  nothing  but  the 
Poncelet  equation. 

Comparing  the  penetration  formula  derived  in  Section  4.2.1 
assuming  an  ideally  plastic  target  material  (equation  4.30b) 
with  any  one  of  the  empirical  formulae,  we  find  that  equation 
(4.42)  is  basically  the  same,  differing  only  slightly  in  the 
exponent  of  the  impact  velocity.  This  exponent  is  2.0  in 
equation  (4.30b)  vs.  1.8  in  (4.42).  Most  surprisingly,  the 
equation  derived  by  assuming  an  ideally  viscous  target  mater¬ 
ial  (equation  4.37)  is  in  its  structure  identical  to  the  empir¬ 
ical  formula  (4.41)  and  also,  except  for  a  constant,,  agrees 
with  the  empirical  formula  (4.40). 

While  these  empirical  formulae  (4.40  and  4.41)  were 
obtained  by  evaluating  data  from  projectiles  striking  rocks  at 
relatively  high  impact  velocities,  the  results  indicate  that 
a  formulation  as  used  in  the  derivation  of  equation  (4.37)  has, 
perhaps,  the  greatest  promise  of  describing  properly  the  ocean 
bottom  penetration  problem.  Particularly  also,  because  the 
lower  velocities  and  the  much  softer  target  materials  of  the 
ocean  bottom  penetration  problem  must  be  a  much  closer  approx¬ 
imation  to  the  theoretical  assumptions  made  in  the  derivation 
of  equation  (4.37)  than  the  conditions  that  prevailed  in  those 
tests . 


■'  J 
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4.2.4  Empirical  Penetration  Results  of  Sandia  Penetrator  j 

I 

The  penetration  into  soils  and  soft  rocks  by  a  specially 
designed  penetrator  dropped  by  an  airplane  from  a  large  height 
has  been  investigated  by  a  research  and  development  program  of 
several  years  by  Sandia  Corporation.  The  device  is  described 
in  detail  in  Section  6.4,  and  a  cross-section  is  shown  in 
Figure  6.30.  The  experience  gained  in  this  development  program 
has  been  evaluated  on  an  empirical  basis  in  the  form  of  a  nomo¬ 
gram  (Figure  4.5).  Young  (1969)  shows  how  the  empirical  par¬ 
ameters  for  developing  this  nomogram  have  been  determined.  The 
nomogram  is  used  as  follows:  For  the  given  diameter  of  the 
penetrator  (for  example,  8  inches),  one  passes  horizontally  to 
the  left  until  the  horizontal  intersects  the  curve  correspond¬ 
ing  to  the  given  weight  of  the  projectile  (1000  lbs.).  From 
there  a  vertical  line  follows  upward  to  the  intersection  with 
the  curve  that  corresponds  to  the  impact  velocity  (1400  fps.). 

From  there  a  horizontal  to  the  right  intersects  with  the  line 
corresponding  to  the  soil  constant  of  the  target  material. 

This  soil  constant  must  be  known  either  from  a  prior  test  or 
it  must  be  estimated  from  experience  with  impacts  on  similar 
materials  (S  =  4.6).  From  this  intersection  a  vertical  line 
downward  follows  to  the  intersection  with  the  line  N  correspond¬ 
ing  to  the  nose  performance  coefficient  N  (N  =  1.11).  Finally, 
from  this  intersection  a  horizontal  line  to  the  left  indicates 
the  penetration  (104  ft.)  of  the  projectile. 

Young  states  that  the  accuracy  of  the  nomogram  depends 
most  strongly  on  the  accuracy  with  which  the  soil  constant  S  is 
determined . 

For  a  series  of  tests  reported  by  him,  the  predicted  pene¬ 
tration  depths  exceeded  actual,  observed  depths  by  20%  on  9% 
of  the  tests  and  by  25%  on  less  than  1.5%  of  the  tests. 


4.2.5  Computer  Programs  to  Calculate  Penetration 

Noh-  (1964)  and  Wilkins  (1964)  have  treated  the  problem 
of  an  object  or  projectile  impacting  upon  and  penetrating  into 
a  target  material  analytically  and  numerically  by  substituting 
the  differential  equations  of  motion  by  finite  difference  equa¬ 
tions  and  using  a  coupled  Euler-Lagrange  (CEL)  mesh.  The  cen¬ 
tral  feature  of  t.he  CEL  calculation  is  the  solution  of  a  time- 
dependent,  two-space  dimensional,  Eulerian  hydrodynamics  prob¬ 
lem  for  a  region  with  an  arbitrary,  polygonal,  approximating 
mesh  having  a  general  moving  (i.e.  time-dependant)  boundary. 

The  resulting  CEL-code  is  a  time-dependent,  two-space 
dimensional  computer  program  (for  a  compressible,  non-heat¬ 
conducting  material)  which  makes  it  possible  to  couple  an 
Eulerian  approximation  of  some  regions  of  a  material  with  a 
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Lagrange  approximation  for  each  of  the  adjacent  regions  of  the 
material.  The  major  problem  in  such  a  method  is  developing  a 
suitable,  Eulerian  formulation.  In  the  CEL  method  the  Eulerian 
boundary  is  defined  by  one  or  more  Lagrange  lines  and,  hence, 
is  defined  by  a  polygonal  traverse.  A  moving  traverse  inter¬ 
secting  the  fixed,  Eulerian  mesh  creates  irregular,  time- 
dependent  boundary  zones  which,  in  general,  can  be  defined  as 
the  union  of  closed  polygons.  The  problem  then,  for  the  Euler¬ 
ian  calculation,  is  one  of  approximating  the  differential  equa¬ 
tions  for  a  time-dependent,  polygonal  region  at  the  boundary- 
in  a  way  that  is  consistent  with  the  equations  in  the  interior. 
An  additional  requirement  demands  overall  consistency  of  the 
difference  approximations  for  both  the  Lagrange  and  Eulerian 
regions . 


The  CEL  method  was  pioneered  at  the  Lawrence  Radiation 
Laboratory,  Livermore,  California  under  sponsorship  of  the 
U.S.  Atomic  Energy  and  has  since  been  applied  to  a  series  of 
problems  from  nuclear  cratering  to  armor  penetration.  See, 
for  example,  Wilkins,  Honodel  and  Sawle  (1967)  and  Wilkins 
(1968).  The  resulting  program  was  called  the  HEMP  Code,  which 
solves  the  following  equations  of  continuum  mechanics  by 
finite  difference  methods. 

The  basic  equations  in  the  HEMP  Code  for  the  calculation 
of  elastic-plastic  flow  are: 

a)  Equations  of  motion  in  cylindrical  coordinates  with 
radial  symmetry  about  the  z-axis: 
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b)  Equation  of  continuity: 


—  =  —  +  JLL  t  — 
V  8z  3r  V 


c)  Energy  equation! 
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d)  Equations  of  State: 

1)  Stress  components: 
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2)  Strain  rate  components; 
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3)  Hydrostatic  pressure: 


■  p  =  a(n-l)  +  bCn-l)"^  +  c(n-l) 


1  P 


V  -  p. 


4)  Von  Mises  yield  criterion: 


2  2  2  2  0 
Si  +  S2  +83^-3  a/  <  0 


where  the  symbols  used  define  the  following; 

r,  z  -  space  coordinates, 
r,  z  -  velocity  in  r  and  z  directions, 
Oj,,a2»P0-  total  stresses, 

Tv,„  -  shear  stress. 
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Sj,,S2»S0  -  stress  deviators, 

SjL,S2»S3  -  principal  stress  deviators, 
y  -  shear  modulus, 

Oq  -  yield  strength  of  the  target  material, 
ej,,e2,e0  “  strains, 

6  -  corrections  for  rotation, 
p  -  hydrostatic  pressure  =  +  00)* 

V  -  relative  volume, 

W  -  internal  energy  per  original  unit  of  volume, 
p  -  density, 

^  2  V  2  A 

q  -  artificial  viscosity:  q  =  Cq  Pq^V^  W  * 

Cq  -  constant, 

A  -  zone  area, 

-  reference  density. 

The  dot  over  any  parameter  signifies  a  time  derivative  along 
the  particle  path. 

Noh  (1964)  describes  the  basic  steps  for  using  this  pro¬ 
cedure  as  follows  (pp.  125-126): 

"...  The  basic  idea  in  the  CEL  code  is  that  the  boundary 
3R  (of  the  region  R  =  U^R^  we  wish  to  approximate)  and  the 
curves  (which  separate  the  subregions  Ri)  are  to  be 
approximated  by  Lagrange  lines.  Thus,  the  moving  boundaries 
( 3R  and  the  D^)  are  made  to  correspond  to  some  Lagrange  lines 
in  the  Lagrange  approximating  grids  . 

"A  subregion  R^  which  is  approximated  by  the  Eulerian 
mesh  will  consequently  have  its  boundary  3Ri  prescribed  by 
the  Lagrange  calculations.  Thus,  the  Eulerian  calculation 
reduces  to  a  calculation  on  a  fixed  mesh  having  a  prescribed, 
moving  boundary  and,  therefore,  constitutes  one  of  the  central 
calculations  in  the  CEL  code. 

"The  CEL  code  consists  of  a  large  rectangular  Eulerian 
mesh  which  we  denote  by  E,  and,  depending  on  the  problem  we 
wish  to  solve,  from  one  to  six  separate  Lagrange  grids  which 
we  denote  by  Li(t),  i  =  1,  2,  ...,  6  [i.e.,  E  and  the  Lj;(t) 
are  the  sets  of  lattice  points  which  define  the  Eulerian  and 
the  Lagrange  grids,  respectively;  also,  the  Lagrange  grids 
are  denoted  as  functions  of  time  since  their  grid  points  move 
with  the  fluids  they  approximate]  . 


A 

Up  to  7000  mesh  points  are  provided  for  in  the  IBM  7094  ver¬ 
sion  of  CEL  and  up  to  20,000  mesh  points  can  be  used  in  the 
IBM  7030  (Stretch)  version  of  the  CEL  code. 


’’The  calculations  that  are  made  at  each  time  step  in 
the  code  are  divided  into  three  main  parts:  Lagrange  calcu¬ 
lations,  Eulerian  calculations,  and  a  calculation  which 
couples  the  Eulerian  and  Lagrange  regions  by  defining  that 
part  of  the  Eulerian  mesh  E  which  is  active  and  by  determin¬ 
ing  the  pressures  from  the  Eulerian  region  which  acts  on  the 
Lagrange  boundaries . 

n 

”We  suppose  that  we  are  at  the  nth  time  step  (t"  =  Z-Atj) 
of  the  CEL  calculation  and  that  we  wish  to  advance  all  quan¬ 
tities  by  one  time  step  to  the  time  .  We  assume  that  the 

state  of  the  fluids  (density,  energy,  etc.)  is  known  at  t^ 
and  also  that  we  know  the  positions  of  the  Lagrange  grids 
[i.e.,  we  have  determined  L^^  =  LjCt’^')]  .  In  addition,  we 
suppose  that  we  have  determined  that  subset  of  the  Eulerian 
mesh  E,  which  is  interior  to,  or  on  the  boundary  of,  those 
subregions  having  an  Eulerian  approximation.  We  denote 
this  subset  of  E  at  t  *  t’^  by  E’^  =  ECt’^).  The  calculations 
for  the  next  time  step  proceed  in  the  following  way. 

"The  first  calculation  uses  the  known  (t  =  t^)  state 
of  the  Lagrange  fluids  and  the  pressures  acting  on  the 
Lagrange  boundaries  to  solve  the  Lagrange  difference  equa¬ 
tions  for  each  of  the  grids  L.^.  The  solution  to  the  dif¬ 
ferences  equations  gives  us  the  t  =  state  of  the 

Lagrange  fluids  and  new  grid  positions  L^”'*’  . 

"Next  it  is  necessary  to  determine  the  set  E""*"^,  and 
this  is  done  by  one  phase  of  the  coupling  calculation  which 
uses  the  nst\’  grid  positions  We  are  then  in  a  position 

to  solve  the  Eulerian  difference  approximation  equations  for 
the  state  of  the  fluid  in  the  region  .  This  is  done 

in  the  Eulerian  phase  of  the  CEL  calculation. 

"Having  determined  the  state  of  the  Eulerian 

region,  the  second  phase  of  the  coupling  calculation  deter¬ 
mines  the  t^"*"^  pressures  which  act  on  the  boundaries  of  the 
Lagrange  grids  .  We  have  thus  advanced  all  of  the 

fluid  quantities  and  grid  positions  to  their  t  =  t”  ^  values 
and  this  then  completes  one  basic  calculational  cycle  (or 
one  basic  time  step  of  the  calculation). 

"The  description  of  the  CEL  calculation  will  be  com¬ 
pleted  by  giving  the  'start  up'  or  initial  (t  =  0)  state  of 
the  fluids  and  the  initial  positions  of  the  mesh  points  for 
the  grids  L^*^,  i  =  1,  2,  ...,  E  and  the  subset  E*^...." 

The  basic  HEMP-CEL  Code  has  been  modified  by  Physics 
International  Company  to  cover  projectile  Impact  at  rela¬ 
tively  high  velocities  onto  target  materials  of  essentially 
zero  strength  such  as  water.  This  modification  was  called 
the  PIPE  Code  and  is  applicable  when  the  projectile  and 
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I 

I  target  deform  rather  than  break  up  and  splatter.  Figure 

I  4.6  shows  the  penetration  calculations  of  the  basic  HEMP 

I  Code  at  four  different  instances  for  a  steel  cylinder  pene-c-’ 

I  trating  into  an  aluminum  plate.  In  Figure  4.7  the  mesh  grid 

I  of  the  PIPE  Code  for  a  conical  projectile  before  and  after 

!  entering  water  at  a  velocity  of  2.44  km/sec  25  ysec  after 

I  contact  is  shown. 

I 

I  The  basic  principle  of  the  HEMP  computational  procedure 

could  also  be  applied  to  penetration  calculations  for  the 
ocean  bottom.  However,  the  existing  codes  will  have  to  be 
modified  to  account  for  the  extremely  low  impact  velocities 
of  the  ocean  bottom  sediments  and  their  special  properties, 
such  as  the  development  of  pore  pressures,  the  increase  of 
strength  with  depth,  and  possible  thixotropic  or  rheological 
stressrstrain  behavior.  Such  a  computer  program  indeed 
should  be  developed. 


4.2.6  Powered  Travel  Penetration 

As  stated  in  Chapter  1,  one  can  make  a  distinction  of 
the  various  kinds  or  types  of  penetration  according  to  the 
forces  acting  during  the  penetration  process.  A  missile, 
ship,  or  submarine  may  still  be  under  power  while  striking 
the  ocean  bottom.  However,  in  most  cases  the  preponderant 
part  of  the  energy  exchanged  during  impact  is  most  likely  to 
come  from  the  kinetic  energy  of  the  missile  rather  than  from 
the  propulsion  system  during  the  short  period  of  impact. 
Therefore,  the  penetration  phenomena  for  free  fall  and 
powered  travel  penetration  will  be  essentially  the  same. 

The  only  possible  difference  may  be  that,  in  powered  travel 
penetration,  when  power  is  supplied  continuously  and  consti¬ 
tutes  an  important  source  of  energy  during  impact,  the 
energy  furnished  during  the  impact  period  must  be  added  to 
the  kinetic  energy  to  account  for  the  total  energy  being 
supplied.  This  can  best  be  done  by  adding  a  kinetic  energy 
correction  in  terms  of  an  equivalent ,  additional  velocity 
to  the  impact  velocity.  If,  for  example,  during  the  impact 
period  T  power  is  supplied  to  the  impacting  object  at  a  con¬ 
stant  rate  w,  the  total  energy  added  to  the  object  during 
penetration  is: 

E  =  wT  =  y  m(AV)^  (4.45) 

where  m  is  the  mass  of  the  object  and  AV  is  the  equivalent 
velocity : 

le. 

AV  =  ^  (4.46) 
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If  AV  is  now  added  to  the  original  impact  velocity  Vq,  : 

the  problem  can  be  treated  as  before  using  | 

vj  =  Vq  +  AV  (4.47)  I 

Alsot  the  object  may  have  a  greater  horizontal  velocity  com¬ 
ponent  than  in  free  fall  penetration.  | 

This  horizontal  velocity  component  or  the  lack  of  align-  | 

ment  of  initial  impact  point  and  center  of  gravity  on  the 
same  impact  trajectory  causes  a  phenomenon  called  whip.  For 
a  cylindrical  object  with  hemispherical  nose,  as  shown  in 
Figure  4.8,  the  impact  decelerations  in  terms  of  the  decel¬ 
erations  of  the  center  of  gravity  c.g.  are  given  as  follows: 

At  the  nose: 

*nose  ~  ^c.g.  ® 

At  the  tail: 

«tail  =  ^c.g,  -  sin  6 

where  6  is  the  whip  acceleration. 

The  resulting  impact  will  be  a  complex  phenomenon  of 
point-  and  broadside  impact  whose  relative  severity  will 
depend  on  the  obliquity  of  the  impact  0  and/or  the  horizontal 
velocity  component . 

It  is  because  of  the  whip  acceleration  and  the  impact 
alignment  problem  that  Schmid  (1966)  suggested  a  spherical 
shape  if  a  remote  sensing  impact  penetrometer  is  to  be 
developed.  This  would  eliminate  one  course  of  spurious 
data  and  results .  While  the  problem  in  water  is  not  as 
critical  as  it  is  in  air,  it  nevertheless  still  exists  and 
a  spherical  or  hemispherical  penetrometer  shape  rather  than 
a  cone  is  strongly  recommended  if  this  shape  can  be  adopted 
without  excessive  penalties  . 
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5 .0  Dynamic  Penetration  -  Penetration  Resistance  -  Bearing 

Capacity 

Dynamic  penetrarion ,  where  specific  forces  are  applied  to 
facilitate  the  penetration  process,  is  most  likely  to  occur 
when  an  object  is  forced  deeply  into  the  sea  floor  sediment, 
that  is  to  say,  when  a  pile,  or  a  caisson,  is  driven  into  the 
ocean  bottom.  Thus,  it  is  necessary  to  review  that  process 
of  penetration  in  which  a  long,  slender  member  is  forced  into 
a  soil  stratum.  This  is  the  problem  of  driving  resistance 
and  bearing  capacity  of  piles.  When  an  object  such  as  a  pile 
penetrates  into  a  foundation  material,  the  force  required  to 
achieve  penetration  is  called  the  ultimate  load  Q  or  the  bear¬ 
ing  capacity.  This  ultimate  load  Q  is  made  up  of  two  compon¬ 
ents:  the  load  transmitted  through  normal  stresses  at  the 

point  of  the  pile  Qp  and  the  load  Qg  transmitted  through  shear 
stresses  along  the  circumferential  surface  area  of  the  shaft 
of  the  pile  (see  Figure  5.1).  If  the  area  of  the  pile  tip  and 
pile  shaft  are  defined  as  Ap  and  Ag  respectively,  the  total 
ultimate  load  can  be  expressed  by: 

Q  =  Qp  +  Qs  =  Ap  Po  +  JVo^k  (S\.l) 

As 

where  p^  (psi  or  tsf)  is  the  bearing  capacity  of  the  soil  at 
the  tip  of  the  pile  and  Tq  is  the  shear  strength  of  the  soil- 
pile  interface.  Both  Pq  as  well  as  are  material  properties 
depending  on  the  shearing  strength  of  the  soil. 


5.1  Point  Resistance 


The  first  solution  for  calculating  p^  was  given  by  Prandtl 
(1920),  (see  Figure  5.2),  applying  the  classical  theory  of 
plasticity  to  the  problem  of  a  rigid  punch  penetrating  into  a 
rigid-plastic  solid.  Reissner  (1924),  Caquot  (1934),  Buisman 
(1935),  Terzaghi  (1943),  Meyerhoff  (1951),  De  Beer  (1948),  and 
Berezantsez  (1952)  used  the  same  approach  and  applied  it  with 
variations  to  the  bearing  capacity  of  soils.  A  recent  review 
of  their  solutions  was  presented  by  Vesic  (1967),  The  general 
formulation  of  the  bearing  capacity  was  given  by  Terzaghi  as: 

Po  =  +  qNq^q  +  j  yBliyQy  (5.2) 

where  c  and  y  are  the  cohesion  (shear  strength)  and  the  unit 
weight  of  the  soil  respectively,  q  is  the  overburden  pressure 
at  the  tip  of  the  foundation,  and  B  the  foundation  width,  N^, , 
Nq,  Ny  are  dimensionless  bearing  capacity  factors  for  an 


infinitely  long  strip  and  ?q,  shape  factors.  The 

coefficients  N  and  ?  are  functions  of  the  friction  angle 

For  a  cohesionless  soil  such  as  sand  or  gravel,  equation 
(5.1)  reduces  to; 


(5.3) 


For  a  frictionless  soil  (^  =  0)  the  coefficients 
0,  and  Cq  =  1  and  equation  (5.1)  reduces  to: 

Po  =  cN.5o  +  q 


=  1,  Ny  = 


(5.4) 


The  coefficients  N  were  evaluated  by  Prandtl-Reissner, 
Terzaghi,  Meyerhof f,  Jaky  and  de  Beer  sometimes  assuming  dif¬ 
ferent  yield  lines  or  shear  patterns. 

Figure  5.3  shows  the  failure  pattern  assumed' by  Reissner 
and  by  Caquo't  and  Buisman.  Figures  5.4  show  the  failure  pat-r 
terns  assumed  by  Meyerhof f,  and  Figure  5.5  shows  the  yield 
line  pattern  of  de  Beer  and  Jaky  where  the  yield  lines  revert 
to  the  shaft  of  the  deep  foundation. 

The  resulting  bearing  capacity  factors  N^,  Nq  and  Ny  as 
obtained  by  Meyerhoff  are  given  in  Figures  5.6,  5.7,  and  5.8. 

Inspection  of  equation  (5.2)  shows  that  the  bearing 
capacity  increases  linearly  with  depth  due  to  the  second  term: 
qNqCq  since  q  is  the  overburden  pressure  prevailing  at  the 
point  of  the  foundation.  For  uniform  soil  conditions  q  =  yd. 
Hence,  for  deep  foundations  this  term  should  furnish  the  major 
part  of  the  total  point  bearing  capacity  Qp .  Vesic  (1967)  com¬ 
pared  the  bearing  capacity  factors  CqNq  obtained  by  various 
investigators  either  by  theoretical  derivation  or  by  empirical 
deduction  with  data  from  experiments  both  for  an  infinitely 
long,  deep  strip  as  well  as  for  a  deep,  circular  footing  or 
pile.  His  results  are  reproduced  in  Figures  5.9  and  5.10 
respectively. 

As  can  be  seen,  a  considerable  spread  of  almost  one-  order 
of  magnitude  exists  between  the  influence  coefficients 
obtained  by  de  Beer-Jaky  and  those  of  either  Terzaghi  or 
Prandtl-Reissner.  Also,  a  comparison  of  Figures  5.7  and  5.9 
shows  that  the  results  of  Meyerhoff  for  Nq  are  identical  to 
those  of  de  Beer. 


5.2  Skin  Resistance  for  Deep  Penetration 


According  to  the  classical  concepts  of  the  sliding  resist¬ 
ance  offered  by  a  surface  in  contact  with  a  soil,  the  total 
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skin  resistance  is  made  up  of  friction  and  adhesion.  •‘Thus, 
for  a  cylindrical  shaft 


Qs  = 


Aga  +  C 


r 

t-n 


a  tan  5dz 


(5.5) 


where  a  is  the  adhesion  between  soil  and  the  foundation  shaft 
which,  for  uniform  soil  conditions,  is  usually  assumed  constant 
throughout  a  soil  layer,  C  is  the  shaft  circumference,  a  is  the 
normal  pressure  on  the  shaft  surface,  and  6  is  the  angle  of 
wall  friction  between  soil  and  shaft.  If  we  assume  the  normal 
stresses  essentially  to  be  "at  rest”  lateral  earth  pressures, 
we  can  calculate  o  as: 


a 


q(z) 


koz 


Yequiv 


(5.6) 


where  kQ  is  the  coefficient  of  lateral  earth  pressure  at  rest 
and  Yeq^j[y  is  the  equivalent  unit  weight  of  the  soil  from  z  » 

0  to  z . 

Thus,  assuming  k^,  <5,  and  a  to  be  constant  with  depth,  the 
total  skin  resistance  is  given  by: 


Qs 


Aga  +  Cd^Yequiv^o  ^ 
1 

Cd(a  +  2  dYequiv^o 


(5.7) 


Now,  the  total  resistance  to  penetration  can  be  formulated 
and  yields: 

Q  ~  ^tip^®^c^c  ^  ^avg^^q^q  ^  2^  YBNyCy) 


+  Cd(a  +  ^equiv'^^o 


(5.8) 


According  to  equation  (5.8)  the  resistance  to  penetration 
should  increase  with  the  square  of  the  depth  of  penetration  d 
for  large  values  of  d. 

The  validity  of  the  formulation  as  derived  above  for  pen¬ 
etration  into  sand  has  been  investigated  recently  by  Vesic. 

He  conducted  large-scale  model  tests  as  well  as  prototype  field 
tests  and  compared  these  results  with  the  theoretical  formula¬ 
tions.  His  findings  are  summarized  below. 
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5.3  Surface  Penetration  into  Sand 


If  we  apply  equation  (5.8)  to  the  penetration  of  a  foot¬ 
ing  at  the  surface  into  a  cbhesionless  material,  then, 
because  c  =  0  and  d  =  0,  the  only  remaining  term  is: 


Q  1 
Po  =  A  =  2 


(5  .9) 


Vesic  found  that  for  this  case  three  distinct  types  of  shear 
failures  could  develop  underneath  the  penetrating  surface 
depending  on  the  relative  density  of  the  sand  Dp.  For  Dp  > 
0.7  shear  failure  would  occur  suddenly  with  a  pronounced 
peak  on  the  load-penetration  diagram  when  the  penetration 
reached  about  7%  of  the  foundation  width.  This  type  of 
failure  was  accompanied  by  the  appearance  of  clearly  defined 
shear  failure  zones  and  bulging  of  the  sheared  mass  of  sand 
very  much  in  accord  with  the  phenomenon  described  as  general 
shear  failure  by  Terzaghi . 

For  intermediate  relative  densities  (0.35  <  Dp  <  0.70) 
no  such  sudden  failure  occurred  and  the  load-penetration 
diagram  gradually^ and  with  irregular  jumps^ moved  up  to 
higher  loads  . 

The  wiggles  in  the  diagram  are  caused  by  small,  local¬ 
ized  shear  failures  that  gradually  develop  starting  at  a 
penetration  of  about  8%  of  the  foundation  width.  At  15% 
penetration  a  shear  zone  boundary  becomes  visible  at  the 
surface.  This  phenomenon  was  described  by  Terzaghi  as  local 
shear  failure. 

On  loose  sand  (Dp  <  0.35)  penetration  occurs  without 
any  bulging  of  the  surface.  The  resistance  increases  as 
penetration  progresses  and  sudden,  small  shears  occur  as 
soon  as  penetration  reaches  about  6  to  8  per  cent,  again 
giving  rise  to  a  wiggly  load-penetration  curve.  This  phe¬ 
nomenon  was  called  punching  shear  in  an  earlier  report  by 
De  Beer  and  Vesic,  The  three  types  of  failure  as  demon¬ 
strated  by  Vesic  are  reproduced  in  Figure  5.11. 

The  magnitude  of  the  measured  resistance  to  penetration 
is  compared-  with  the  theoretically  predicted  values  in  Fig¬ 
ure  5.12.  Theoretical  and  observed  values  reported  by  Vesic 
for  the  bearing  capacity  coefficient  ( Caquot-Kerisel , 

1953)  are  shown  as  a  function  of  the  dry  unit  weight.  Accord¬ 
ing  to  Vesic,  the  experimental  values  are  generally  1.2  to  4 
times  higher  than  the  corresponding  theoretical  values  . 

This  is  in  accord  with  the  experience  of  other  investigators 
(Hansen  and  Odgard,  1960;  De  Beer  and  Vesic,  1958;  and 
L’Herminier  et  al . ,  1961),  although  a  satisfactory  explana¬ 
tion  for  this  phenomenon  is  not  yet  available. 


5 .4  Penetration  into  Sand  at  Shallow  Depths 

The  same  three  characteristic  types  of  failure  that- 
occur  in  surface  penetration  also  are  observed  at  shallow 
depths.  As  the  depth  of  embedment  d  increases,  the  relative 
density  Dj,  at  which  the  general  shear  failure  changes  into 
local  shear  and  then  into  punching  shear  also  increases . 

The  approximate  limits  of  Dp  where  these  changes  take  place 
are  shown  in  Figure  5.13.  This  figure  shows  relative  depth 
d/B  versus  relative  density  Dp  for  the  three  types  of  fail¬ 
ure  . 


From  Figure  5.13  we  can  note  that  there  are  certain- 
depths  below  which  only  a  ;punching  type  failure  occurs.  For 
a  circular  foundation  this  is  around  d/B  =  4.4;  for  a  long, 
rectangular  foundation  it  occurs  beyond  d/B  =  9.4. 

The  limits  delineated  in  Figure  5.13  also  depend  on  the 
compressibility  of  the  soil.  In  general,  the  more  compress¬ 
ible  a  soil,  the  smaller  will  be  the  relative  depth  at  which 
say,  general  shear  failure  changes  into  local  shear  failu^re. 
Therefore,  some  materials  may  only  exhibit  punching  shear 
failure . 


5 .5  Penetration  into  Sand  at  Large  Depths 

According  to  equation  (5.8),  the  penetration  resistance 
for  a  pile  in  sand  should  increase  continuously  at  least 
linearly  with  depth  d  and  as  d^  for  large  values  of  d. 

Practical  observations  show  that  at  shallow  depths, 
penetration  resistance  indeed  does  increase  linearly  with 
depth  up  to  about  a  relative  depth  of  d/B  =  4.  However, 
thereafter,  for  larger  depths,  the  bearing  capacity  or  pene¬ 
tration  resistance  increases  at  a  smaller  rate  until,  at 
about  d/B  =  15,  any  increase  in  penetration  resistance  with 
depth  ceases .  This  certainly  is  in  conflict  with  any  theor¬ 
etical  prediction  of  bearing  capacity.  However,  the  obser¬ 
vation  is  confirmed  by  several  investigators.  Kerisel 
claims  this  is  due  to  a  complex  dependence  of  Nq  on  the 
friction  angle  (ji,  d/B  and  B.  Vesic  attempts  to  explain  it 
by  redefining  the  quantity  q  as  qf,  the  effective  normal 
stress  at  failure  acting  on  a  plane,  hprizontal  element  next 
to  the  foundation  base.  For  depths  d/B  >  15  this  stress  qf 
remains  constant  while  at  d/B  <  4,  q^  is  equal  to  the  over¬ 
burden  q.  Whatever  the  reason,  the  practical  fact  and  exper 
ience  remains  that  penetration  resistance  and  bearing  capa¬ 
city  in  a  homogeneous  sand  increase  first  linearly  with  the 
depth  of  embedment  up  to  about  d/B  =  4.0,  then  it  increases 
less  rapidly  until  at  about  d/B  =  15  no  further  increase  in 
the  penetration  resistance  occurs.  This  latter  fact  can. 
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perhaps,  also  be  explained  by  the  facts  that  (a)  for  larger 
.values  of  d/B  only  a  localized  punching  type  of  failure 
occurs  (see  Figure  5-13),  and  (b)  in  the  vicinity  of  the 
pile  point  actual  tensile  stresses  are  produced.  While  the 
soil  immediately  underneath  the  point  is  severely  compressed, 
a  certain  stress  relief  due  to  the  tensile  stresses  in  the 
vicinity  is  possible.  Figure  5.14  shows  the  strain  distri¬ 
bution  around  a  pile  point  in  deep  penetration  as  reported 
by  Robinsky  and  Morrison  (1964). 

Vesic  derived  an  expression  for  the  bearing  capacity 
coefficient  for  punching  shear  failure,  using  a  shear  fail¬ 
ure  pattern  shown  in  Figure  5.15,  and  obtained: 

Nq  =  exp(3.8  ^  tan  (|))tan^(45  +  y)  ■  (5.10) 

A  comparison  between  these  theoretical  values  and  experi¬ 
mental  values  observed  with  circular  model  foundations  is 
shown  in  Figure  5.16.  The  shape  factors  c  (equation  5.2) 
used  for  evaluating  the  test  results  were  ?  =  3.0.  The 
results  show  good  agreement  with  the  values  of  equation 
(5.10)  at  low  and  intermediate  relative  densities  while  at 
high  relative  densities  the  test  results  are  in  better 
agreement  with  the  values  for  Nq  obtained  by  the  Prandtl- 
Reissner  theory.  The  bearing  capacity  coefficients  Nq  based 
on  general  shear  as  given  by  Meyerhof,  De  Beer,  and  Jaky  are 
about  one  order  of  magnitude  too  high  when  compared  with 
experimentally  observed  values.  Another  fact  that  supports 
this  observation  is  the  experience  with  actual  pile  driving. 
If  the  general  shear  bearing  capacity  factors  were  correct, 
it  would  be  practically  impossible  to  drive  piles  into  sand 
at  large  depths.  This,  however,  is  done  quite  frequently. 
Since  ocean  bottom  soils  usually  are  very  loose,  with  low 
relative  densities,  bearing  capacity  factors  Mq  obtained  from 
equation  (5.10)  appear  more  appropriate  for  ocean  bottom 
penetration . 


5 .6  Surface  Penetration  into  Frictionless  Material 

Sandy  and  silty  soils  are  terrigenic  ocean  bottom  soils 
and  will,  therefore,  be  encountered  mostly  close  to  shore. 

The  pelagic  bottom  siediments  on  the  other  hand,  mostly  oozes, 
and  terrigenic  clays,  may  be  found  farther  offshore  and  in  the 
deep  ocean.  Pelagic  sediments  may  be  considered  essentially 
frictionless,  having  only  a  cohesion  shear  strength  compon¬ 
ent.  For  these  materials  the  original  Prandtl  solution  may 
be  applied  directly  (Figure  5.2).  Applying  equation  (5.8) 
to  a  frictionless  material  at  the  surface  ( d  =  0 ,  <})  =  0 , 
hence,  Ny  =  0),  we  find  it  reduced  to: 

Q  —  ^tip  ^^c^c 


(5.11) 
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The  value  of  can  be  taken  directly  from  Prandtl's 
problem  shown  in  Figure  5-2  and  is  given  by: 

Nq  =  (2  +  tt)  =  5.14 

The  value  of  the  shape  factor  for  both  a  square  or  a  circular 
cylinder  is  usually  taken  as  =  1.3. 

Hence j  the  total  penetration  resistance  for  a  square  or 
circular  object  penetrating  into  the  surface  of  the  ocean 
bottom  is 


Q  =  Ap  s  5-14  X  1.3  =  6.7  s  Ap  (5.11a) 

where  s  is  the  shear  strength  of  the  ocean  bottom  soil  and  Ap 
is  the  projected  area  of  penetration.  One  of  the  main  prob¬ 
lems  in  using  this  formula  will  be  in  the  determination  of 
the  proper  value  of  s . 


5 .7  Penetration  into  Frictionless  Material  at  Depth  d  >  0 

For  values  of  d  >  0,  we  find  Nq  =  1.0,  ^q  =  1,  and 
equation  (5.8)  becomes: 

Q  =  A^^p  (cN^Cc  +  q)  +  Cda  (5.12) 

We  are  concerned  here  with  The  penetration  resistance 
and  not  with  bearing  capacity  for  a  static  load.  Therefore, 
any  adhesion  is  most  unlikely  to  develop  during  the  driving 
process  because  of  remolding.  Thus,  the  second  term  of 
equation  (5.12)  may  be  neglected. 

According  to  Henky’s  theorem  for  the  theory  of  plas¬ 
ticity,  the  bearing  capacity  coefficient  for  a  long  strip 
penetrating  into  plastic  materials  at  d  >  B  is  given  by: 

Nc  =  2  +  2tr  =  8.28  (5.13) 

The  values  are  also  given  by  Figure  5.6  for  ij)  =  0 . 

Now,  the  penetration  resistance  may  be  written  as: 

Q  =  A^£p  (c  8.28  X  1.3  +  q) 

=  Atip  (10-8  =  +  Vavg'3)  (S-lt) 

Practical  experience  indicates  that,  in  general,  the 
bearing  capacity  factor  as  given  by  equation  (5.13)  is  a 
good,  reliable  value  for  the  penetration  resistance.  While 
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values  of  as  low  as  5  and  as  high  as  20  have  been 
reported,  the  preponderance  of  experimental  and  field  data 
supports  the  value  as  given  by  equation  (5.13)  or  by  Figure 
5.6  at  <p  =  0.  Also,  it  is  not  at  all  clear  how  much  of  the 
extreme  deviations  from  equation  (5.13)  are  attributable  to 
errors  in  the  determination  or  deviation  from  the  assumed 
value  for  c.  It  appears  then  that  these  values  can  be  used 
with  reasonable  confidence  for  the  determination  of  ocean 
bottom  penetration  for  frictionless  soils. 


5 . 8  Penetration  Dynamics 


If  an  object  is  to  penetrate  into  the  ocean  bottom  the 
penetration  resistance  as  developed  in  the  soil  and  calcu¬ 
lated  by  the  prior  articles  of  this  Chapter  has  to  be  over¬ 
come  .  This  can  be  accomplished  by  the  application  of  inter- 
mitteiit,  high  energy  impact  blows  as  is  done  in  conventional 
pile  driving;  it  may  be  accomplished  by  a  high,  statically 
applied  jacking  force  which,  of  course,  requires  correspond¬ 
ing  counterweights  or  reactions;  or  it  may  be  achieved  by 
continuously  pulsating,  vibratory  forces  (vibratory  driving) 

For  all  cases  of  penetration,  the  driving  force  F  must 
be  greater  than  Q  and,  if  so,  the  motion  of  the  penetrating 
object  could  be  given  by  Newton's  equation: 


mx  =  F  +  W 


(5  .15) 


where  m  is  the  mass  of  the  object  and  W  its  weight.  Since 
both  F  and  Q  may  be  functions  of  time  (5.15)  can  also  be 
written : 


X  = 


F(t)  -  Q(t) 
m 


+  g 


(5.16) 


In  order  to  solve  this  equation,  the  time  variations  of  F(t) 
and  Q(t)  will  have  to  be  specified.  For  impact  type  pile 
hammers  F(t)  depends  on  the  energy  of  the  hammer  and  the 
cushioning  of  the  blow.  These  factors  will  be  discussed 
briefly  in  Chapter  6  in  the  sections  discussing  penetration 
equipment . 

If  the  penetration  is  achieved  by  jacking,  F(t)  is,  of 
course,  the  jacking  force. 

For  vibratory  driving,  if  the  driving  force  is  gener¬ 
ated  by  counterrotating  eccentrics,  the  dynamic  force  is 
given  by: 


F(t)  =  m 


sin  ojt 


(5.17) 


) 
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where  mo  is  the  mass  of  the  eccentrics,  r^  the  eccentricity, 
and  0)  is  the  angular  velocity  of  rotation.  Thus,  equation 
(5,16)  becomes: 


nio^'o 


O 

Sin  (Jit  - 


Q(t) 


(5.18) 


The  principles  of  vibratory  pile  driving  are  briefly 
discussed  in  Section  5.9. 

The  resisting  force  Q(t)  may  also  be  a  time  dependent 
function.  As  in  Chapter  4,  one  can  attempt  a  formulation 
of  Q(t)  according  to  Poncelet*s  law: 

Q(t)  =  Qo  +  Qix  +  Q2(x)^  +  .  .  .  (5.19) 

The  static  components  Qq  of  the  penetration  resistance 
for  the  various  cases  was  discussed  in  the  preceding  articles 
of  this  Chapter. 

Very  little  information  is  available  on  the  viscous 
component  of  equation  (5.19).  For  sands  at  high  vibra¬ 

tory  energy  input  experimental  data  as  well  as  practical 
experience  suggest  that  this  second  term  may  be  negative. 

The  vibrations  in  sand  cause  a  softening  or  liquefaction  of 
the  soil  around  the  pile  such  that  penetration  is  enhanced. 
This  phenomenon  is  the  cause  for  the  remarkable  success  of 
vibratory  pile  driving. 

Cohesive  soils  such  as  pelagic  clays  and  oozes  are 
likely  to  have  the  rheological  properties  of  a  Bingham 
body.  That  is  to  say,  they  will  have  a  definite  yield  or 
shear  strength  and  beyond  that  a  small  viscous  component 
increasing  with  the  penetration  velocity  x.  The  shear 
strength  Tq  and  viscosity  of  kaolinite  and  bentonite  soil 
pastes  were  investigated  by  Wenz  (1963)  who  found  that  the 
shear  strength  was  very  sensitive  to  the  water  content,^ 
and  that  the  additional,  viscous  strength  increase  due  to  x 
was  small  for  kaolinite  clays  and  absolutely  negligible  for 
bentonite  clays  . 

Recently,  Goble  and  Scanlan  (1967)  and  Tomko  (1968) 
have  attempted  to  use  Poncelet's  law  (equation  5.19)  to 
determine  experimentally  the  static  bearing  capacity  Qq . 

The  basic  principle  chosen  is  to  substitute  equation  (5.19) 
into  equation  (5.15): 

mx  =  F(t)  +  W  -  (Qq  +  Q^x  +  Q2X^  +  ...)  (5.20) 
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If  is  measured  by  an  accelerometer  and  F(t)  by  a 
set  of  strain  gages*  say,  at  the  top  of  the  pile,  then  the 
driving  resistance  Q(t)  can  be  evaluated  in  principle  since: 

Q(t)  =  (Qo  +  Qix  +  Q2x2  +  ...)  =  F(t)  +  W  -  mx 

(5 .20a) 

Since  the  weight  and  mass  of  a  pile  are  generally  known,  all 
terms  on  the  right-hand  side  of  equation  (5.20a)  are  known 
or  measured  experimentally  and  Q(t)  thus  can  be  evaluated. 

The  static  bearing  capacity  is  particularly  conven¬ 
ient  to  determine  by  considering  the  measured  data  at  an 
instant  tQ  where  x  =  0.  Then  equation  (5.20a)  reduces  to: 

Qo  -  W  =  F(tQ)  -  m  xCt^T  (5.20b) 

The  practical  use  of  this  principle  for  field  applica¬ 
tion  requires  dynamic  instrumentation  that  at  this  stage  of 
development  is  probably  still  too  sensitive  and  perhaps  too 
complicated  for  general  field  application  on  a  normal  con¬ 
struction  site.  Such  instrumentation  has  to  be  rugged,  yet 
sensitive  and  reliable,  to  serve  under  the  varying  field 
conditions  of  the  construction  industry.  Tomko  (1968) 
reports  and  evaluates  a  series  of  tests  using  such  instru¬ 
mentation  and  reports  moderate  success.  His  final  conclu¬ 
sions  are  perhaps  worth  recording  here:  "The  present  study 
has  resulted  in  an  intense  familiarization  with  the  rapidly 
passing  response  profiles  occurring  over  a  few  tens  of 
milliseconds  after  the  hammer  impact.  That  these  profiles 
are  both  believable  and  physically  explainable  is  in  the 
author's  opinion,  to  be  considered  a  step  forward  in  pile 
dynamic  action." 

Nevertheless,  even  though  this  method  of  assessing  the 
bearing  capacity  of  a  pile  cannot  yet  rely  on  wide  exper¬ 
ience  and,  at  this  time,  is  not  yet  strongly  underpinned  by 
empirical  verification,  it  appears  extremely  interesting 
since  it  may  make  it  possible,  at  last,  to  measure  the 
static  bearing  capacity  of  a  pile  without  resorting  to  a 
pile  load  test.  Since  load  tests  require  large,  often 
expensive  reaction  force  -  particularly  in  marine  applica¬ 
tions  -  the  continued  development  and  full  proof -testing  of 
this  method  is  clearly  warranted.  With  the  ever-increasing 
pile  dimensions  and  pile  loads  in  offshore  platform  con¬ 
struction  and  the  absence  of  established  and  empirically 
verified  relationships  between  driving  resistance  and  the 
ultimate  bearing  capacity,  this  approach  acquires  added  sig¬ 
nificance  and  ought  to  be  given  high  priority. 


Practical  experience  in  the  construction  industry  has 
shown  that  pile  driving,  and  especially  also  extracting,  can 
be  done  very  effectively  and  efficiently  by  pile  vibrators. 
Vibratory  pile  driving  was  initiated  in  the  Soviet  Union  in 
the  late  1940s  and  has  found  ever  wider  applications  from 
year  to  year.  A  wide  variety  of  pile  vibrators  has  been 
developed  (Schmid  and  Hill,  1966)  and,  especially  during 
the  last  two  years,  vibrators  powered  by  hydraulic  motors 
have  come  on  the  market.  The  pile  vibrators  operate  on  the 
principle  that,  due  to  the  high  vibrations,  the  skin  fric¬ 
tion  along  the  length  of  the  piles  is  reduced  to  a  negligible 
amount  since  the  soil  particles  adjacent  to  the  pile  are  kept 
in  continuous  motion.  The  remaining  resistance  to  penetra¬ 
tion  is  practically  made  up  entirely  of  the  point  resistance 
which  is  overcome  by  a  high,  dynamic  force  F(t)  given  by 
equation  (5.17),  and  a  dead  load  W  consisting  of  the  weight 
of  the  pile,  pile  vibrator  and  any  bias  surcharge  that  may 
be  added.  If  the  penetration  resistance  of  the  pile  point 
is  known  from  a  relationship  as,  for  example,  given  by  equa¬ 
tion  (5.14),  the  required  dynamic  force  of  the  pile  vibrator 
can  be  determined  from  equation  (5.16). 

There  are  basically  two  different  methods  of  vibrating 
piles:  by  low  frequency-rigid  body  vibrations  and  by  high 

frequency-longitudinal  resonance  vibrations.  The  former  are 
employed  by  most  commercial  pile  vibrators  both  in  the  U.S. 
and  abroad  (Schmid  and  Hill,  1964).  The  resonance  vibration 
principle  is  based  on  a  U.S.  patent  issued  to  A.  Bodine  and 
the  resulting  machine  was  originally  called  a  ’’Sonic  Pile 
Driver.”  The  principle  of  longitudinal  resonance  is  without 
question  applicable  to  driving  piles.  However,  practical 
experience  has  shown  that  the  high  power  necessary  for 
resonance  pile  driving  (the  power  requirements  increase  with 
the  cube  of  the  frequency)  as  well  as  the  exacting  tolerances 
demanded  by  such  a  high  frequency  machine  and  its  maintenance 
problems  make  this  method  impractical  and  uneconomical. 
Especially  since  almost  all  the  advantages  of  vibrodriving 
can  also  be  reasonably  obtained  at  lower  frequencies  and  at 
much  lower  costs.  For  these  reasons  the  development  of  the 
resonant  pile  vibrator  was  plagued  by  numerous  disappoint¬ 
ments  (Jochums  and  Denny,  1961)  and  the  commercial  rights 
and  development  of  the  machine  had  a  checkered  corporate 
history.  The  last  owner  and  operator,  the  Resonant  Pile 
Corporation,  is  now  defunct. 

The  limitations  of  vibratory  pile  driving  are  deter¬ 
mined  by  the  fact  that  the  maximum  available,  dynamic  driv¬ 
ing  force  F(t)  is  limited: 

2 

max  F(t)  =  mQrQO) 


Jt 
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where  iHq  is  the  mass  of  the  eccentrics*  r^  is  their  eccen¬ 
tricity  and  u)  is  their  angular  velocity. 

Since  the  operating  frequency  is  usually  limited,  the 
maximum  driving  force  of  any  pile  vibrator  has  an  upper 
limit.  An  impact  hammer,  in  contrast,  has  a  given  impact 
energy  Eq  available .  In  easy  driving  the  set  per  blow  As 
is  large  and  hence,  the  average  driving  force  P  =  Eq/As  is 
relatively  small.  As  the  pile  "fetches  up"  and  encounters 
heavier  driving  resistance,  the  set  As  becomes  smaller  and, 
thus,  the  available  driving  force  increases.  Therefore,  an 
impact  hammer  has  the  ability  to  furnish  a  higher  punching 
force  in  hard  driving,  albeit  at  a  considerably  reduced 
driving  rate.  A  vibratory  pile  driver,  under  tough  driving 
conditions,  just  stops  and  fails  to  move  the  pile  at  all. 
This  happens  when  the  total  weight  plus  bias,  plus  max  F(t) 
is  smaller  than  the  total  resistance  to  penetration  as  given 
by  equations  (5.8)  or  (5.12).  This  situation  occurs  most 
frequently  in  dense  clays,  hardpan,  glacial  till  or  in  very 
compact  and  cemented  sands  and  gravels. 
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6 .0  Penetration  Instruments  and  Equipnientt  The  State  of 
the  Art 


There  are  a  series  of  activities  that  may  be  considered 
penetration  into  the  ocean  bottom  that  have  been  going  on 
for  some  time  and  it  is  appropriate  to  review  the  current 
state  of  the  art  and  the  available  equipment  for  these  activ¬ 
ities  . 

Probably  the  most  immediately  thought  of  penetration 
activity  is  the  coring  and  sampling  of  ocean  bottom  sedi¬ 
ments.  A  wide  variety  of  new  samplers  and -corers  has  been 
developed  and  put  to  use  recently  with  great  ingenuity.  They 
are  reviewed  and  discussed  below  in  Section  6.1. 

The  anchoring  of  ships,  buoys,  mines  and  fixed  install¬ 
ations  often  requires  the  penetration  of  anchors  or  anchor¬ 
ing  devices  into  the  ocean  bottom.  Although  the  penetration 
of  most  vehicle  anchors  is  quite  shallow,  their  state  of 
development  is  reviewed  briefly  in  Section  6.2, 

One  of  the  most  effective  and  permanent  methods  of 
anchoring  a  structure  is  by  driving  piles  into  the  ocean 
bottom.  This  method  is  extensively  employed  in  fixing  off¬ 
shore  platforms  for  oil  drilling  and  production,  offshore 
light-ships  and  radar  towers,  etc.  Therefore,  a  review  of 
the  current  status  of  the  offshore  pile  driving  technology  is 
also  presented  in  this  Chapter  in  Section  6.3. 

Finally,  there  is  an  interesting  instrument  developed 
recently  that,  while  not  yet  used  extensively  for  ocean  bot¬ 
tom  penetration,  nevertheless  has  been  used  there  occasion¬ 
ally  or  has  the  clear  potential  of  so  being  used.  We  think 
here  of  the  Sandia  Corporation  terradynamic  penetrator.  This 
instrument  is  suitable  either  for  penetration  into  the  ocean 
bottom  as  an  exploratory  tool  or  as  an  anchoring  device. 


6  . 1  Penetration  of  Sediment  Samplers  and  Corers 

A  great  variety  of  sediment  samplers  and  corers  have 
been  developed  in  recent  years  and  a  thorough  review  of  all 
of  them  would  be  beyond  the  scope  of  this  study.  However, 
it  is  obvious  that  the  principles  discussed  in  earlier  chap¬ 
ters  apply  fully  to  the  problem  of  penetration  of  sediment 
samplers  and  corers.  For  example,  the  terminal  velocity  for 
a  free  fall  corer  clearly  can  be  calculated  by  the  method 
demonstrated  in  Section  3.9.  The  depth  of  penetration  can 
be  assessed  by  one  of  the  methods  discussed  in  Section.  4.2 
and/or  Chapter  5,  bearing  in  mind,  however,  that  for  corer 
penetration  the  area  of  penetration  is  that  of  the  projected 
area  of  the  instrument  minus  the  area  of  the  resulting  core. 
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It  is  also  clear  that,  for  submarine  samplers,  the  same 
principles  govern  as  for  ordinary  soil  samplers  (Hvorslev, 
1949),  namely,  that  the  area  ratio  of  the  sampler  Cg  should 
be  reduced'  to  a  minimum.  The  area  ratio  is  defined  as: 


where  is  the  outside  diameter  of  the  sampl'-.r  and  Dg  is 
the  entrance  diameter  of  the  core.  Kvorslev  also  made 
recommendations  regarding  the  inside  clearance  ratio  of 
the  sampler,  stating  that 


should  not  exceed  1.5  per  cent,  Dg  being  the  inside  diameter 
of  the  core  barrel  (see  Figure  6.1a).  The  reason  for  this 
rule  is  that  too  large  a  ratio  would  permit  excessive  expan¬ 
sion  of  the  soil  sample,  while  too  small  a  ratio  would  cause 
increased  disturbance. 

The  problem  of  sample  expansion  may  be  even  more  severe 
for  a  marine  sediment  sample  taken  at  large  water  depth 
because  of  the  expansion  of  the  pore  fluid.  A  detailed 
analysis  of  this  problem,  however,  is  beyond  the  scope  of 
this  study. 

To  increase  depth  of  penetration,  while  keeping  sample 
disturbance  at  a  minimum,  we  feel  that  the  terminal  velocity 
should  be  kept  at  some  optimum  range  of,  say,  15  feet  per 
second  and  any  increased  penetration  be  obtained  at  the 
cost  of  increased  mass  rather  than  increased  velocity.  Also, 
one  of  the  most  critical  parameters  for  good  sample  recovery 
is  the  perfect  alignment  of  the  corer  at  impact.  For  this 
reason  the  marine  sampler  should  be  designed  in  such  a  way 
that  the  vertical  distance  between  the  submerged  center  of 
gravity  and  the  center  of  buoyancy  is  a  maximum. 

Without  claiming  to  present  a  complete  review  of  all 
submarine  sediment  samplers  and  corers,  the  more  representa¬ 
tive  ones  are  reviewed  in  the  following  section. 


6.1.1  The  Gravity  Corer 

This  instrument  is  a  long,  cylindrical  tube  weighted 
by  heavy  blocks  to  provide  a  large  penetration  force  and  a 
high  impact  energy.  A  typical  section  of  an  older  gravity 
corer  is  shown  in  Figure  6.1  as  reported  by  Emery  and  Dietz 
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(1941).  This  corer  was  used  on  the  research  vessel  E .  W. 
Scripps  and  yielded  cores  with  an  average  length  of  6.3 
feet,  and  a  maximum  length  of  over  16  feet.  The  operation 
cf  the  corer  is  described  as  follows: 

"After  the  instrument  is  lowered  into  the  water,  it  is 
dropped  slowly  until  within  about  300  feet  of  the  bottom  as 
indicated  by  the  ship's  echo  sounding  device  and  is  then 
allowed  to  run  rapidly  with  only  a  slight  brake  pressure  to 
insure  a  vertical  descent.  Velocities  of  12  to  21  feet  per 
second  are  attained  just  before  striking  bottom,  regardless 
of  the  length  of  cable  paid  out.  In  depths  shallower  than 
3000  feet  a  slackening  of  cable  indicates  that  the  instru¬ 
ment  has  hit  bottom-  Between  depths  of  3000  and  6000  feet, 
a  pointer  attached  to  an  accumulator  spring  has  been  found 
useful,  but  at  greater  depths  and  in  rough  weather  there  is 
some  difficulty  in  determining  the  instant  of  impact.  When 
bottom  is  reached,  the  winch  is  stopped  immediately  to  pre¬ 
vent  kinking  of  the  cable." 

A  more  modern  gravity  corer  is  the  Sphincter  120  mm 
corer  (Kermabon,  Blavier,  Cortis  and  Delauze,  1966)  which 
seals  the  cores  within  a  watertight  container  during 
retrieval  and  also  has  an  electrical  release  system.  Its 
details  are  shown  in  Figure  6.2. 

A  similar  arrangement  in  which  a  pilot  weight  releases 
the  gravity  corer  to  obtain  a  higher  free  fall  energy  than 
would  be  possible  by  paying  out  the  cable  is  reported  by 
Hvorslev  and  Stetson  (1946)  and  is  shown  in  Figure  6.3. 

Recently,  unattached  sediment  corers  have  been 
developed  (Moore,  1961;  Sachs  and  Raymond,  1965).  The 
operation  of  such  an  unattached  instrument  involves  (i) 
the  launching  of  an  over-ballasted,  deep,  submersible  float 
to  which  the  instrument  is  attached,  (ii)  the  release  of 
the  ballast  at  some  predetermined  depth,  and  (iii)  the 
recovery  of  the  pay-load  at  the  surface.  Particularly  in 
deep  water,  unattached  instruments  can  often  be  operated 
more  effectively  and  with  less  difficulty  than  wire-line 
instruments,  and,  in  some  instances,  they  are  usable  where 
wire-line  instruments  are  not. 

Recently,  Moore  (1961)  solved  the  problem  of  withdraw¬ 
ing,  with  little  force,  a  short  core  from  the  bottom,  and 
this  made  possible  the  successful  development  of  a  free, 
unattached  sediment  sampler. 

The  resulting  sampling  instruments  are  shown  in  Fig¬ 
ures  6.4  and  6.5. 


% 


J 


ir  -r^-Ji^asiXZi^sBasiBmsBaKwsss^^ 


70 


6.1.2  Vibratory  Impact  Coring 

A  quite  recent  development  is  the  coring  of  ocean 
bottom  sediments  by  vibrations.  While  vibratory  soil  samp¬ 
ling  on  land  has  been  going  on  for  about  twenty  years, 
especially  in  the  Soviet  Union,  it  is  only  recently  that 
ocean  bottom  sediments  are  being  sampled  with  the  help  of 
vibrators . 

Figure  6.6  shows  a  vibratory  core  sampler  developed 
by  the  Water  Resources  Division,  U.S.  Geological  Survey, 
Portland,  Oregon  for  a  cooperative  study  with  the  U.S. 
Atomic  Energy  Commission  (Prych  and  Hubbel,  1966).  This 
sampler  has  been  used  successfully  in  water  depths  up  to 
20  meters  and  in  mean  flow  regimes  up  to  1.5  m/sec  in  the 
Columbia  River  estuary.  Ordinarily,  it  is  used  off  the  aft 
end  of  a  12  m  commercial  fishing  vessel  and  achieves  total 
penetrations  of  180  cm  (6  ft.)  by  inducing  longitudinal 
vibrations  with  a  vibrating  time  of  less  than  one  minute. 

Another  vibratory  sampler  of  ocean  bottom  sediments 
that  is  capable  of  applying  longitudinal  as  well  as  tor¬ 
sional  vibrations  is  described  by  Menard  (1965).  Details 
of  his  device  and  method  are  reproduced  in  Figures  6.7  and 
6.8.  Menard  claims  that  this  equipment  is  capable  of  oper¬ 
ating  in  water  depths  up  to  500  fathoms. 

A  list  of  recent  developments  in  ocean  bottom  sampling 
appeared  in  the  1968  Handbook  on  Under  Sea  Technology  and 
is  reproduced  below  as  Table  6.1. 


6.2  Penetration  of  Anchors 


There  are  various  types  of  anchors  available  for  hold¬ 
ing  a  ship  or  object  in  a  relatively  fixed  position.  These 
anchors  and  several  new  anchor  types  developed  by  the  Naval 
Facilities  Engineering  Command  were  discussed  by  Tudor 
(1967).  A  short  discussion  of  them  is  presented  below. 


6.2.1  Gravity  or  Dead  Weight  Anchors 

These  are  the  most  elemental  type  of  anchors.  They 
are  massive  and  have  a  low  holding  power  to  weight  ratio. 
Because  of  their  mass  they  are  'expensive  and  difficult  to 
handle  and  transport.  They  also  have  a  tendency  to  move 
laterally  at  small  lateral  loads  and  are  not  very  reliable 
if  the  ocean  bottom  is  sloping  or  unstable.  Their  penetra¬ 
tion  characteristics  are  negligible. 
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Table  6.1. 
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6 .2 . 2  The  Standard  Drag-Type  or  Fluke-Type  Anchors 


These  anchors,  such  as  the  Navy  standard,  stockless 
anchor,  the  Admiralty  anchor,  or  the  Danforth  anchors,  pro¬ 
vide  primarily  a  horizontal  holding  power.  Large,  flat 
flukes  penetrate  into  the  ocean  bottom  by  being  dragged 
along  the  bottom  and  gain  their  holding  power  from  the 
shear  strength  of  the  bottom  sediment.  Details  of  such  a 
f-luke-type  anchor  are  shown  in  Figure  6.9  which  shows  the 
BUDOCKS  STATO  mooring  anchor  (Towne,  1959).  Superior  hold¬ 
ing  power  in  mud  bottom  is  claimed  for  this  anchor.  The 
penetration  of  these  anchors  is  limited  to  the  fluke  length. 
Often,  a  large,  horizontal  distance  is  needed  to  set  this 
kind  of  anchor  adequately,  requiring  large  amounts  of  line 
and  gear.  Typical  ratios  H/d  of  length  of  anchor  cable  to 
depth  of  water  are  2.0  to  3.0. 


6.2,3  Explosive  Anchors 

Explosive  anchors  stem  from  a  relatively  recent  and 
intensive  development.  They  are  small,  lightweight,  rapidly 
installed,  and  are  adaptable  to  a  variety  of  ocean  bottom 
conditions.  Development  has  been  carried  out  primarily  by 
the  Concept  Engineering  Company,  by  Aerojet-General  Corpor¬ 
ation,  and  by  the  U.S.  Army  Engineering  Research  and  Develop¬ 
ment  Laboratory  (ERDL).  Explosive  anchors  are  currently 
being  used  in  shallow  depths  of  200-300  feet  but  there  is 
no  obstacle  that  would  preclude  their  use  in  depths  to  1000 
feet  or  more. 

Figure  6.10a  shows  an  explosive  embedment  anchor 
developed  by  the  U.S.  Army  Engineers  Research  and  Develop¬ 
ment  Laboratory  (ERDL).  This  anchor,  according  to  Tudor, 
is  being  modified  by  NCEL  for  use  in  water  depths  up  to 
6000  feet.  Figure  6.10b  shows  a  cross-section  of  another 
explosive  self-embedment  anchor  developed  by  the  Magnavox 
Company.  It  consists  of  a  flotation  assembly,  electronics 
package,  the  anchor,  and  cable  pack.  When  the  system 
enters  water,  the  assemblies  separate  and  the  cable  pays 
out  from  the  back  of  the  anchor  package  during  descent, 
preventing  cable  drag.  When  the  anchor  package  hits  bot¬ 
tom,  the  unit  fires,  driving  the  anchor  into  the  ocean 
floor.  At  the  same  time  the  cable  is  locked  up,  prevent¬ 
ing  further  payout.  The  anchor  is  tied  to  its  package 
housing  by  a  cable  long  enough  to  ensure  against  breakage 
due  to  premature  anchor  snubbing.  As  the  cables  pull  taut, 
the  spring-loaded  anchor  flukes  dig  in  and  spread  o.ut . 

System  safety  is  assured  by  design.  The  fuse  is  armed  by 
hydrostatic  pressure,  and  the  arming  mechanism  is  locked 
until  the  system  enters  water  and  descends  to  a  preset 
depth.  The  use  of  mechanical  devices  also  eliminates 
unwanted  initiation  by  stray  electromagnetic  signals. 
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For  all  explosive  anchors*  the  explosive  charge  must 
be  properly  designed  so  that  the  required  penetration  is 
achieved.  If  the  driving  energy  is  too  little,  the  desired 
penetration  and,  hence,  holding  power  will  not  be  attained. 

If  it  is  grossly  in  excess,  the  anchor  may  be  damaged.  The 
required  explosive  energy  can  be  calculated  for  any  given 
penetration  from  a  penetration  formula,  say,  equation  (4.26): 


req 


=  P 


X  =  A  p 
req 


V 

o  "^req 


(6.1) 


where  A  is  the  cross-sectional  area,  p^  the  bearing,  capacity 
and  is  the  required  depth  of  penetration  determined  by 

the  desired  pull-out  force. 

Tudor  reports  recorded  holding  power  for  the  ERDL 
explosive  anchor  in  sand  and  clay  of  60,000  lbs.  for  the 
50  kip  rated  anchor  and  of  200,000  lbs.  for  the  anchor 
rated  at  300  kip.  In  soft  mud  these  two  anchors  held 
15,000  lbs.  and  58,000  lbs.  respectively. 


6.2.4  Free-Fall  Embedment  Anchors 


Like  the  free-fall  corer  or  sampler,  this  anchor 
achieves  its  penetration  by  the  energy  of  its  impact  .when 
hitting  the  bottom.  Thus,  the  determination  of  the  .strik¬ 
ing  velocity  must  follow  the  principles  and  methods  described 
in  Chapter  3  and  the  penetration  phenomenon  is  controlled  by 
the  relationships  discussed  in  Chapter  4. 

One  such  typical  free-fall  anchor  and  its  placement  is 
shown  and  discussed  by  Smith  and  Dantz  (1963).  Tudor  (1967) 
also  shows  a  free-fall  embedment  anchor  with  a  holding  power 
to  weight  ratio  of  6:1.  Terminal  velocities  obtainable 
appear  to  be  in  the  range  of  20  to  25  fps. 


6.2.5  Umbrella  Pile  Anchors 

These  anchors  were  designed  to  be  driven  or  vibrated 
into  the  sea  bottom  like  a  pile  and  then  open  up  to  obtain 
greater  bearing  capacity  or  holding  power  and/or  greater 
resistance  to  uplift  or  break-out.  A  typical  example  is 
shown  by  Tudor  (1967)  who  also  gives  bearing  or  holding 
capacities  that  can  be  achieved  with  this  anchor.  Umbrella 
pile  anchors  are  used  where  conventional  anchors  cannot  be 
dragged,  c  do  not  provide  sufficient  embedment,  or  where 
the  bedrock  is  too  deep  for  obtaining  point  support  of  bear¬ 
ing  piles.  The  Vertohold  Embedment  Anchor  of  Edo-Western 
Co.  is  an  umbrella  type  anchor  that  is  driven  into  the  bot¬ 
tom  soil  by  an  explosive  charge  fired  either  upon  bottom 
contact  or  by  remote  control. 


74 


Holding  capacities  of  220  kips  in  bearing  and  of  300 
kips  in  uplift  have  been  reported  for  such  umbrella  type 
anchors  at  an  embedment  of  20  ft.  in  sand.  When  embedded 
in  blue  mud  interspersed  with  fine  sand  and  shells,  a 
capacity  in  excess  of  100  kips  has  been  recorded.  In 
general,  the  holding  power  of  such  umbrella  pile  anchors 
should  be  proportional  to  the  projected  area  of  the  opened 
umbrella  and  to  the  square  of  the  depth  of  embedment . 


6.2.6  U.S.  Patents  on  Sea  Bottom  Anchors 


The  problem  of  anchoring  a  device  in  the  ocean  bot¬ 
tom  has  attracted  the  interest  and  imagination  of  a  wide 
variety  of  scientists,  engineers  and  inventors,  as  a  look 
at‘  the  list  of  U.S.  patents  given  in  Table  6.2  will  reveal: 


Table  6.2 


List  of 

Recent  U.S.  Patents 

on  Sea 

Bottom  Anchors 

U.S.  Patent  # 

Subject 

Inventor 

Date 
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Thomason 

10/27/64 

et 

al 

• 

3,170,433 

Imbedment  Anchor 

J. 

T. 

Gardiner 

2/23/65 

Assembly 

Re.  25,764 

Embedment  Anchor 

R. 

S. 

Robinson 

4/13/65 

3,187,705 

Dynamic  Anchor 

R. 

B. 

Costello 

6/8/65 

et 

al 

• 

3,207,115 

Explosive -Ope rated 

M. 

H. 

Anderson 

9/21/65 

Anchor  Assembly 

3,280,782 

Marine  Anchor 

C. 

H. 

Holm 

10/25/66 

3,291,092 

Mooring  Apparatus 

P. 

V  . 

Halberg 

12/13/66 

et 

al 

• 

3,311,080 

Pressure  Actuated 

V. 

C. 

Anderson 

3/28/67 

et 

al 

. 

3,315  ,629 

Underwater  Anchor 

J. 

A. 

Standridge 

4/25/67 

Gun  Device 

et 

al 
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Table  6.2  (continued) 

U.S.  Patent  §  Subject 

Inventor 

Date 

3,330,338 

Anchor  and  Method 
of  Installing 

J. 

R. 

Dozier 

7/11/67 

3,370,566 

Embedment  Device 

J  . 

A  . 

Dorr  et 

al . 

2/27/68 

3,371,643 

Hydraulically 
Actuated  Driver 

W. 

H. 

Dunham 

3/5/68 

3,399,646 

Submarine  Anchor 
Assembly 

R  . 

P  . 

Vincent 

9/3/68 

3,411,473 

Deep  Water  Anchor 

6. 

E  . 

Mott  et 

al. 

•l<l/19/68 

3,417,724 

Vibratory  Sea 
Anchor  Driver 

T. 

B  . 

Edwards 

12/24/68 

Nearly  all  of  these  patents  concern  mechanisms  and 
details  for  explosive  anchors  describing  and  claiming 
specific  features  for  arming,  igniting,  and  achieving  pene¬ 
tration  and  anchorage  for  these  devices. 


6 . 3  Penetration  of  Piles  into  the  Ocean  Bottom 

Piles  have  been  driven  into  the  ocean  bottom  for  hun¬ 
dreds  of  years.  These  installations  usually  were  for  near¬ 
shore  facilities  such  as  docks,  wharfs,  quays  and  piers,  and 
the  piles  are  usually  driven  from  above  the  water. 

Within  the  last  two  or  three  decades,  however,  an  entire 
new  industry  has  developed  concerned  with  the  exploration 
and  production  of  natural  oil  and  gas  far  offshore  on  the 
continental  shelves  .  Most  of  this  activity  takes  place  on 
platforms  anchored  to  the  ocean  bottom  by  long,  penetrating 
piles  and  some  of  the  technical  developments  of  this  new 
industry  are  of  interest  here. 


6.3.1  Offshore  Construction  and  Its  Technoloj 


Beginning  in  rather  shallow  water  offshore  Louisiana 
shortly  before  and  after  the  Second  World  War,  offshore  drill¬ 
ing  and  oil  production  progressed  at  a  dizzying  pace  such  that 
by  the  end  of  1957  approximately  1200  producing  oil  and  gas 
wells  and  500  nonproductive  wells  had  been  drilled  on  the. 
northern  shelf  of  the  Gulf  of  Mexico.  Since  then  "offshore 
fever"  has  spread  all  over  the  globe  and  by  1969  offshore  oil 
exploration  and  production  is  carried  out  among  other  places 
off  the  coasts  of  California,  Alaska,  West  Africa,  Australia, 
Indonesia,  Brazil,  Nicaragua,  Peru,  Trinidad,  Nova  Scotia, 
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Spain,  France,,  Tunisia,  New  Zealand,  Borneo,  Egypt,  and  in 
the  Arabian  Gulf,  the  Gulf  of  Suez,  the  Sea  of  Japan,  the 
Adriatic,  the  Java  Sea,  the  Nigerian  Gulf,  the  Caspian  Sea, 
Lake  Maracaibo,  Lake  Erie  and  in  the  North  Sea.  By  1967, 
after  20  years  of  operation,  well  over  2000  platforms  had 
been  installed  and  over  9000  wells  had  been  drilled. 

Expenditures  in  the  offshore  industry  at  present  are 
in  excess  of  one  billion  dollars  per  year  and  the  world 
petroleum  industry  is  expected  to  invest  at  least  25  bil¬ 
lion  dollars  off  the  coasts  of  more  than  60  nations  during 
the  next  decade  to  explore,  develop,  and  produce  offshore 
oil  fields. 

At  a  cost  of  several  million  dollars  apiece,  the  cur¬ 
rent  inventory  of  the  industry  includes  150  mobile  offshore 
drilling  rigs  and  120  self-contained  drilling  platforms. 
Table  6.3  shows  the  extent  of  recent,  new  construction  of 
mobil  rigs  supported  on  piles.  In  addition,  during  the  same 
period  covered  by  the  Table,  twelve  drill  ships  and  ten 
semi-submersibles  were  under  construction.  Drillships  and 
submersibles  do  hot  utilize  piles  that  penetrate  into  the 
ocean  bottom. 

With  all  this  activity  in  ocean  bottom  pile  driving, 
it  is  proper  and  necessary  in  the  context  of  this  study  to 
review  the  state  of  the  art  of  pile  driving  and  penetration 
into  the  ocean  bottom  as  currently  practiced  by  the  offshore 
industry.  Before  doing  this  it  may  be  useful  to  present 
some  general  background  information  on  offshore  platforms 
and  structures  and  the  associated  engineering  problems. 

The  current  technology  of  offshore  structures  has 
limited  the  exploration  and  exploitation  of  reserves  to 
the  relatively  shallow  depths  of  the  continental  shelf. 

Even  so,  the  structural  design  of  production  platforms  has 
imposed  great  challenges  on  the  engineer. 


6.3.2  Types  of  Offshore  Structures 

The  basic  objective  of  the  offshore  structure  is  to 
drill  and  support  several  oil  wells  in  the  ocean  floor.  To 
meet  this  end,  several  different  types  of  structures  are 
available.  The  selection  of  one  structure  over  another 
would  necessitate  a  performance  comparison  of  at  least  the 
following  items; 

(1)  Stability  on  location 

(2)  Ease  of  movement 

(3)  Minimum  dependence  on  other  craft 


Table  6.3  Mobil  Rig  Construction  for  13  Months  Ending  June  1968 

Water 

Depth 

Rig  Type _ Name  Owner  Builder  Completion  Capacity 


Propulsion- 

Assisted 

J  ackup 

Rig  58 

Penrod 
Drilling 
Co . 

Propulsion- 

Assisted 

Jackup 

Rig  59 

Penrod 
Drilling 
Co . 

Self- 

Elevating 

Drillship 

Un¬ 

named 

The 

Offshore 
Co . 

Jackup 

Maersk 

Explor¬ 

er 

Zapata 

Offshore 

Jackup 

Rig  59 

Trans¬ 

world 

Drlg  Co. 

Jackup 

Gulf- 

tide 

Husky 

Oil 

( Alberta) 

Jackup 

Chazar 

USSR 

Jackup 

Petro- 
bras  I 

Petro- 

bras 

Jackup 

Sed- 
neth  II 

Sea  Drig 
Nether¬ 
lands 

Jackup 

Norsmec 

I 

Norsmec 

Jackup 

Un¬ 

named 

Japan 

Ocean 

Drlg . 

Workover- 

Jackup 

Dresser 
Rig  I 

Dresser 

Offshore 

Workover- 

Jackup 

Dresser 
Rig  II 

Dresser 

Offshore 

Workover- 

Jackup 

Dresser 
Rig  III 

Dresser 

Offshore 

LeTourneau 

Late  1968 

200 

LeTourneau 

1969 

200 

Upper 

Clyde  Ship¬ 
builders  , 
Ltd. 

Winter 

1969 

250 

LeTourneau 

May  1967 

125 

Bethlehem 

May  1967 

125 

John  Brown 

6  Co . 

( Clydebank ) 
Ltd. 

August 

1967 

200 

IHC  Holland 

August 

1967 

200 

Companhia 

Commercio 

E  .Navagacao 

Imminent 

100 

Rotterdam 

Dockyards 

May  1968 

200 

Norsmec 

June  1968 

150 

Seeking 

Contract 

McDermott 

Shipyard 

May  1967 

60 

McDermott 

Shipyard 

May  1967 

60 

McDermott 

Shipyard 

May  1968 

70 

78 


Table  6.3  (continued) 


Rig  Type 

Name 

Owner 

Builder 

Completion 

Water 

Depth 

Capacity 

Workover- 

Jackup 

Dresser 
Rig  IV 

Dresser 

Offshore 

McDermott 

Shipyard 

Imminent 

70 

Workover- 

Jackup 

Topper 

II 

Crestwave 

Offshore 

LeTourneau 

March  1968 

125 

Workover- 

Jackup 

Cowrie 

I 

Shell  (U. 
K.)  Ltd. 

IHC 

Holland 

May  1968 

30 

Workover- 

Jackup 

Ranger 

I 

Walker- 

Huthnance 

Bethlehem 

Imminent 

70 

(4)  Towability 

(5)  Ease  of  equipment  and  personnel  transfer 
under  adverse  sea  conditions 

(6)  Efficiency  of  drilling  operation 

(7)  Storage  capacity 

The  structures  available  may  be  grouped  as  follows: 

(1)  Mobile  platforms 

(2)  Fixed  platforms 

Mobile  platforms,  possessing  the  ability  to  move  to 
different  locations,  were  conceived  as  an  exploratory  tool. 
Fixed  platforms  would  subsequently  develop  those  areas  which 
were  found  to  contain  sufficient  petroleum  reserves  to  justify 
permanent  installations.  In  practice,  the  mobile  platforms 
have  been  found  capable  of  economical  production  of  complete 
oil  fields.  Both  types  will  be  described  briefly. 


6.3.3  Mobi le . Platforms 


Jack-Ups 

This  structure  is  the  most  widely  used  mobile  platform 
in  use  today.  In  concept,  the  structure,  consisting  of  a 
hull  and  three  or  more  legs,  is  towed  out  to  sea,  and  the 
legs  are  jacked  down  far  enough  to  raise  the  hull  out  of  the 
water  a  predetermined  amount.  (See  Figures  6.11,  6.12,  and 
6.13.)  The  structure  must  fulfill  the  following  functions: 


(1)  Position  of  the  platform  must  not  be  unduly 
influenced  by  the  action  of  the  sea  and  wind. 


(2)  Platform  must  be  seaworthy. 

(3)  Raising  and  lowering  of  legs  must  be  safe 
and  efficient. 

(4)  Withdrawal  of  legs  which  have  penetrated 
the  sea  bottom  must  be  possible  without 
endangering  the  structure. 

The  jack-up  legs  may  be  circular,  triangular . or  square  in 
cross-section  and  are  usually  of  tubular  design.  The  bot¬ 
tom  of  the  legs  are  usually  fitted  with  a  watertight  can  in 
order  to  increase  the  bearing  area  and  prevent  excessive 
penetration  into  the  ocean  floor.  Difference  in  bottom 
conditions  will  determine  the  shape  of  this  can.  Water 
jets  are  often  used  to  settle  the  legs  in  the  ocean  floor 
or  free  a  wedged  leg  before  rig  transfer. 


Submersible  Jack-Ups 

This  is  the  further  refinement  of  the  jack-up  for 
operation  in  soft  bottom  conditions.  The  jack-up  legs  are 
attached  to  a  large,  submerged  barge  and  do  not  penetrate 
into  the  ocean  bottom.  (See  Figures  6.14  and  6.15.)  The 
submerged  hull  can  permit  a  fixed  condition  to  be  assumed 
for  the  column  legs. 


Semi-submersibles 


Semi-submersibles  can  work  as  bottom-supported  units 
or  floaters.  (See  Figures  6.16  and  6.17.)  The  bottom  hull 
is  attached  to  the  working  hull  by  large  diameter,  tubular 
columns.  Stability  is  maintained  by  the  fact  that  the 
center  of  gravity  is  below  the  surface  of  the  water.  Posi¬ 
tion  may  be  maintained  by  anchors.  The  rig  is  extremely 
difficult  to  tow,  but  is  capable  of  production  in  1500-foot 
water  depths  . 

The  area  requirements  for  the  bottom  can  of  the  jack¬ 
up  legs  or  the  submerged  hull  of  a  submersible  jack-up  - 
once  the  strength  of  the  bottom  sediment  is  known  or  esti¬ 
mated  -  can  be  calculated  by  the  methods  of  Section  4.2.1 
and  4.2.2. 


6.3.4  Fixed  Platforms 


Fixed  platforms  (see  Figures  6.18 
of  three  main  parts:  template  (jacket), 
The  template  is  fabricated  from  trussed, 
a  shipyard  and  towed  to  sea.  At  arrival 


and  6.19)  consist 
piles,  and  hull, 
tubular  steel  at 
on  location,  the 


jacket  is  lifted  or  launched  in  the  water  with  its  hollow 
legs  upright.  The  structure  is  then  leveled.  Piles  (con¬ 
sisting  of  large  diameter  steel  tubes)  are  driven  through 
the  template  to  the  desired  depth.  The  tops  of  the  piles 
are  subsequently  welded  to  the  template.  Finally,  the  hull, 
consisting  of  trussed  or  plate  girders,  is  welded  to  the  top 
of  the  piles.  The  disadvantage  of  this  unit  is  the  lack  of 
mobility.  Platforms  have  been  designed  to  work  in  340-foot 
depths  and,  at  present,  a  design  for  a  platform  in  400  feet 
of  water  has  been  completed  and  is  in  the  process  of  bidding. 


Structural  Evolution  of  the  Fixed  Platform 


Current  versions  of  fixed  platforms  have  been  pro¬ 
foundly  influenced  by  structural  improvement  of  the  earliest 
designs.  At  the  beginning  of  the  oil  industry's  quest  for 
underwater  deposits,  the  structures  were  simple.  Rigs  were 
supported  by  200-300  timber  piles  of  a  maximum  length  of 
about  90  feet.  The  decks  were  large  enough  to  accommodate 
the  normal,  land-based  drilling  apparatus.  Structural  design 
was  limited  to  a  normal  design  for  the  weight  and  configura¬ 
tion  of  the  deck  versus  the  bearing  capacity  of  the  piles. 

The  effects  of  wind  and  waves  were  acknowledged  by  a  liberal 
amount  of  bracing  above  the  water  line,  but  insufficient  know 
ledge  of  these  forces  prevented  a  rigorous  analysis. 

In  retrospect  we  must  admire  the  offshore  pioneer  and 
his  timber  pile  rig.  Admittedly,  many  losses  did  occur,  but 
pre-war  daring  provided  a  basis  for  post-war  burgeoning.  In 
1946,  steel  was  substituted  for  timber  and  the  size  of  the 
deck  was  reduced  to  a  minimum.  The  drilling  operation  was 
supported  by  a  tender  which  was  anchored  near  the  rig.  Acces 
was  maintained  between  the  rig  and  the  tender  by  a  ramp  or 
"widow  maker."  During  heavy  seas,  operations  had  to  be  sus¬ 
pended  while  the  tender  backed  away  from  the  drilling  plat¬ 
form.  This  discontinuity  in  operation  led  by  evolutionary 
change  to  the  "self-contained  platform."  (The  tender  oper¬ 
ation  has  not  been  discontinued.) 

This  concept  placed  all  of  the  operation  back  on  one 
deck,  and  forced  the  structural  designer  to  more  sophisti¬ 
cated  methods  of  analysis.  Proceeding  with  a  minimum  of 
data,  an  imperfect  knowledge  of  the  theory  of  wave  motion, 
and  a  few  newspaper  accounts  of  the  magnitude  of  hurricanes, 
the  early  structures  were  designed  and  built.  According  to 
the  "eternal  amazement"  of  one  authority,  some  of  them  are 
still  in  use  after  more  than  20  years  of  service  . 

The  configuration  and  size  of  the  structure,  the  choice 
of  structural  materials  and  the  varied  type  of  drilling 
requirements  resulted  in  manj'  different  designs.  The  final 


approach,  as  noted  in  th^  description  of  fixed  platforms, 
was  to  build  a  template  out  of  tubular  steel,  float  this 
template  out  to  sea,  sink  it  to  the  boi'+pm  at  the  desired 
location  and  drive  piles  through  the  tubes,.  The  original 
designs  called  for  many  tubes  -  30  to  40  -  about  12  to  16 
inches  in  diameter  and  12  to  15  feet  apart.  After  the  pil¬ 
ing  had  been  set  in  place,  the  structure  was  placed  on  top 
and  drilling  could  begin.  The  relatively  close  spacing  of 
these  piles  i^rmed  a  wall  which  storm  waves  would  strike 
with  enormous  impact,  overturning  moment,  and  uplift.  Many 
structures  based  on  this  concept  were  lost.  The  idea  was 
further  refined  by  decreasing  the  number  of  tubes  in  the 
template  and  increasing  their  diameter.  Consequently, 
modernized  designs  called  for  eight  to  ten  30-inch  piles 
of  about  300-ton  capacity  apiece  at  a  spacing  of  36>  feet. 
The  fewer  number  of  piles  allowed  waves  to  pass  through  the 
structure  with  far  less  wave  impedance.  This,  in  turn^ 
meant  that  trussed  spans  were  required  to  carry  the  deck 
loads  between  the  fewer  supporting  piles. 

Of  course,  bigger  piles  meant  bigger  hammers,  leads, 
cranes,  fabrication  techniques  and  barges.  The  require¬ 
ments  of  design  spawned  a  new,  heavy  industry. 


6.3.5  Pile  Driving  in  Offshore  Construction 

The  evolution  of  offshore  platform  design  and  con¬ 
struction  was  reviewed  recently  by  Lee  (1968).  Figures 
6.20  and  6.21  show  the  development  in  the  design  from  1947 
to  1967  and  the  records  of  water  depth  for  fixed  platform 
construction . 

What  this  development  has  meant  for  the  problem  of 
driving  the  piles  was  discussed  and  shown  by  McClelland, 
Focht  and  Emrich  (1967,  Figures  6.22  through  6.25). 

Piles  presently  used  in  offshore  construction  are 
typically  30  to  48-inch  diameter.  Penetrations  to  produce 
the  required  ultimate  capacities  are  frequently  in  the 
range  of  250  to  350  feet  -  although  these  vary  widely. 

Since  the  piles  generally  extend  some  distance  above  the 
water  level,  total  pile  lengths  of  400  to  600  feet  occur 
frequently . 

McClelland  et  al.  show  some  typical  pile  dimensions 
(Figure  6.22)  and  the  range  of  pile  loads  (Figure  6.23). 

The  installation  procedure  is  shown  in  Figure  6.24. 
When  considering  the  problem  of  offshore  platform  construc¬ 
tion  one  must  bear  in  mind  that  the  installation  of  a  pile, 
above  or  below  water,  not  only  consists  of  driving  the  pile 
home.  Before  this  can  be  done  the  pile  has  to  be  properly 
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positioned  and  has  to  be  guided  throughout  the  driving  pro¬ 
cess.  Also,  the  pile  driver  has  to  rest  on  or  be  attached 
or  clamped  to  the  pile.  It  is  for  these  reasons  that  we 
foresee  future  underwater  construction  in  water  depths  up 
to,  perhaps,  1000  to  2000  feet  to  follow  very  much  the 
technology  developed  by  the  offshore  platform  construction 
industry  and  suggested  by  Figures  6.24A  through  C: 

(1)  A  template  towed  (horizontally)  to  location 
and  positioned  by  properly  sequenced  flooding 
of  selective  compartments . 

(2)  The  driving  of  piles  through  the  template 
sleeves  by  pile  drivers  that  follow  into  the 
sleeves  below  the  water. 

(3)  Connecting  the  driven  piles' and  the  lower  part 
of  the  template  by  grouting  the  pile  into  the 
sleeve  . 

(4)  Cutting  off  the  upper  section  of  the  template, 
refloating  it  and  returning  it  to  port  for 
reuse.  The  bottom  part  of  the  template  either 
would  completely  contain  the  fixed  structure  to 
be  installed  on  the  ocean  bottom,  or  it  would 
have  a  template  device  for  easily  docking  the 
facility  onto  the  pile  trestle  resulting  from 
this  construction. 

One  of  the  major  problems  in  offshore  pile  penetration 
has  been  the  availability,  or  rather  the  lack  of  the  neces¬ 
sary  pile  hammers.  The  difference  between  a  normal  pile 
driving  job  on  land  and  a  typical  offshore  job  is  enormous 
(see  for  example.  Figure  6.25  comparing  the  load  range  of 
conventional  pile  load  tests  on  land  to  a  typical  offshore 
pile  load).  The  biggest  pile  hammers  that  were  on  the 
market  until  just  a  few  years  ago  had  a  maximum  energy  of 
about  60,000  ft .  lbs  . /blow .  •*' 

The  past  few  years,  however,  have  seen  the  advent  of 
significantly  heavier  equipment.  Delmag  has  introduced  a 
new  D-44  Diesel  pile  hammer  with  a  maximum  energy  of  87,000 
ft. lbs. /blow  (Figure  6.26).  The  Vulcan  Ironworks  have 
developed  the  040  single  acting  steam  hammer  with  a  ram 
weight  of  40,000  lb.,  a  drop  height  of  3  ft.,  and,  hence,  a 
maximum  energy  of  120,000  lbs.  (Figure  6.27).  And  the  MKT 
Corporation  has  developed  a  special  offshore  single  acting 
steam  hammer,  the  OS-60,  with  a  ram  weight  of  60,000  lb.,  a 
fall  of  3  ft.,  and,  thus,  a  total  energy  of  180,000  ft. lbs/ 
blow  (Figure  6.28). 

The  most  important  new  development,  however,  for  pile 
driving  and  penetration  into  the  ocean  bottom,  in  our 
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opinion,  appears  to  be  the  development  of  a  hydraulically- 
powered  pile  vibrator. recently  introduced  by  .the  MKT  Division 
of  Koehring  Company,  Dover,  Mew  Jersey:  the  MKT  V-10  (Figure 
6.29).  This  vibratory  pile  driver  has  a  maximum  dynamic 
force  of  112,000  lbs.,  a  total  weight  of  9,500  lbs.,  and 
operates  at  frequencies  from  1500  to  1850  cpm  at  a  total 
power  of  110  h.p.  It  is  sealed  against  submersion  .in  water 
although  operation  at  large  water  depth  undoubtedly  will 
require  special  gaskets  and  packing  seals,  etc. 

The  adv:antage  '  of  hydraulic  power  for  an  underwater 
pile  driver  is  that  the  prime  mover  (e.g.  electric  or  diesel 
motor  and  hydraulic  pump)  can  remain  above  the  water  on  a 
barge  or  ship  and  the  driving  tool,  powered  through  a  hydrau¬ 
lic  hose  umbilical  cord,  is  driven  by  hydraulic  motors.  This 
is  a  closed  system  requiring  no  exhaust  under  water.  Such  a 
system  does  not  have  the  hazard  nor  the  sealing-connection 
problems  associated  with  high  power  electric  cables.  Also, 
since  pressure  in  the  hydraulic  hose  as  well  as  the  ambient 
pressure  in  the  ocean  increases  linearly  with  water  depth 
and,  in  both  cases,  almost  at  the  same  rate,  the  pressure 
difference  and,  hence,  the  circumferential  stress  in  the  hose- 
remains  practically  constant  throughout  the  entire  depth. 

The  only  depth  limitations,  in  principle,  appear  to  be  the 
friction  losses  in  the  hydraulic  hoses  between  the  driving 
tool  and  the  prime  mover. 

As  the  water  depths  in  which  piles  have  to  be  driven 
continue  to  increase,  it  is  quite  clear  that  sooner  or  later 
piles  will  have  to  be  driven  not  from  a  platform  or  from  a 
template  above  water  but  from  below  the  water  surface.  The 
method  of  providing  the  required  power  for  the  prime’  mover 
for  such  an  underwater  job  then  becomes  critical. 

A  limited  underwater  driving  capability  exists  today 
with  existing  equipment  by  the  use  of  an  MKT  single  acting 
pile  driver  that  can  be  modified  for  underwater  operations . 
This  type  of  pile  driver  is  capable  of  driving  in  water 
depths  of  150  to  200  feet.  The  deepest  operation  of  such 
an  underwater  pile  driver  so  far  was  accomplished  during 
the  construction  of  the  pier  foundations  for  the  Narragan-  . 
sett  Bay  Bridge  in  Newport,  Rhode  Island  where  12-inch 
steel  H-piles  were  driven  at  water  depths  of  up  to  162  feet. 
The  operation  of  the  pile  driver  at  these  depths  required 
compressed  air  at  125  psi  and  2400  cfm.  In  most  cases, 
however,  where  piles  are  driven  into  the  ocean  bottom,  driv¬ 
ing  still  is  carried  out  above  water  level. 

There  are  a  number  of  possible  alternatives  for  power¬ 
ing  an  underwater  pile  driving  operation  (see,  for  example. 
Publication  No.  1702,  National  Academy  of  Sciences)  and  the 
possibilities  and  limitations  should  be  explored  by  a 
detailed  feasibility  and  preliminary  design  study  particu¬ 
larly  if  the  present  state  of  the  art  is  to  be  advanced 
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significantly.  But  we  believe  that  a  hydraulically-powered 
vibrator  (or  impact  hammer)  is  uniquely  suited  for  under¬ 
water  application  where  the  prime  mover  remains  above  water 
and  power  is.  supplied  by  hydraulic  hoses.  A  pile  driving 
tool  such  as  the  V-10  should  be  able  to  operate  easily  with 
but  minor  modifications  in  water  depths  up  to  1000  feet. 

The  only  limitation  for  extending  operational  capabilities 
to  even  larger  depths  appears  to  be  the  friction  losses  in 
the  hydraulic  fluid  hose  lines  between  prime  mover  and 
driving  tool.  Clearly,  at  some  water  depth  of,  say,  2000 
feet  or  more  it  will  become  advantageous  to  dispense  with 
any  umbilical  cord  and  develop  an  integral  system  of  prime 
mover  and  driving  tool  both  operating  under  water.  The 
development  costs  for  such  a  system,  however,  will  be  con¬ 
siderable.  It  will  certainly  be  many  multiples  of  the  costs 
of  adapting  a  hydraulic  vibrator  for  underwater  operation. 

As  explained  in  Section  5.9,  pile  vibrators  have,  as 
of  this  date  (1969),  definite  limitations  in  driving  capa¬ 
city  in  stiff,  dense  clays,  hardpan,  and  very  compact  or 
cemented  sands  and  gravel.  Vibrators  normally  outperform 
impact  hammers  under  medium  and  light  driving  conditions  by 
increasing  driving  speeds  by  a  factor  of  at  least  ten.  If 
soil  conditions  prevail,  however,  as  those  cited  above, 
vibrators  may  not  be  able  to  drive  piles  at  all.  This  is 
due  to  the  still  relatively  low  power  tools  that  are  avail¬ 
able.  Maximum  power  ratings  of  most  commercial  vibrator 
lines  on  the  market  are  in  the  100  to  140  h.p.  range.  How¬ 
ever,  second  generation  machines  are  under  development  that 
will  significantly  increase  the  available  power  and,  hence, 
the  available  driving  force . 


6.4  The  Sandia  Penetrator 


An  interesting  new  arrival  under  exploration  already 
for  a  number  of  years  is  the  earth-penetrating  projectile 
developed  by  Sandia  Corporation  (Colp,  1965,  1967,  1968; 
Thompson  et  al.,  1965;  Thompson,  1966;  Young,  1969).  A 
section  of  a  typical  penetrator  is  shown  in  Figure  6.30. 

This  tool  is  essentially  a  solid  steel  cylinder  with  a 
conical  or  ogive  nose  that  is  given  very  high,  up  to  ball¬ 
istic,  vertical  velocities  by  dropping  it  at  large  heights 
from  an  airplane  or  helicopter.  Near  vertical  position  at 
release  and  aerodynamic  fins  assure  nearly  vertical  align¬ 
ment  at  impact.  The  afterbody  of  the  penetrator  holds  a 
piezo-electric  accelerometer  giving  the  deceleration  of  the 
instrument  through  a-  trailing  wire.  The  device  has  a 
potential  both  as  a  penetrometer  measuring  the  deceleration 
experienced  as  it  traverses  soil  layers  of  various  strengths 
and  depth  as  well  as  for  the  embedment  of  an  anchor.  The 
use  of  the  device  as  a  soil  exploration  instrument  is  demon¬ 
strated  by  Figures  6.31  through  6.34  which  are  taken  from 


various  reports  by  Colp  and  Thompson  cited  above .  In  par¬ 
ticular,  Figure  6.32  shows  how  the  instrument  might  be  used 
for  exploring  the  properties  of  ocean  bottom  sediments.  It 
appears  that  the  method,  however,  is  rather  expensive, 
requiring  the  use  of  an  aircraft  or  helicopter  and  entailing 
the  loss  of  the  penetrometer  for  every  point  where  the 
properties  are  measured.  Also,  it  does  not  lend  itself  for 
use  at  large  water  depths  because  the  loss  of  energy  may 
become  significant.  The  long  trailing  wire  required  to  take 
out  the  signal  may  present  problems.  Therefore,  instruments 
tion  to  preclude  the  use  of  trailing  wire  may  be  necessary. 

The  penetrator  could  also  be  used  as  a  device  for  the 
emplacement  of  soil  and  rock  anchors  or  instrument  stations 
in  the  ocean  bottom.  A  penetration  nomogram  based  on  the 
pant  experience  with  such  devices  is  shown  in  Figure  4.5 
and  is  discussed  in  Section  4.2.3  (Young,  1969). 
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7 .0  Immediate  Use  Recommendations  for  Predicting  Object 
Penetration  into  the  Ocean  Bottom 


Recommendations  for  calculating  and  predicting  object 
penetration  into  the  ocean  bottom  depend  very  much  on  the 
conditions  and  circumstances  of  each  particular  problem. 
There  are  a  number  of  important  parameters  that  control  the 
penetration  process.  It  is  of  great  significance  how 
detailed  and  reliably  these  parameters  are  known  or  can  be 
determined . 

The  controlling  parameters  according  to  our  preceding 
considerations  clearly  are: 

1.  The  impact  velocity; 

2.  The  mass,  geometry  and  structure  of  the  penetrat¬ 
ing  object; 

3.  The  impact  configuration  (i.e.  impact  trajectory 
and  geometry); 

4.  The  properties  of  the  ocean  bottom. 

The  impact  velocity  can  be  calculated  fairly  well  by 
the  methods  discussed  in  Chapter  3. 

Which  one  of  the  various  methods  should  be  used  to 
calculate  penetration  depends  now  largely  on  the  reliabil¬ 
ity  with  which  the  four  impact  parameters  are  known.  Clear¬ 
ly,  it  makes  not  much  sense  to  spend  an  excessive  amount  of 
effort  to,  say,  determine  the  impact  velocity  within  a 
fraction  of  one  per  cent  if  the  soil  properties  can  only  be 
ascertained  betv/een,  say,  ±50  per  cent  and  vice  versa. 

Similarly,  it  does  not  make  sense  to  go  through  complex 
calculations  that  try  to  compute  penetration  or  velocities 
to  a  fraction  of  a  per  cent  such  as  equations  (4.13)  or 
(4.21)  if  the  properties  or  parameters  that  have  to  be  used 
in  these  equations  are  not  known  or  have  to  be  guessed. 

Then  again,  some  approaches  or  methods  might  be  more 
suitable  than  others  for  the  particular  set  of  parameters 
known.  At  any  rate,  to  be  most  successful  with  any  calcula¬ 
tion,  one  should  strive  to  have  a  balance  in  the  reliability 
or  the  probability  with  which  all  the  data  entering  into  the 
computation  are  known. 

In  the  absence  of  reliable  information  one  can  always 
work  from  two  extreme  assumptions  and,  narrow  the  solution 
down  to  the  more  or  most  probable  middle  range. 


Except  for  equations  (4.13)  and  (4.21),  none  of  the 
equations  of  Chapters  4  and  5  are  complicated'  or  cumbersome 
to  solve  and  one  would  probably  do  well  to  make  trial  calcu¬ 
lations  with  several  formulae. 


7 . 1  Recommendations  for  Free  Fall  Penetration 

For  an  essential  rigid  body  penetration  of  a  more  or 
less  bulky  (spherical  or  eliptical)  object,  one  should  try 
first  equations  (4.17)  and  (4.23)  if  the  impact  velocity  is 
very  low  such  that  one  would  have  to  assume  the  phenomenon 
to  be  that  of  increasing  area  penetration. 

If  the  impact  parameters  are  such  that  a  constant  area 
penetration  is  indicated,  one  should  try  equations  (4.14), 
(4.27),  and  (4.37).  Also,  equations  (4.43)  and  (4.44) 
would  be  possible  if  a  reasonable  way  exists  to  assess  the 
value  of  the  constants  K. 

As  far  as  equation  (4.43)  is  concerned,  it  is  entirely 
empirical  and  the  coefficients  will  have  to  be  determined 
from  several  series  of  tests .  Any  extrapolation  far  beyond 
the  range  of  the  tests  will  be  hazardous.  Equation  (4.44) 
is  semi-empirical  due  to  the  fact  that  its  structure  basic¬ 
ally  follows  the  Poncelet  solution.  But  again,  the  actual 
coefficients  must  be  determined  from  experimental  results. 
The  difference  is  that  hez’e  extrapolation  may  perhaps  be 
carried  farther  with  higher  confidence  .  Since  experimenta¬ 
tion  increases  rapidly  in  cost  with  the  scale  of  the  exper¬ 
iment,  the  major  part  of  the  experimental  work  could  be  done 
on  small-scale  models  and  only  a  few  experiments  would  be 
necessary  at  large-scale  to  check  the  validity  over  a  wide 
range  . 

The  parameters  for  the  rheological  models,  for  example 
Gj^  and  m  of  equation  (4.23),  can  be  obtained  by  a  curve 
fitting  technique  comparing  curves  obtained  experimentally 
with  those  calculated  theoretically  by  varying  combinations 
of  the  parameters  and  (Tsai  and  Schmid,  1969). 

If  an  impact  penetrometer,  for  example,  furnishes  a 
certain  signal  trace  giving  the  deceleration  as  a  function 
of  time,  a  corresponding  x  versus  time  curve  can  be  calcu¬ 
lated  assuming  estimated  values  for  Gj  and  ni.  If  the  two 
curves  do  not  match  within  acceptable  limits,  a  new  set  of 
parameters  Gj^  and  m  is  assumed  and  the  theoretical  curve 
is  recomputed.  The  comparison  is  repeated  until  all  dis¬ 
crepancies  are  within  acceptable  limits.  This  process  can 
be  computerized  facilitating  the  evaluation  of  a  large 
amount  of  data  and  permitting  a  close  match. 
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7 .2  Recommendations  for  Dynamic  (Driven)  Penetration 

For  dynamic,  driven  penetration  into  sand  or  other 
:cobas.lonless  bottom  sediment,  the  penetration  resistance 
and  bearing  capacity  is  best  calculated  by  using  equation 
(5.8)  with  the  bearing  capacity  factor  given  by  equa¬ 
tion  (5.10).  For  penetration  into  clay  or  cohesive 
materials,  equations  (5.12)  or  (5.14)  are  indicated. 

These  formulae  are  also  indicated  when  the  bearing 
capacity  of  a  driven  pile  has  to  be  calculated.  The  only 
difference  is  that  the  skin  friction  is  more  pronounced 
than  during  driving. 


7.3  Recommendations  for  Powered  Travel  Penetration 


A  missile  or  vehicle  may  penetrate  into  the  ocean 
bottom  while  still  under  power.  For  such  a  case  the  fol¬ 
lowing  procedure  for  calculating  penetration  is  recommended; 

The  power  furnished  by  the  power  plant  during  the 
impact  period  is  considered  to  be  constant  and  equals  q 
( ft . lbs . /sec . )  .  The  duration  of  the  impact  period  T  can  be 
calculated  from  equation  (4.30).  The  energy  supplied  by  the 
power  plant  during  impact  then  is  q*T.  This  energy  can  be 
converted  i.nto  an  equivalent  velocity  uj  to  be  added  to  the 
impact  velocity  Uq  as  a  velocity  correction.  Thereafter, 
the  problem  can  be  treated  as  a  problem  of  free  fall  pene¬ 
tration  (see  Section  7.1)  with  impact  velocity  U2  and: 


where  m  is  the  mass  of  the  missile . 


7 .4  Recommendations  for  Use  of  Sandia  Nomograph 

If  information  on  the  soil  parameter  S  for  the  ocean 
bottom  sediment  is  available  from  prior  tests  or  experience, 
the  penetration  of  a  Sandia  penetrator  may  be  calculated  by 
use  of  the  nomograph  of  Figure  4.5  as  described  in  Section 
4.2.3.  For  this  purpose,  however,  the  impact  velocity  Uq 
must  be  calculated  from  the  impact  velocity  Uj^  at  water  sur¬ 
face  impact  minus  the  change  in  velocity  Au  due  to  the  drag 
loss  between  the  water  surface  and  the  ocean  bottom.  The 
phenomenon  is  the  same  as  that  described  by  equations  (3.22) 
to  (3.25),  (assuming  a  constant  drag  coefficient),  and  the 
impact  velocity  Uq  at  the  ocean  bottom  would  be  given  by: 


-Kh 

U2  =  ui  e  +  Uq 


where  K  is  given  by: 


V  1  o  A 

^  =  2  ^D  tn  Pf 


h  is  the  water  depth,  and  Uq  is  the  free  fall-terminal 
velocity  in  sea  water  for  the  penetrator.  Hence,  the 
bottom  impact  velocity  U2  decreases  exponentially  with 
water  depth  and  approaches  asymptotically  the  terminal 
velocity  of  the  penetrator'  Uq . 

As  a  consequence,  in  deep  water,  little  can  be  gained 
by  dropping  the  penetrator  from  the  air  that  could  not  be 
gained  by  dropping  it  from  a  ship.  By  evaluating  the  pen¬ 
etrator  record  and,  hence,  using  the  penetrator  as  a  pene¬ 
trometer,  a  measure  of  the  strength  of  the  soils  traversed 
can  be  obtained.  However,  considerable  additional  research 
work  appears  to  be  required  to  get  good  quantitative  data 
and  results  for  deep  penetrations. 


7 . 5  Sample  Calculations 

To  indicate  in  detail  the  use  of  the  methods  and  .equa¬ 
tions  recommended,  a  few  sample  calculations  may  be  appro¬ 
priate  . 

To  calculate  free  fall  penetration  we  consider  the 
hypothetical  case  of  a  submarine,  its  hull  intact,  sinking 
to  the  ocean  bottom.  We  assume  the  submarine  to  be  of  .the 
class  Guppy  II  with  a  length  of  306  ft.,  a  beam  of  27  ft., 
a  submerged  displacement  of  2,445  tons,  and  a  minimum  and 
maximum  surface  displacement  of  1500  tons  and  1840  tons 
respectively  (Muga,  1968).  The  first  parameter  to  be 
obtained  is  the  impact  velocity. 

We  may  consider  two  cases; 

case  a):  The  hull  sinking  with  longitudinal  axis  hori¬ 
zontal  due  to,  say,  a  general  flooding  of  the  entire  hull. 

case  b):  The  hull  sinking  vertically  due  to  preferential 
flooding  at  the  bow  and  a  remaining  air  pocket  at  the  stern 
of  the  hull,  thus  giving  a  center  of  buoyancy  ab,ove  the 
center  of  gravity  of  the  partially  flooded  hull.  For  this 
case  we  assume  a  degree  of  flooding  of  75%. 
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The  terminal  velocity  is  calculated  from  (3.21)  as 
follows : 


Total  volume  (submerged)  of  submarine:  2,445  m^ 
Total  gross  weight  of  submarine:  1,500  tons 

O 

Air  space  in  ballast  tanks:  945  m 

Volume  of  steel:  Vg  =  1,500/7.9  =  190  m^ 

O 

Available  space  in  hull:  1500  -  190  =  1,310  m 

Consider  20%  of  hull  space  filled  by 
light  weight  equipment. 

Weight  of  water  flooding  hull  (case  a): 

Wf  =  0.8  •  1310  •  1.0  =  1,.050  tons 

Weight  of  water  flooding  hull  (case  b): 

Wf  =  0.8  •  0.75  •  1310  •  1.0  =  785  tons 

Average  density  of  sinking  object: 


,  1500  +  945  +  1050  .  .  „  ,3 

case  a):  Pq  =  - 27445 -  "  1.425  t/m 


case  b) :  Pq  " 


1500  t  945  +  785 
2,445 


=  1 .325  t/m' 


Assuming  drag  coefficients  for  both  cases  from  Table 
3.7  (since  L/D  -  10),  and  respective  drag  areas: 


case  a):  Ct,  =  0.7  ;  A  =  -r  •  27  •  306  =  5500  ft  =  512  m' 


case  b):  Cn  =  0.5  ;  A 


Ti  =  575  ft^  =  53.4  m 


The  terminal  velocities  Uq  then  are  given  by: 


2  .445  •  0.395 


case  a):  u^  =  1.41  q  .7  '•"Tl^""  l“3"o'  =  2.28  m/sec 


2  445  •  0  295  ^ 

case  b):  u^  =  1.41  oTs'*'-  53.4  ITolo  =  7. 20  m/sec 


The  formula  that  is  applicable  to  this  problem  is 
equation  (4.27)  giving  the  maximum  penetration: 


2  _  1  ""o 


Sax  “  2Apq  ^o  ■  g  2Ap 


Cjj-values  are  extrapolated  beyond  the  available  data  due 
to  the  large  Reynolds  numbers . 
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Assuming  a  bearing  capacity  for  the  ocean  bottom  of  1.5  psi 
=  0.105  kg/cm^,  and  introducing  the  values  found  previously 
(using  metric  units  throughout),  we  obtain: 


case  a): 

X  = 
'^max 


3.495  •  10^  ♦(2.28)^  «  10*^ 

981  •  2  •  512  •  10*+  •  0.105  =  ^^2  cm 


5.95  ft. 


Also,  the  impact  period  is  given  by: 


T  = 


APo 


228 


981 


512 


3.495  ' 
•  10** 


10' 


•  0.105 


=  1.51  seconds 


Since  this  penetration  is  relatively  small  and  to  check 
by  a  second  calculation  whether  this  result  is  in  the  right 
order  of  magnitude,  we  may  try  equation  (4.17)  which  is 
based  on  an  “increasing  area  penetration"  formulation.  How¬ 
ever,  this  equation  was  derived  for  the  penetration  of  a 
sphere  and,  hence,  it  can,  at  best,  give  us  a  rough  approx¬ 
imation  provided  we  convert  the  penetration  area  into  an 
equivalent  circular  area: 

R^tt  =  512  m2,  giving  R  =  12.7  m  =  1270  cm 


Now,  using  (4.17)  we  get: 


^max  ~ 


228 


3,495  »  10^ 
6.24  •  1270  » 


=  470  cm  =  14.4  ft. 


and  the  impact  period: 


3 . 2  seconds  . 


This  second  calculation  shows  that  both  equations  are  con¬ 
sistent  and  essentially  in  agreement.  The  second  results 
indicate  that  if  a  mass  equal  to  that  of  the  submarine  had 
a  spherical  shape  and  would  arrive  with  the  same  impact 
velocity,  its  penetration  would  be  about  240  per  cent  of 
the  one  calculated  by  equation  (4.27). 


case  b): 

Using  again  equation  (4.27)  with  the  impact  velocity 
and  penetration  area  values  for  this  case  yields: 


3  .495  »  10^  »  (7.20)^,  «  10*^ 


981  •  2  •  53,4  •  10 


0  .105 


=  16,400  cm  =  540  ft 


and  the  impact  period  for  this  case  yields: 


720  •  3.495  »  10 
981  •  53.4  •  10‘^  •  0.105 


=  50.2  seconds 


While  these  values  are  probably  somewhat  excessive  due  to  the 
fact  that  no  consideration  was  given  to  the  increase  in  shear 
strength  and,  hence,  bearing  capacity  Po  with  depth,  we  feel, 
nevertheless,  that  they  are  entirely  realistic  values  and 
must  be  within  the  right  range  . 

The  actual  attitude  of  a  submarine  during  free  fall  to 
the  ocean  bottom  most  probably  will  neither  conform  to  case 
a)  nor  to  case  b) .  It  will  have  some  intermediate  attitude 
which,  given  the  usual  ignorance  about  the  details  of  a  mis¬ 
hap,  is  difficult  to  assess.  Nevertheless,  the  sample  calcu¬ 
lations  show  the  limits  that  may  be  expected.  Together  with 
some  calculations  at  several  different  yaw  angles  one  could 
probably  construct  a  diagram  for  probable  penetration  vs  . 
yaw  angle  at  impact.  The  final  discovery  photographs  of  the 
submarine  Thresher  (Buchanan,  1965)  showed  that  only  the 
tail  section  and  the  sail  projected  above  the  ocean  bottom 
sediment  a  distance  of  about  twelve  and  eight  feet  respec¬ 
tively.  Since  the  sediment  shear  strength  at  the  Thresher 
disaster  site  was  about  0.3  psi,  thus  giving  an  ultimate 
bearing  capacity  Po  =  1»5  psi,  these  observations  confirm 
the  general  correctness  of  the  results  obtained. 

In  general,  we  wish  to  emphasize  that  possibly  none  of 
the  various  equationsi  or  methods  will  turn  out  to  yield  the 
one,  only  and  always  correct  result. 

There  will  most  likely  be  varying  and  changing  success 
with  several  of  the  proposed  approaches,  and  in  the  absence 
of  any  extensive  proof  testing,  it  is  difficult  to  even 
hazard  a  guess  which  one  is  likely  to  be  most  successful. 

Also,  successful  prediction  depends  very  much  on  the 
success  with  which  the  properties  or  the  parameters  that 
enter  into  any  formula  can  be  measured  or  otherwise  ascer¬ 
tained.  Hence,  the  quality  of  any  penetration  prediction 
really  cannot  be  assessed  until  some  practical  experience 
has  been  gained. 

There  is,  nevertheless,  a  tremendous  possibility  and 
opportunity  for  improving  our  potential  prediction  capabil¬ 
ities  and  their  reliability  by  perhaps  at  least  one  order 
of  magnitude  or  more , 


The  details  of  this  opportunity  are  outlined  in  Chap¬ 
ter  8 . 
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8 .0  Recommendations  for  Future  Research  and  Development 

To  advance  our  present  capabilities  of  predicting,  as 
well  as  achieving  penetration  into  the  ocean  bottom,  a  num¬ 
ber  of  research  and  development  projects  can  and  should  be 
initiated.  To  increase  the  reliability  of  penetration  pre¬ 
dictions,  those  parameters  that  are  the  "weak  links”  in  the 
logic  of  the  prediction  process  should  be  attacked  first . 
The  two  weakest  links  are  clearly  the  reliability  of  the 
strength  assessment  of  the  ocean  bottom  soils  and  the  reli¬ 
ability  of  the  various  formulae  discussed  or  derived  in 
Chapters  4  and  5 . 

There  exists  very  little  evidence  on  the  ”in  situ” 
strength  of  ocean  bottom  soils  and  the  evidence  on  the 
relationship  between  the  ”in  situ”  shear  strength  and  the 
strength  obtained  from  measurements  on  core  samples  con¬ 
sidering  the  data  of  Vey  and  Nelson  (1967) 'is  inconclusive. 

It  is  not  that  we  can  claim  an  abundance  of  data  from 
cores.  Most  of  the  strength  information  available  is  that 
obtained  from  core  samples  reported  by  Richards  (1961), 
Moore  (1962),  and  Keller  (1967).  Considering  the  vastness 
of  the  ocean  floor,  the  volume  of  information  is  nothing  to 
be  complacent  about.  Therefore,  there  exists  an  urgent 
need  to  develop  apparatus  and  instrumentation  to  measure 
the  shear  strength  of  ocean  bottom  soils  in  situ  in  a 
manner  that  is  fast,  reliable  and  cheap  such  that  many 
measurements  can  be  performed  over  a  small  area  in  one 
trip.  Only  then  will  it  become  practically  a  routine  for 
any  oceanographic  activity  to  carry  out  also  a  bottom 
strength  determination.  The  work  reported  by  Kretschmer 
(1967)  is  considered  truly  a  pioneering  effort  in  this 
direction . 

Preferably,  however,  rather  than  work  with  a  static- 
l.oad  type  plate  bearing  test,  one  might  investigate  a 
dynamic  impact  test  such  as  that  used  by  Scott  (1967)  or 
Schmid  (1966).  It  is  believed  such  instrumentation  will 
lead  to  cheaper,  more  manageable  equipment.  Any  such 
development  of  measurement  instrumentation  must  account 
for  the  fact  that  the  water  content  of  the  ocean  bottom 
soils  usually  is  far  above  the  liquid  limit  and,  therefore, 
the  conceptual  formulation  of  the  stress-strain  behavior 
should  be  one  that  is  based  on  the  principles  of  visco¬ 
elasticity  rather  than  on  those  of  elasto-plasticity .  The 
treatment  of  most  problems  in  conventional  soil  mechanics 
is  based  on  the  elasto-plastic  failure  concept.  This  is, 
perhaps,  an  admissible  simplification  for  most  conventional 
foundation  problems.  It  becomes  an  entirely  inadequate 
oversimplification,  however,  for  ocean  bottom  soils.  The 
data  reported  by  Kretschmer  (1967)  on  average,  constant 
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penetration  rates  and  stress  relaxation  for  in  situ,  sea 
floor  plate  bearing  tests  clearly  support  this  observation. 
Thus.,  it  appears  that  an  approach  such  as  suggested  by 
Schmid  and  .Kitago  (1965)  or  Tsai  and  Schmid  (1969)  is 
necessary  for  the  evaluation  of  in  situ  shear  strength 
data . 


There  are  two  recommendations  for  future  research  and 
development  that  can  be  deduced  from  these  observations, 
namely: 

1)  Research  on  the  visco-elastic  stress-strain  char¬ 
acteristics  of  typical  ocean  bottom  soils  should  be  carried 
out  with  soil  constituents,  porosities,  void  ratios,  water 
contents  and  pore  pressure  environments  as  prevail  on  the 
ocean  bottom. 

2)  Apparatus  and  instrumentation  should  be  developed 
that  will  permit  rapid,  reliable,  and  cheap  routine  deter¬ 
mination  of  the  ocean  bottom  sediment  stress-strain  and 
strength  characteristics .  The  author  believes  that  the 
most  promising  approach  in  this  direction  would  be  the 
development  of  an  impact  penetrometer  for  deep  underwater 
application  using  accelerometers,  such  as  described  by 
Scott.,  Schmid,  or  the  Sandia  device.  However,  we  do  not 
wish  to  exclude  other  methods  that  may  have  promise  such 
as,  e.g.,  an  in  situ  vane  shear  test. 

The  second  weak  link  in  the  reliability  of  the  pene¬ 
tration  px-ediction  capability  is  the  validity  of  the  vari¬ 
ous  formulae  reviewed  in  Chapter  4.  Again,  two  recommenda¬ 
tions  are  suggested  to  improve  existing  capabilities, 
namely: 

3)  A  computer  program  development  should  be  initiated 
patterned  after  the  HEMP  Code  and  the  coupled  Euler-Lagrange 
(CEL)  formulation  as  explained  in  Section  4.2.5.  Such  a 
program  should  be  developed  with  the  view  in  mind  that  any 
plausible,  rheological  stress-strain  and  failure  character¬ 
istic  that  might  result  from  work  as  recommended  in  1)  can 
be  accommodated. 

4)  A  carefully  controlled,  model  test  program  (on  a 
laboratory  model  scale)  of  idealized  objects  (preferably 
spheres)  free-falling  through  water  onto  ocean  bottom  model 
test  beds  should  be  carried  out  to  test  the  validity  of  the 
various  formulae  both  for  a  larger  number  of  small-scale 
and  some  larger-scale  model  tests.  The  smaller-scale  model 
tests  will  permit  the  testing  under  a  variety  of  conditions; 
the  larger-scale  model  tests  should  provide  sufficient  data 
for  extrapolation  to  check  the  scaling  laws  and  similitude 
considerations.  Perhaps  later,  also,  some  data  for  proto¬ 
type  confirmation  may  becoiil?  available. 
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At  any  rate,  the  complete  absence  of  data  on  the  reli¬ 
ability  of  the  available  formulae  for  assessing  ocean  bottom 
penetration  should  be  corrected.  It  is  believed  best  to  be 
accomplished  by  model  scale  experiments.  We  wish  to  note, 
however,  that  we  believe  maximum  benefit  will  result  from 
such  an  experimental  program  only  after  recommendations  1) 
and  3)  have  been  carried  out  and  the  experimenters  have  the 
benefit  of  the  results  obtained  there. 

Alternately,  if  time  is  pressing,  the  experimental  work 
might  be  carried  out  concurrently  with  the  work  recommended 
under  1)  provided  it  is  done  by  the  same  group  or  at  least' 
carried  out  in  closest  cooperation  between  both  research 
task  forces. 

After  a  clearer  picture  is  obtained  on  the  stress- 
strain  and  strength  characteristics  of  ocean  sediments  from 
programs  initiated  to  follow  up  recommendations  1)  through 
4),  it  might  also  be  useful  to  reexamine  the  wealth  of 
acoustical  data  that  is  available  for  ocean  bottom  soils 
and  attempt,  perhaps,  a  better  correlation  between  such 
data  and  stress-strain  or  strength  parameters  of  the  sedi¬ 
ments.  Prior  work  in  this  direction  reported  by  Nafe  and 
Drake  (1963)  and  Horn  et  al .  (1968),  while  not  completely- 
encouraging,  nevertheless  gives  an  indication  of  some 
promise  that  possible  relationships  might  be  explored. 

Horn  et  al .  show  that  some  causal  relationship  exists 
between  wave  velocity  on  the  one  hand  and  porosity,  wet 
density,  and  moisture  content  on  the  other.  No  such  clear- 
cut  indication  is  given  for  the  shear  strength  data.  Faas 
(1968)  reports  an  inverse,  linear  relationship  between  the 
reflection  of  an  acoustic  impulse  from  the  ocean  bottom  and 
the  porosity  of  the  bottom  sediment.  This  relationship  has 
been  shown  to  be  consistent  for  sediments  taken  from  differ¬ 
ent  oceans  and  different  environments.  For  these  reasons 
we  believe  it  worthwhile  to  take  another  look  at  the  data 
to  see  whether  some  useful  correlation  can  be  established 
between'  the  wealth  of  acoustic-seismic  data  and  strength 
parameters.  For  example,  a  relationship  should  exist  be¬ 
tween  observed  dispersion  of  acoustic  waves  and  the  viscos¬ 
ity  of  the  sediment.  Thus,  we  recommend  that 

5)  The  relationships  between  acoustic  properties  and 
data  of  ocean  bottom  sediments  and  their  stress -rstrain- 
strength  parameters  be  investigated  particularly  in  the 
light  of  experiences  gained  under  work  done  according  to 
recommendation  1). 

Finally,  with  respect  to  efforts  for  increasing  our 
capabilities  of  achieving  deep  ocean  bottom  penetration, 
we  feel  that  deep  penetrations  in  deep  waters  eventually 
will  become  an  operational  necessity  and  it  will  be  useful 
to  have  an  understanding,  at  least,  of  the  problems  with 


which  one  will  be  confronted .  Any  deep  penetration  by 
piles,  for  example,  in  the  deep  ocean  will  have  to  be  done 
by  underwater  pile  driving  since  the  present  state  of  the 
art  is  pushing  close  to  the  limits  of  technological  feasi¬ 
bility  as  far  as  driving  piles  from  above  the  water  is  con¬ 
cerned.  Piles  having  a  length  of  400  to  600  feet  being 
driven  above  the  water  in  water  depths  of  up  to  350  feet  is 
the  present  advanced  state  of  the  art.  Any  significant 
increase  in  water  depth  and,  consequently,  in  pile  length 
will  result  in  piles  either  so  slender  that  they  will  buckle 
under  the  driving  loads  or  so  big  and  bulky  as  to  make  them 
unwieldy  and  unmanageable  even  with  the  heaviest,  water¬ 
borne  equipment.  Hence,  the  only  reasonable  solution  even¬ 
tually  will  be  to  drive  the  piles  underwater. 

There  are  a  number  of  possible  approaches  to  the  prob¬ 
lem  of  driving  piles  below  the  water.  Whichever  method  is 
best  clearly  will  depend  upon  the  water  depth  at  which  an 
underwater  pile  driver  has  to  operate.  The  various  possi¬ 
bilities  may  be  grouped  into  three  classes: 

A)  Atmospheric  exhaust  tools. 

B)  Closed  systems  (hydraulic  or  electric)  with 
power  supplied  through  an  umbilical  cord  from 
a  surface  barge  or  tender. 

C)  Self-contained  systems  with  power  package  and 
fuel  supply  on  board  for  deep  submergence 
operations  . 

A)  Atmospheric  exhaust  tools  are  steam  or  compressed 
air  hammers  with  exhaust  to  the  atmosphere  through  exhaust 
hoses.  This  is  the  current  state  of  the  technology.  The 
prototype  tool  is  essentially  a  conventional,  MKT  single- 
acting  pile  hammer  modified  with  pressure  seals  to  operate 
under  water.  As  mentioned  before,  the  maximum  water  depth 
at  which  it  has  been  used  in  the  past  is  162  feet.  This, 
perhaps,  can  be  pushed  to  300+  feet.  But  the  method  becomes 
rapidly  inefficient  at  water  depths  in  excess  of  200-250 
feet  due  to  the  pressure  losses  and  the  large  diameter  hoses 
that  will  be  required.  The  long  hoses  present  quite  a  hand¬ 
ling  problem  during  operations.  Nevertheless,  achieving 
operational  capability  within  this  water  depth  range  (say, 
up  to  300  feet)  probably  can  be  accomplished  at  a  relatively 
low  cost . 

B)  Closed  systems  are  hydraulically-powered  pile 
drivers  with  a  prime  mover  (hydraulic  pump)  on  a  barge  driv¬ 
ing  hydraulic  motors  or  pistons  on  the  driving  tool.  The 
prime  mover  and  driving  tool  are  connected  by  an  "umbilical 
cord"  in  the  form  of  hydraulic  hoses.  Since  such  a  system 
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could  be  operated  at  pressures  of  1,500,  3,000  or  5,000  psi 
respectively,  the  corresponding  sizes  of  the  required  pumps, 
motors  and/or  pistons  become  progressively  smaller.  Also, 
since  internal  pressure  in  the  hose  as  well  as  external, 
ambient  pressure  increases  linearly  and  at  almost  or,  if 
sea  water  is  used  as  hydraulic  fluid,  at  exactly  the  same 
rate  with  water  depth,  no  serious  depth  limitation  would 
exist  because  of  the  pressure  capacity  limits  of  the  hose. 

The  limitations,  again,  would  come  from  the  friction  losses 
in  the  hoses  and  the  handling  problems  due  to  the  hose 
length.  We  believe  that  these  limitations  may  present  ser¬ 
ious  handicaps  if  such  a  system  is  to  operate  at  water 
depths  beyond  1,500  to  2,000  feet. 

An  alternate  method  would  be  to  supply  power  from  the 
surface  by  an  electric  generator  set  and  electric  cables. 

This  method,  however,  will  present  even  more  serious  prob¬ 
lems  of  sealing  the  electrical  motors  and  cable  connections. 

It  has  the  additional  disadvantage  of  being  quite  hazardous 
in  case  the  cable  gets  sheared  or  damaged.  Also,  the  elec¬ 
tric  motors  are  much  more  bulky  than  hydraulic  motors  for  the 
same  power  output. 

As  of  this  instance,  there  is  one  driving  tool  on  the 
commercial  market  that  is  suitable  for  adaptation  to  under¬ 
water  environment  and  the  pressure  ranges  that  appear  prac¬ 
tical  for  this  type  of  driving  system.  This  is  the  MKT  V-10 
hydraulic  pile  vibrator.  Operating  specifications  for  this 
machine  are  given  in  Figure  6.29. 

Sufficient  experience  exists  in  vibratory  pile  driving 
that  one  can  safely  say  it  has  been  proven  out  technologically. 
As  discussed  in  Section  5.9,  however,  pile  vibrators  in  the 
power  range  commercially  available  at  this  time  (1969)  still 
have  definite  limitations  with  respect  to  the  maximum  driving 
force  available  under  tough  driving  conditions.  Fortunately, 
there  are  vibrators  now  under  development  with  a  significantly 
stepped  up  maximum  dynamic  force  that  should  put  them  into 
the  same  range  as  conventional,  heavy,  single-acting  hammers. 

Since  hydraulically-driven  vibrators  appear  to  be  ideally 
suited  for  underwater  application  and,  moreover,  since  such 
vibrators  are  already  developed  commercially  and,  hence,  may 
be  adapted  to  underwater  use  at  relatively  low  costs  of 
development,  we  recommend  that 

6)  the  adaptation  of  these  tools  for  underwater  opera¬ 
tion  be  investigated. 

The  limits  of  suitability  should  be  explored,  as  well  as 
the  design  changes  that  would  become  necessary  to  achieve 
underwater  operational  capability,  and  the  limits  up  to  which 
an  underwater  hydraulic  vibrator  powered  by  hydraulic  pumps 
above  water  is  practicable. 


Another  similar  tool  in  this  group,  equally  suitable 
for  underwater  pile  driving  applications,,  is  a  hydroacous¬ 
tic  oscillator  developed  by  the  Electronics  Division  of 
General  Dynamics  Corporation  as  an  outgrowth  of  underwater 
Sonar  development  (Anon.,  1966).  The  principle  of  this 
tool  is,  essentially,  a  hydraulic  piston  that  serves  simul¬ 
taneously  as  a  valve  and  thus  represents  a  machine  with 
only  one  moving  part  (see  Figure  8.1).  This  feature  makes 
it  extremely  attractive . 

The  tool  could  be  used  as  a  linear  oscillator  or,  with 
appropriate  anvil  coupling,  as  an  impact  tool  (Figure  8.2). 
The  configuration  and  dimensions  of  such  a  machine  depend 
largely  on  the  frequency  of  the  hydraulic  piston-valve  and 
the  required  power  output.  Figure  8.3  shows  some  typical 
dimensions  for  such  a  tool  operating  at  a  1500  psi  pressure 
differential  and  various  frequencies  and  power  outputs. 

The  disadvantages  of  the  tool  lie  in  its  extremely 
narrow  frequency  range  --  the  tool  is  essentially  a  single¬ 
frequency  machine.  Also,  the  machine  is  not  yet  developed 
commercially  and,  hence,  there  may  be  considerable  develop¬ 
ment  costs.  There  is  no  question,  however,  that  the  system 
has  promise  for  underwater  applications. 

C)  The  self-contained  system  is,  essentially,  a  closed 
system  having  prime  mover,  fuel  supply,  and  driving  tool 
aboard  and,  hence,  such  a  system  requires  nothing  more  than 
a  cable  connection  for  lifting  and  lowering.  Thus,  it  has 
the  largest  depth  capability  of  all  possible  systems.  There 
are  no  problems  of  friction  losses  nor  those  of  handling 
long  lengths  of  hoses.  The  power  for  the  required  intensi¬ 
ties  and  duration  of  operations  is  most  practically  obtained 
by  providing  chemical  steam  through  burning  hydrogen  and 
oxygen  or,  preferably,  hydrogen  and  hydrogen  peroxide.  The 
actual  driving  too,  then,  can  be  an  especially  sealed,  con¬ 
ventional  type  steam  hammer.  A  fairly  detailed  concept  for 
such  a  system  has  been  proposed  by  the  MKT  Corporation 
(Anon.,  1967).  Again,  there  will  be  considerable  costs 
associated  with  such  a  development  but,  eventually,  it  will 
have  to  be  undertaken. 

For  these  reasons  we  recommend  that 

.7)'  Feasibility  studies  on  tne  design  and  development 
of  underwater  pile  driving  systems  possible  within  the 
present  and  future  state  of  the  art  should  be  carried  out. 
Such  studies  might  reasonably  begin  with  assessing  the 
underwater  capabilities  achievable  with  slight  modifica¬ 
tions  of  presently  commercially  available  vibratory  equip¬ 
ment  and  proceed  from  there  with  establishing  the  design 
envelopes  for  underwater  pile  drivers  of  increasing 


capabilities  both  with  respect  to  water  depth  and  opera¬ 
tional  capacity.  Concurrent  with  such  studies  or,  perhaps, 
even  preceded  by  it,  there  should  be  a  study  made  on  pro¬ 
jected  U.S.  Navy  requirements  for  operational  capabilities 
on  the  ocean  floor.  As  difficult  as  this  may  appear  at 
this  time,  we  believe  much  can  be  learned  from  a  record  of 
the  growth  of  past  requirements.  This  would  also  permit  a 
reasonable  phasing  of  any  development  work  that  will  be 
required  as  a  consequence  of  the  study. 

We  conceive  these  studies  as  a  two-phase  or  two¬ 
pronged  program.  One  phase  would  investigate  the  techni¬ 
cal  possibilities  within  a  short  or  intermediate  time  range 
by  maximum  use  of  existing  equipment  and  technology  to 
advance  present  capabilities  at  relatively  low  cost.  We 
consider  this  a  high  priority,  high  payoff  item  capable  of 
significantly  advancing  the  state  of  the  art  at  a  relatively 
low  cost.  The  second  phase  should  be  a  long  range  look  at 
technological  capabilities  and  possibilities,  particularly 
assessing  overall  development  and  operation  costs  for  the 
various  possible  alternates  such  as,  for  example,  the  self- 
contained  systems  with  great  depth  capabilities. 

Since  many  future  deep  sea  installations  are  very 
likely  to  be  founded  on  piles  penetrating  deep  into  the 
ocean  bottom  sediment,  which  may  be  driven  by  vibrations, 
we  recommend  that 

8)  The  bearing  capacity  of  vibro-driven  piles  be  estab¬ 
lished  experimentally  both  in  an  ocean  environment  and  in 
a  near  land  or  land  foundation  situation.  There  is  only 
very  limited  data  available  both  from  prototype  and  model 
tests  on  the  bearing  capacity  of  vibro-driven  piles.  Also, 
some  driving-bearing  capacity  relationship  would  be  highly 
desirable  not  only  as  a  construction  quality  control  device 
but  also,  ideally,  as  a  direct  relationship  similar  to  the 
various  pile  driving  formula  used  for  conventional  pile 
foundations .  Although  the  validity  of  many  pile  driving 
formulae  is  subject  to  grave  doubts  and  reservations,  the 
possibility  of  some  relationship  might  be  explored. 

Here,  particularly,  the  new  developments  reported  by 
Goble  et  al.  (1967)  and  discussed  in  Section  5.8  look  most 
promising . 

Another  possibility  is  the  use  of  the  one-dimensional 
wave  equation  (Lowery  et  al.,  1967). 

Also,  birring  capacity  formulae  have  been  derived  for 
terrestial  foundations  and  applications.  Their  applicabil¬ 
ity  in  the  marine  environment  at  the  high  porosities  and 
high  pore  pressures  prevailing  in  the  marine  environment 
should  be  reviewed. 
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Explosive  anchors  are  an  effective  and  relatively 
cheap  method  of  achieving  penetration  into  the  upper  stra¬ 
tum  of  the  ocean  bottom.  Hence,  we  recommend 

9)  Continued  development  and  proof-testing  of  explo¬ 
sively  driven  anchors,  umbrella  piles  and  similar  devices. 

The  Sandia  penetrator  has  applicability  as  a  penetra¬ 
tion  device  in  shallow  water  depths  only.  As  we  have  shown 
its  bottom  impact  velocity  decreases  exponentially  with 
water  depth  and  asymptotically  approaches  that  of  its  ovm 
terminal  velocity  in  water.  We,  therefore,  see  its  useful¬ 
ness  as  a  penetration  device  severely  limited  to  shallow 
water  depths.  However,  we  believe  its  potential  as  a  pene¬ 
trometer  should  be  explored.  Its  use  could  be  conceived 
as  a  free  falling  penetrometer  (although  we  believe  a  spher 
ical  shape  to  be  preferable  for  such  a  purpose)  or,  more 
promising,  as  a  deep  penetration  device  that  is  discharged 
from  a  frame  lowered  to  the  bottom  (by  an  explosive  or  by 
a  hydraulically-powered  charge)  into  the  deeper  layers  of 
the  ocean  bottom  sediment.  Therefore,  we  recommend  that 

10)  The  use  of  a  Sandia-type  penetrometer  for  measur¬ 
ing  the  penetration  resistance  at  large  depth  in  the  ocean 
bottom  sediment  be  explored. 

In  summary,  then,  we  recommend  the  following  future 
work  to  advance  the  state  of  the  art  of  predicting  as  well 
as  achieving  ocean  bottom  penetration: 

1)  Investigation  of  the  viscoelastic  stress-strain 
and  strength  characteristics  of  ocean  bottom 
sediment  soils  . 

2)  Development  of  reliable  apparatus  for  quick  and 
cheap  routine  determination  of  the  strength  of 
ocean  bottom  soils. 

3)  Development  of  a  computer  program  similar  to 
the  HEMP-CEL  Code  for  calculating  penetration 
that  is  capable  of  accommodating  any  visco¬ 
elastic  stress-strain-strength  law. 

4)  Experimental  model  studies  of  actual  free-fall 
penetration  phenomena  and  their  comparison  with 
available  formulae. 

5)  Review  of  the  possible  relationships  of  acousti¬ 
cal  data  to  -Strength  and  viscosity  properties  of 
ocean  bottom  soils. 

6)  Feasibility  studies  for  the  development  of  under¬ 
water  pile  drivers  at  various  depth  and 
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operational  capability  levels,  both  with  a  short 
range  and  longer  range  objective.  These  studies 
should  be  coupled  with,  or  preceded  by,  a  study 
of  U.S.  Navy  operational  requirements  projections. 

7)  Establishment  of  the  bearing  capacity  of  vibro- 
driven  piles  both  in  an  ocean  environment  and  on 
land,  and  development  of  some  relationship  .between 
driving  parameter(s)  and  bearing  capacity  both  for 
vibro-and  impact-driving. 

8)  Continued  development  and  proof-testing  of  explo¬ 
sive  anchoring  devices. 

9)  Study  of  the  application  of  a  Sandia  type  pene- 
trator  lowered  in  a  frame  to  the  ocean  bottom  and 
then  discharged  into  the  deeper  sediment  layers 
as  a  penetrometer  to  measure  the  strength  of  the 
deeper  layers  . 


Princeton,  New  Jersey 
March  3i  ,  1969 
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Figure  2.1 

Shear  Strength  Distribution  of  Bottom  Sediments 
in  the  North  Atlantic  (after  Keller,  1967) 


Figure  2.2 

Shear  Strength  Distribution  of  Bottom  Sediments 
in  the  North  Pacific  (after  Keller,  1967) 
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figure  2 . 3 
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aeita  front  sediments 


Figure  3.3  Drag  coefficients  of  two-dimensional  forms 
of  thickness-chord  ratio  d/c  =  0.25. 


STREAMLINE  PROFILE  LAMINAR  FLOW  PROFILE 


Figure  3.4  Typical  streamline  and  laminar  flow  profiles 
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Figure  3.5  Drag  coefficients  of  various  blunt  forms. 
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Figure  3.8.  Drag  coefficient  Cp  *  aFp/Afu^  of  the 
sphere  in  terms  of  the  Reynolds  num*> 
bers.  Curve  1:  Cj  *  24/R,  theory  of 
very  slow  motion.  (A  *  cross  section 
area  of  sphere.) 
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Figure  3.9 


Drag  Coefficients  C 
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Figure  3.12  Terminal  Velocity  for  Cylinder 

L/D  =  10  falling  with  Axis  Horizontal 


AVERAGE  OBJECT  DENSITY  (fi/ee) 

Figure  3.13  Terminal  Velocity  u^  for  Cylinder 

L/D  =  10  falling  with  Axis  Vertical. 


Test  Series  Ho.  38 

Target  Material:  Clayey,  Sandy  Silt,  Undisturbed 
Calibration  Factors:  Penetrometer:  K  =  8.6  mV/g 

r^,:  K  =  31  mV/g;  r2 :  K  =  34  mV/g 


Drop  Height:  H  =  275  cm 


Oscilloscope  Calibration: 


Penetrometer: 


s  200  mV/cm 
Mt  =  2  msek/cm 


J 


Figure  4.3 


Typical  Impact  Penetrometer  Data 
(after  Schmid,  1966) 


Figure  4 . 4 

Penetration  of  Steel  Spheres  into  Target  Materials 

(after  Kornhauser) 


Curve  Projectiles 


Target 


1 

2 


4 

5 

6 

7 

8 
9 

10 


4.3-in.  cannon  balls 
Many  sizes,  high  ballistic 
density 

0.356-in.  steel  balls 
0.356-in.  steel  balls 
Many 
Many 

Many 

3/16-in.  steal  balls 
9/32-in.  and  7/16-in. 

steel  balls 
Many  high  ballistic  ’ 
density 

Many  high  ballistic 
density 


Sand  and  gravel 
Sand 

Loose  sand 

Weakly  cemented  sand 
Rock  -  250  psi  crushing  strength 
Rock  -  2,000  psi  crushing  ^ 

strength 

Rock  -  10,000  psi  crushing 
strength 

Red  sandstone  -  2,780  psi 
Granite  -  16,500  psi 
Soft  rock 

Rock  or  sand  -  hypervelocity 
cratering 


11 


PENETRATION  NOMOGRAM 


Figure  4.5 


Penetration  Nomogram  for  Sandia  Penetrator. 
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Figure  4.6 


HEMP  Code  Calculation  of  Steel  Cylinder 
Impacting  Upon  Aluminum  Plate 
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Figure  4.8 

Acceleration  Components  for  Oblique  Impact 
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aring  Capacity  Factor  Nq  as  a  Function  Bearing  Capacity  Factor  as  a  Function 

Friction  Angle  4  (after  Meyerhof).  Friction  Angle  (after  Meyerhof). 


ANCLE  OF  INTEimAL  FMCTKM  ^  « 

Figure  5.9 

Comparison  of  Theoretical  Bearing  Capacity 
Factors  Nq  for  Strip  Foundation. 
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Figure  5.12 

Theoretical  and  Observed  Bearing  Capacities 
of  Foundations  at  Surface  (after  Vesic). 
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Figure  5.15 

Analysis  of  Punching  or  Local 
Shear  Failure  (after  Vesic). 


K.  ton  S  onit  N 


Figure  5.16 

Theoretical  Bearing  Capacity  Factors  Nq  at  Greater  Depth 
and  Comparison  of  Model  Test  Results  (after  Vesic,  1967) 


5/Ur  STEEL  ROD  WELDED  TO  SHATT 


weight  of  corer 


weight  of  corer 


(b)  (c) 

Figure  6 . 2 

(a)  The  ‘’Sphincter”  120  mm  Corer 

(b)  and  (c)  Electrical  Release  Mechanism 


•SUfH  in  operation  of  ftee<orer. 


Figure  6.5 

Diagram  of  free-corer  with  polyethylene  bottle  (a) 
which,  when  filled  with  gasoline,  acts  as  pressure- 
proof  buoyant  chamber.  Core  barrel  attached  to 
buoyant  chamber  fits  loosely  inside  of  expendable 
weight-casing  (b). 


Figure  6.6(a) 

Assembly  drawing  of  vibrator  housing  and 
bore  barrel^  A,  electrical  connector;  B, 
switch;  C,  piston  wire;  D,  plastic  liner; 
E,  locking  pin  socket;  F,  piston;  G, 
removable  filler. 


Figure  6.6(b) 

Schematic  diagrams  of  core  sampler  in 
operation.  A,  descending;  B,  coring;  C, 
extracting;  D,  ascending. 


Principle  of  Vibratory  Coring.  1)  hammer, 
2)  anvil,  3)  excentric,  4)  spring. 


Figure  6.7(b) 

Coring  Assembly.  1)  Diesel  engine,  2)  hydrau 
lie  pump,  3)  mud  pump,  4)  hosing,  5)  vibro- 
hammer . _ 6)  bnig-h  mA-hn-w 
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Figure  6.9 

Details  of  BUDOCKS  STATO  mooring  anchor 


Figure  6.11 

A  platforim  being  toued  to  sea. 


Figure  6.14 
Mat-type  jack-up  rig 


Figure  6.15  Schematic  of  mat-type  jack-up. 
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Stmi-iubmenibh  Slaflo  (ttmbU  white  afloat)  mobtla  drilling  platform  will  drill 
for  ShtU  and  Es$o  In  North  Sam  when  delivered  by  m  British  sMpyard 
nest  yaor.  Schematic  view  shows  large  working  area  and  arrangement 
of  cranes  that  will  enabU  highly  eflSdent  drUltng  and  loading  of 
supplies.  Overall  height  of  unit.  Including  derrick.  Is  275  ft. 


Figure  6.17 
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Platform  Construction 
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Figure  6.22 

Example  of  design  and  assembly  schedule 
for  offshore  pile  (after  McClelland  et  al . ) 


REQUIRED  ULTIMATE 
AXIAL  CAPACITIES* 
IJ  -  to  3000T 
P,  “tolOOO^ 

HORIZONTAL  LOADS* 
R  -  to  ISO’^ 


Figure  6.23 

Range  of  pile  loads  used  in  current  offshore 
structure  design  (after  McClelland  et  al . ) . 
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ULTIMATE  PILE  LOAD,  TONS 

Figure  6 .25 

Example  of  the  load  range  represented  by 
one  major  source  of  documented  pile  load 
tests,  as  compared  to  typical  offshore 
pile  load  (after  McClelland  et  al . )  . 


.le  Hammer  with 


Figure  6.27 

Vulcan  040  Offshore  Single-Acting  Steam  Hammer 
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HYDRAULICALLY 

POWERED 


it  Vibration  Complatsly  iaolatad  From  Crana 

*  Safa  Hydraulic  Powar 

*  Narrow  14“  Width 

it  Not  Vunarabla  to  Damage 
it  Quick  Starting  &  Stopping 
it  Quick  Clamping  &  Undamping 
it  No  Overload  on  Starting  &  Driving 
it  Balanced 

it  Intamal  Force  Lubrication  System 

*  Saalad  Against  Submersion  in  Watar 

it  Capable  of  Adding  Non-Vibratory  Weight 
it  Interchangeable  Jaws 


SPECIFICATIONS 


Frequency  . 

Available  Horsepower  .  .  .  . 

Weight . 

Maximum  Lifting  Fpree  .  .  . 
Start  or  Stopping  Time  .  ,  . 
Clamp  Open  or  Closing  Time  . 
Max.  Hydraulic  Preswre  .  .  . 
Dynamic  Force  (At  1850  CPM) 
Clamping  Force  ...... 


1650  <  1850  CPM 
.  .  .  110 
.  .  .  .  5TON 
.  .  .  .  20TON 
5  -7  SEC. 
2  -3  SEC. 
.  .  3000  PSI. 

112,000  LBS. 
150,000  LBS. 


Figure  6.29 


V-10  Hydraulic ■ Pile 
Corporation)  and  its 
Parameters . 
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Performance 
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Division  of  Koohring  Compony 
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The  Earth-Penetration  Event 
(sketches,  no  scale) 
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MOTE:  Scal«  of  record  distorted  for  clarity. 


Figure  6.32 


Example  of  Offshore  Probe  using  Projectiles 


Figure  6.33 

Acceleration  Signatures  for  Nominal 
Penetrometer . 


w 


INFERRED  PROFILE  FROM 
PROJECTILE  DATA 


MATERIAL 

MATERIAL 

MATERIAL 

MATERIAL 

MATERIAL 

MATERIAL 

MATERIAL 


DENSE;  GRANULAR;  EITHER  LARGE  GRAINS,  OR  HIGHLY  ERRATIC. 
MEDIUM  DENSE  TO  DENSE;  GRANULAR,  SAND  SIZES;  SEEMS  UNIFORM. 
SOFT;  COHESIVE. 

SIMILAR  TO  B  . 

SIMILAR  TO  A  . 

ROCK,  MODERATE  TO  WEAK  STRENGTH. 

UNCERTAIN  IF  ROCK  AT  LEVEL  SHOWN. 


EXAMPLE  ONLY,  NOT  CHECKED  FOR  KINEMATIC  CONSISTENCY. 


Figure  6.34 

Example  of  Soil  Investigation  Using  High-Speed 

Projectiles 


HYDROACOUSTIC  OSCILLATOR 


LOW  PRESSURE 


ACOUST 


Unclassified  _ 


Stvutilv  Cldssific^ation 


DOCUMENT  CONTROL  DATA  •  R  &  D 


tSi'i'unty  ot  title,  Itotly  ot  ttbsttuct  aiiu  mUeKinti  4»nnQfittl(m  niunt  be  entefetl  when  the  ovcrMlt  report  In  elataHleU) 


t  ORIGINATING  AcTivtxy  (Cofpotste  author) 


REPORT  SECURITY  CLASSIFICATION 


Werner  E.  Schmid*  P.E. 
Consulting  Engineer 
Princeton,  New  Jersey  08540 


3.  REPORT  TITLE 

PENETRATION  OF  OBJECTS  INTO  THE  OCEAN  BOTTOM  (The  State  of  the  Art) 


4.  DESCRIPTIVE  NOTES  (Typo  ot  repott  and  Inclusive  dotoo) 

Final  Report  25  June  1968  -  31  March  1969 


5.  AUTHORIS)  (Fieot  nomo,  mtddio  Initiot,  loot  nomo) 

Werner  E.  Schmid 


«.  REPORT  DATE 

March  1969 


•a.  CONTRACT  OR  ORANT  NO. 

N62399-68-C-00M4 

b.  PROJECT  NO. 

Y-F015-21-02-005A 


10.  OISTRIB^m^  STATEMEI^T^l 


7*.  TOTAL  NO.  OF  PACES 


7A.  NO.  OF  REFS 


Sa.  ORIGINATOR'S  REPORT  NUMSERIS) 


OTHER  REPORT  NOISI  (Any  other  nimbnt,  Ihml  msy  bo  asafdnarf 
Ihio  topotl) 


I 

,/V  96.  OTHER  REPORT  6 

thio  ropoit) 

()o\CR  69.030 


Naval  Civil  Engineering  L'aboratory 
Port  Hueheme,  California  93041 


13.  ABSTRACT 


-  The  state  of  the  art  of  predicting  as  well  as  achieving  the 
penetration  of  objects  into  the  sediments  of  the  ocean  floor  is 
presented.  Both  free-fall  as  well  as  driven  penetration  is  reviewed. 
The  information  known  about  ocean  bottom  soils  is  briefly  summarized. 

For  free-fall  penetration  the  impact  velocity,  attitude  and 
direction  are  most  important  and  methods  for  calculating  these 
parameters  are  demonstrated. 

Impact  penetration  relations  are  then  discussed  based  on 
theoretical  as  well  as  empirical  relations. 

The  theoretical  considerations  valid  for  driving  objects  into  the 
bottom  sediments  are  reviewed  and  the  equipment  and  machinery  available 
at  present  (1969->-  for  achieving  penetration  or  for  driving  piles  into 
the  ocean  bottom  are  presented.  Recommendations  for  immediate  use  and 
for  continuing  research  are  made. 


,1473 


S/N  0 lot.  807- 6801 


Unclassified 

Security  Classirication 


u 


nclassified _ 

Security  CUmiftcaHon 


L.4. 


KCV  WOMOt 


penetration,  penetration 

penetrometers 

penetrators 

free  fall  penetration 

driven  penetration 

offshore  pile  driving 

sediment  samplers 

corers 

anchors 

explosive  anchors 
ocean  sediments 
sediment  properties 
sediment  shear  strength 
bearing  capacity 
terminal  velocity 
free  fall  trajectory 
hydrodynamic  drag 


DD  .'.r..l473 

(PAGE  2) 


NOLC 


link  C 


equations 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 


X 

X 

X 


X 

X 

X 

X 


Unclassified 

Security  Classification 


